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ABSTRACT 


We  model  the  problem  of  minimum-cost  procurement  and  allocation  of  anti-ship 
cruise  missiles  to  naval  combat  ships  as  a  two-period  chance-constrained  program  with 
recourse.  Discrete  scenarios  in  two  periods  define  “demands”  for  missiles  (i.e.,  targets 
and  number  of  missiles  required  to  kill  those  targets),  which  must  be  met  with  acceptable 
probabilities.  After  the  first  combat  period,  ships  may  replenish  their  inventories  from  a 
depot,  if  the  depot’s  inventory  suffices.  A  force  commander  assigns  targets  to  ships 
based  on  missile  load-outs  and  target  demands. 

The  deterministic-equivalent  integer  program  solves  too  slowly  for  practical  use. 
We  propose  a  specialized  decomposition  algorithm,  implemented  in  MATLAB,  which 
solves  the  two-period  model  via  a  series  of  single-period  problems.  The  algorithm  yields 
optimal  solutions  for  a  wide  range  of  missile-allocation  directives,  and  usually  near- 
optimal  solutions  otherwise.  We  exploit  the  fact  that  each  single-period  problem  is  a 
probabilistic  integer  program,  whose  solution  must  be  a  p-efficient  point  (PEP)  of  that 
period’s  demand  distribution.  Our  algorithm  uses  PEP-enumeration  techniques 
developed  by  Beraldi  and  Ruszczynski,  and  a  specialized  algorithm  from  Kress,  Penn  and 
Polukarov.  The  algorithm  solves  real-world  problem  instances  in  a  few  minutes  or  less. 
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EXECUTIVE  SUMMARY 


We  present  models  for  the  procurement  and  allocation  of  anti-ship  cruise  missiles, 
a  difficult  problem  facing  many  modern  navies,  which  must  plan  for  many  possible 
combat  scenarios.  The  number  of  missiles  required  by  each  ship  over  an  entire  conflict 
may  exceed  its  capacity,  so  a  ship  may  need  to  return  to  a  depot  at  a  port  to  load  more 
missiles  between  successive  periods  of  combat.  We  consider  two  periods  of  combat, 
“period  I”  and  “period  II,”  and  seek  to  determine  minimum-cost  initial  ship  load-outs 
plus  depot-level  inventory,  while  ensuring,  with  sufficiently  high  probability,  that  all 
ships  can  satisfy  their  assigned  missions  for  each  combat  period.  Each  mission  represents 
a  target  and  its  associated  “demand,”  which  is  the  number  of  missiles  required  to 
successfully  engage  that  target. 

We  formulate  this  problem  as  a  two-period,  three-stage,  chance-constrained 
inventory  model,  which  we  denote  by  CCNIM  (Chance-Constrained  Naval  Inventory 
Model).  We  assume  that  all  targets  will  be  in  range  of  each  combat  ship  in  the  fleet,  and 
the  force  commander  can  assign  any  target  to  any  ship;  however,  the  force  commander 
assigns  at  most  one  target  to  each  ship,  and  does  so  based  on  the  available  missile 
inventories  and  the  target  demands.  We  further  assume  that  ships  have  bounded 
capacities,  carry  no  safety  stocks,  are  not  recalled  to  port  to  offload  excess  inventories 
prior  to  period  II,  and  there  is  no  direct  transshipment  of  missiles  between  ships.  A  “cost 
ratio,”  the  cost  of  a  missile  stored  in  the  depot  for  potential  use  in  period  II  divided  by 
those  initially  allocated  to  the  ships,  reflects  operational  preferences  rather  than  actual 
monetary  costs. 

To  solve  CCNIM,  we  first  formulate  it  as  FFAM  (Fully-Flexible  Assignment 
Model),  a  mixed-integer  program  with  enumerated  scenarios.  Because  we  expect  the 
force  commander  to  follow  some  assignment  heuristic  when  assigning  targets  to  ships, 
we  add  constraints  to  FFAM  that  assign  targets  with  larger  demands  to  ships  with  larger 
inventories.  We  refer  to  this  assignment  policy  as  the  “MAP,”  and  show  that  assigning 
targets  according  to  this  policy  satisfies  any  single-period  scenario  that  can  possibly  be 
satisfied  with  the  existing  ship  inventory  levels.  Consequently,  it  is  the  best  myopic 

xvii 


policy  a  force  commander  can  adopt.  We  also  prove  that  the  inventories  will  be 
sufficient  to  sustain  the  two  periods  of  combat  for  an  arbitrary  assignment  plan  if 
transshipment  is  allowed. 

Clearly,  robustness  can  always  be  achieved  by  procuring  large  quantities  of 
missiles.  However,  if  targets  are  assigned  according  to  the  MAP  in  every  scenario,  then 
the  number  of  missiles  we  must  allocate  to  cover  the  random  demands  with  the  required 
probability  of  success  is  minimal  in  any  particular  combat  period.  Empirically,  we 
observe  that  enforcing  the  MAP  rarely  increases  the  cost  of  an  optimal  solution  in 
FFAM.  We  adopt  CCNIM-MAP,  a  version  of  CCNIM  that  includes  the  MAP,  as  our 
baseline  model. 

Although  enforcing  the  MAP  in  FFAM  improves  solution  times  significantly, 
FFAM  remains  too  difficult  to  solve  for  cases  of  practical  size.  For  example,  an  instance 
of  FFAM,  involving  six  ships  and  five  scenarios  in  each  period,  generates  10,305 
equations  and  14,672  variables  in  the  integer  programming  model.  It  can  be  solved  in 
less  than  an  hour  only  when  the  cost  ratio  is  very  large  or  very  small.  We  therefore 
require  better  solution  methods  to  solve  CCNIM-MAP. 

We  propose  CCNIM-dc,  a  decomposition  algorithm  to  solve  CCNIM-MAP. 
The  solution  is  provably  optimal  for  every  cost  ratio  not  greater  than  1,  and  if  the  cost 
ratio  is  sufficiently  greater  than  1 .  In  other  cases,  the  solution  is  not  provably  optimal, 
but  CCNIM-dc  also  provides  a  lower  bound  on  the  optimal  cost,  so  an  optimality  gap 
can  be  calculated.  In  most  of  the  cases  examined,  the  relative  optimality  gap  is  only  a 
few  percent.  However,  the  optimality  gap  cannot  be  reduced  by  further  computation.  If 
the  gap  is  deemed  too  large,  less  efficient  techniques  must  be  used  to  solve  that  instance. 

We  exploit  the  fact  that  the  problem  of  allocating  sufficient  missiles  to  satisfy  a 
single  period  of  combat,  when  the  MAP  is  enforced,  is  an  instance  of  a  probabilistic 
integer  program  (PIP),  as  described  by  Beraldi  and  Ruszczyhski.  The  feasible  region  of  a 
PIP  has  a  special  structure  that  is  defined  by  a  set  of  “p-efficient  points”  (PEPs),  a 
concept  developed  by  Dentcheva,  Prekopa  and  Ruszczyhski.  CCNIM-dc  enumerates  a 
relatively  small,  cost-ratio-dependent  set  of  period-I  allocations,  and  calculates  the 
minimum  depot  inventory  that  follows  from  each  one.  Minimum  depot  inventories  are 


found  either  by  enumerating  PEPs,  or  by  an  algorithm  proposed  by  Kress  et  al.  for  a 
ground-combat  version  of  CCNIM.  CCNIM-dc  is  not  based  on  linear  or  integer 
programming,  and  can  be  implemented  using  a  standard  programming  language. 

We  implement  CCNIM-dc  in  MATLAB™  version  7.0  on  a  personal  computer 
with  a  2.8  GHz  Intel  Pentium  IV  processor.  Its  solution  times  are  faster  than  those  of 
FFAM  by  several  orders  of  magnitude.  CCNIM-dc  solves  in  a  few  milliseconds 
instances  that  FFAM  cannot  solve  in  an  hour,  and  can  solve  instances  of  practical  size  in 
a  reasonable  length  of  time.  For  example,  cases  with  as  many  as  25  ships  and  25 
scenarios  in  each  period  usually  solve  in  less  than  100  seconds,  and  never  require  more 
than  1,000  seconds. 
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I.  INTRODUCTION 


Two  of  the  key  questions  that  military  logisticians  ask  are:  How  much  ordnance 
should  be  procured,  and  how  should  it  be  distributed  to  combat  units  given  the 
uncertainty  of  combat?  This  dissertation  develops  models  and  solution  methods  to  help 
answer  instances  of  these  questions  in  the  context  of  naval  surface  warfare. 

A.  A  GENERAL  SUPPLY-CHAIN  FRAMEWORK  FOR  COMBAT  UNITS 

Combat  forces  routinely  use  a  variety  of  ordnance  to  suppress  enemy  forces, 
destroy  them,  or  disrupt  their  weapons.  Sustaining  a  combat  force  requires  plans  to 
prescribe  the  total  amount  of  ordnance  of  each  type  that  should  be  procured,  the  initial 
allocation  plan  to  combat  units,  and,  by  implication,  the  amount  of  ordnance  to  be  stored 
in  one  or  more  depots  for  distribution  at  a  later  time.  An  optimal  plan  must  identify  a 
minimum-cost  package  of  ordnance  that  meets  operational  requirements  with  a 
sufficiently  high  probability  of  success.  Uncertainty  arises  because  a  variety  of  combat 
scenarios  may  arise,  each  with  potentially  distinct  demands  for  ordnance. 

We  define  a  combat  scenario,  or  simply  scenario,  as  a  set  of  demands  for 
ordnance  in  a  single  period  of  combat.  A  demand  denotes  the  amount  of  ordnance  that 
must  be  fired  by  a  shooter  (friendly  unit)  to  successfully  prosecute  its  assigned  mission  in 
that  scenario.  A  mission  typically  consists  of  one  or  more  targets  (enemy  units)  that  must 
be  prosecuted.  Because  the  amount  of  ordnance  fired  in  a  successful  engagement  of  a 
specific  target  is  essentially  random  [Ancker  1982,  pp.  1-8],  we  may  represent  a  single 
envisioned  set  of  targets  by  several  different  scenarios.  Note  that  we  only  consider 
demands  for  a  single  type  of  ordnance,  such  as  artillery  shells  or  anti-ship  cruise  missiles 
(ASCMs).  The  duration  of  a  combat  period  is  not  specified  in  days,  but  rather  is  the  time 
between  successive  opportunities  to  replenish  supplies.  In  the  naval  context,  for 
example,  a  combat  period  can  last  several  days,  and  is  delimited  by  “lulls”  in  the  combat, 
during  which  some  or  all  of  the  ships  in  a  fleet  could  return  to  port  for  resupply. 

For  a  force  to  successfully  prosecute  a  scenario,  we  assume  that  all  demands  must 
be  met  using  the  units’  available  supply  of  ordnance.  We  refer  to  successfully  prosecuted 


1 


scenarios  as  satisfied  scenarios.  If  any  unit  cannot  satisfy  the  demand  for  its  assigned 
mission,  the  scenario  is  considered  unsatisfied.  It  is  possible  to  relax  this  assumption  and 
require  that  only  a  fixed  fraction  of  the  demands  be  met  to  satisfy  a  scenario.  We  do  not 
consider  this  model  variant  in  this  dissertation,  however. 

The  planning  horizon  may  include  several  periods  of  combat  in  succession,  and 
the  probability  that  a  particular  scenario  in  period  t>l  occurs  can  be  conditional  on  the 
scenario  that  actually  occurred  in  period  t-1.  We  refer  to  the  sequence  of  single-period 
scenarios  for  t  =  l,...,T  as  a  compound  scenario.  Note  that  in  most  stochastic¬ 
programming  literature,  what  we  call  a  “scenario”  is  referred  to  as  a  “node  in  a  scenario 
tree,”  and  our  “compound  scenario”  is  simply  a  “scenario.”  Our  terminology  is  more 
convenient  when  considering  multi-period  problems  as  a  series  of  single-period 
problems,  a  recurring  idea  in  this  dissertation. 

We  assume  that  the  scenarios,  with  their  respective  probabilities,  have  been 
established  by  planners  as  an  extension  of  the  war  plans.  The  scenarios  can  represent 
different  hypothesized  battles  as  well  as  uncertain  demands  within  each  battle.  The 
duration  of  the  conflict  we  plan  for  is  usually  long  enough  that  combatants  have  the 
opportunity  to  replenish  their  supplies,  often  more  than  once  [Rabinovitch  1997,  p.  252]. 

Each  combat  unit  has  some  inventory  capacity  associated  with  it.  The  capacity 
may  be  “strict,”  as  is  the  case  for  missiles  carried  aboard  combat  ships,  or  “elastic,”  as  is 
the  case  for  many  ground  units,  which  can  increase  the  number  of  supply  vehicles 
traveling  with  the  combatants.  A  unit’s  total  requirement  for  ordnance  to  satisfy  its 
forecasted  assigned  missions  in  a  compound  scenario  may  exceed  its  carrying  capacity. 
In  that  case,  the  unit  may  be  forced  to  replenish  its  supply  after  some  period  of  combat. 

Even  if  a  unit’s  carrying  capacity  is  sufficient  to  meet  the  largest  conceivable 
compound  demand,  it  may  be  preferable  not  to  load  the  unit  to  that  level.  Eollowing  the 
realization  of  a  scenario  in  some  period,  the  probabilities  of  scenarios  in  subsequent 
period  are  updated.  It  is  possible  then  that  high-demand  scenarios  no  longer  seem  likely, 
and  the  requirements  for  certain  units  can  be  lowered.  We  wish  to  avoid  situations  in 
which,  after  a  period  of  combat,  some  units  hold  excess  ordnance  that  may  be  needed  by 
other  units  in  the  next  period.  Replenishment  from  central  locations  is  advantageous 
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then,  since  it  utilizes  risk  pooling  (e.g.,  Simchi-Levi  et  al.  [2000,  pp.  56-60])  to  reduce 
the  total  amount  of  ordnance  that  must  be  deployed  in  the  first  combat  period.  Note  that 
risk  pooling  refers  to  the  risk  of  shortages  due  to  high  demands,  and  not  to  the  risk  of 
losing  ordnance  associated  with  disabled  units.  (It  is  clear  that  reducing  the  amount  of 
deployed  ordnance  reduces  the  amount  of  ordnance  that  can  be  lost  to  enemy  actions,  but 
we  do  not  address  such  issues  in  this  dissertation.) 

One  of  the  common  features  of  our  problem  and  stochastic  supply-chain  models 
in  the  civilian  sector  is  recourse.  After  a  period  of  combat,  we  can  resupply  ships  in 
order  to  improve  our  chances  for  success  in  subsequent  periods.  This  action  constitutes 
recourse,  as  defined  in  the  stochastic -programming  literature  (e.g.,  Birge  and  Louveaux 
[1997,  pp.  84-100,  122-127]).  Thus,  the  model  we  eventually  create  will  be  a  multi-stage 
stochastic-programming  model  with  recourse. 

B.  CHANCE-CONSTRAINED  PROGRAMMING 

One  of  the  key  features  that  distinguishes  the  problem  of  procuring  and  deploying 
ordnance  from  other  supply  problems  is  the  singularity  of  war  [Kress  2002,  p.  242].  This 
contrasts  with  civilian  supply-chain  planning  models  that  typically  apply  to  repetitive 
scenarios  over  relatively  long  time  horizons  [Diwekar  2002].  Because  we  expect  the 
planned  supply  chain  to  be  tested  in  war  only  once,  albeit  across  multiple  time  periods, 
our  models  will  incorporate  probability  requirements,  i.e.,  probabilistic  constraints,  or 
chance  constraints.  (For  example,  see  Birge  and  Louveaux  [1997,  pp.  103-108].)  Multi¬ 
stage  supply-chain  models  using  chance  constraints  to  guarantee  service  levels  do  exist 
(e.g.,  Chames,  Cooper  and  Symonds  [1958],  Murr  and  Prekopa  [1996],  Bassok  et  al. 
[2002],  Paschalidis  et  al.  [2004]),  but  a  typical  multi-stage  stochastic  program  would 
evaluate  the  effect  of  uncertainty  through  an  objective  function  involving  expected  costs 
and/or  penalties  (see  Porteus  [1990]  and  the  references  therein)  as  opposed  to  using 
chance  constraints. 

The  first  mathematical-programming  model  to  incorporate  chance  constraints  is 
attributed  to  Charnes  et  al.  [1958],  who  model  the  multi-period  planning  of  heating-oil 
production  under  uncertain  demand.  They  introduce  multiple  individual  chance 
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constraints  to  guarantee  each  period’s  demand  is  met  with  the  required  probability. 
Miller  and  Wagner  [1965]  introduce  joint  chance  constraints  as  a  theoretical  extension  of 
the  work  of  Chames  et  al.  Since  the  1960s,  the  theoretical  understanding  of  chance 
constraints  has  increased  significantly  (see  Prekopa  [1995]  and  the  references  therein), 
and  several  specialized  algorithms  for  solving  problems  formulated  with  such  constraints 
have  been  proposed.  However,  most  of  this  research  has  been  restricted  to  models 
involving  continuous  distributions  and  decision  variables 

When  the  probability  distributions  are  discrete,  or  the  decision  variables  integer, 
chance-constrained  models  become  significantly  harder  to  solve  [Beraldi  and 
Ruszczyhski  2001],  and  there  is  relatively  little  literature  dedicated  to  such  models. 
Dentcheva,  Prekopa  and  Ruszczyhski  [2000]  provide  an  extensive  analysis  of  the 
properties  of  the  feasible  region  of  the  discrete  probabilistic  linear  program  (DPLP)  with 
a  joint  chance  constraint  involving  the  discrete  random  vector  ^ : 


(DPLP)  min  c^x  (1.1) 

s.t. 

Pr{rx>^}>p  (1.2) 

Ax>b  (1.3) 

x>0  (1.4) 


To  analyze  this  problem,  the  authors  use  the  concept  of  a  p-efficient  point,  first 
introduced  by  Prekopa  [1990]: 

Definition:  Let  pe[0,l],  and  let  F  denote  the  probability  distribution  for  the  n- 
dimensional,  discrete  random  vector  % .  A  point  v  e  M"  is  called  a  p-efficient  point  (PEP) 
of  F  ,  if  F(v)  >  p  and  there  is  no  y  <  v ,  y  such  that  F  (y)  >  p  .  | 

Dentcheva  et  al.  prove  that  there  is  a  positive  and  finite  number  of  PEPs 
associated  with  any  discrete  random  vector,  and  that  the  feasible  region  for  the 
probabilistic  constraint  in  any  DPEP  can  be  reformulated  as  a  union  of  cones  emanating 
from  those  PEPs.  In  particular,  let  v^  ,  je  J  ,  be  all  the  PEPs  of  F  and  let  K.  =  v^  -l-M" 

denote  the  cone  emanating  from  v^ ;  then  constraint  (1.2)  is  equivalent  to 
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(1.5) 


j<^J 

We  can  solve  a  DPLP  by  enumerating  all  its  PEPs,  and  by  then  processing  the 
simple  linear  programs  that  result  from  considering  each  cone  separately.  Prekopa, 
Vizvari  and  Badics  [1996]  and  Beraldi  and  Ruszczynski  [2001]  develop  algorithms  for 
enumerating  PEPs;  Beraldi  and  Ruszczynski  [2002]  provide  additional  details.  If  %  has 
many  dimensions,  the  algorithms  may  be  computationally  expensive,  so  the  authors 
propose  iterative  methods  that  generate  new  PEPs  as  needed,  en  route  to  finding  an 
optimal  solution. 

Replacing  constraints  (1.4)  with  the  integrality  restriction  xeZ"  creates  a 

variation  of  a  DPLP  that  Beraldi  and  Ruszczynski  [2001]  call  a  probabilistic  integer 
program  (PIP).  These  authors  use  the  PEPs  of  E  in  a  specialized  branch-and-bound 
procedure  to  solve  the  PIP.  Instead  of  generating  PEPs  one  at  a  time,  as  Dentcheva  et  al. 
do,  they  begin  with  an  upper  bounding  vector  on  the  PEPs  as  the  root  of  a  PEP  search 
tree.  They  define  the  level  of  an  integer  vector  as  the  sum  of  its  components,  and  use  a 
backwards  enumeration  procedure  to  generate  PEP  candidates,  i.e.,  nodes  in  the  search 
tree,  one  level  at  a  time.  Eor  each  candidate,  they  generate  an  integer  program  (IP)  in 
which  the  probabilistic  constraint  is  replaced  by  some  lower  bound  on  the  value  of  any 
future  node  emanating  from  that  branch.  They  then  solve  a  linear  relaxation  of  that  IP, 
and  are  able  to  prune  the  search  tree  if  the  objective  value  is  higher  than  that  of  the  best 
known  feasible  solution.  Leaves  in  the  next  level  are  generated  for  branches  that  are  not 
pruned.  This  approach  works  better  than  that  proposed  by  Dentcheva  et  al.  for  integer 
decisions  variables. 

In  another  approach  to  solving  PIPs,  Ruszczynski  [2002]  reformulates  the 

probabilistic  constraints  as  set-covering  constraints,  in  a  model  that  chooses  which  of  the 

discrete  scenarios  are  to  be  satisfied  and  which  are  not.  A  partial  ordering  is 

defined  on  the  scenarios,  based  on  the  component-wise  values  of  their  respective  right- 

hand-side  realizations:  s<s'  ^  .  If  the  ordering  is  consistent,  we  can  generate 

constraints  of  the  form  ,  where  z^  =  I  if  scenario  5  is  not  satisfied  and  0  otherwise. 

The  PIP  can  be  viewed  as  a  knapsack  problem,  where  we  choose  the  scenarios  that  are 
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not  satisfied  subject  to  a  maximum  aggregate  probability  threshold.  Ruszczynski 
proposes  a  solution  method  that  iteratively  generates  and  lifts  cover  cuts  until  a  solution 
to  the  PIP  is  found. 

C.  INITIAL  FORMULATIONS 

This  dissertation  presents  work  that  extends  research  initiated  in  Kress,  Penn  and 
Polukarov  [2004],  hereafter  referred  to  as  “KPP.”  These  authors  formulate  a  two-level, 
two-period,  chance-constrained  inventory  model  (CCIM)  with  recourse,  in  a  military 
setting.  The  KPP  version  of  CCIM,  which  we  refer  to  as  CCGIM  (“G”  stands  for 
“ground”),  captures  the  basics  of  ground  combat,  in  which  units  are  assigned  different 
sectors  of  a  battlefront,  and  each  must  meet  the  opposition  in  its  sector.  We  refer  to  the 
demands  in  each  scenario  as  being  rigidly  assigned,  because  once  a  scenario  unfolds, 
each  combat  unit  is  tasked  with  a  specific  mission,  in  its  assigned  sector,  and  the  demand 
for  ordnance  induced  by  that  mission.  Each  combat  unit  carries  its  own  ordnance  and 
attempts  to  prosecute  its  mission  using  that  ordnance.  Additional  ordnance  is  stored  at  a 
central  depot  that  can  ship  it  to  units  whose  supplies  need  to  be  replenished  for  the  second 
period  of  combat. 

KPP  assume  that  the  units  must  rely  on  their  initial  stocks  to  satisfy  any  demands 
in  the  first  period  (period  I),  and  that  each  unit  maintains  sufficient  safety  stock  with 
which  to  carry  out  any  mission  even  when  its  nominal  supply  cannot.  However,  if  any 
unit  must  use  its  safety  stock,  the  scenario  is  considered  unsatisfied.  Following  period  I, 
units  replenish  their  supplies  (including  their  safety  stocks,  if  used)  to  acceptable  levels  in 
view  of  the  second  period’s  (period  II)  projected  scenarios.  Replenishments  are  made 
from  the  central  depot,  or  by  inter- unit  transfers.  Thus,  excess  ordnance  can  be 
redistributed  following  period  I.  KPP  assume  the  cost  of  capacity  at  the  combat  units  is 
linear  in  the  amount  of  supply,  and  this  cost  is  factored  into  the  price  of  ordnance  initially 
allocated  to  the  units.  Hence,  it  is  reasonable  to  assume  that  the  cost  of  carrying  ordnance 
with  the  ground  forces  is  greater  than  that  of  storing  it  in  the  depot. 

Under  the  conditions  described  above  (possible  consumption  from  the  safety 
stock,  possible  transshipment  during  resupply,  and  greater  cost  for  ordnance  at  the 
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combat  units  than  at  the  depot),  CCGIM  can  be  decomposed  into  two  separate  single¬ 
period  problems  and  solved  sequentially.  The  period-I  problem  allocates  a  minimum 
amount  of  ordnance  so  as  to  be  able  to  meet  uncertain  demands  in  period  I  with  a 
specified  probability.  The  period-I  allocations  are  used  as  parameters  in  the  period-II 
problem  that  minimizes  the  amount  of  ordnance  to  be  added  to  the  units’  inventories  from 
the  depot.  Each  of  those  problems  is  a  special  case  of  a  combinatorial  problem  that  KPP 
call  the  Minmax  Subset  Problem  (MSP).  KPP  provide  an  empirically-efficient,  exact 
algorithm  to  solve  the  MSP  (see  Kress  et  al.  [2004]),  and  prove  that  solving  the  two  MSP 
problems  sequentially  solves  CCGIM. 

We  present  CCNIM  (“N”  stands  for  “Naval”),  a  second  variant  of  CCIM  that 
models  consumption  of  ASCMs  in  naval  surface  warfare.  ASCMs  constitute  the  major 
weapon  system  of  modern  navies  that  do  not  rely  on  airpower  (for  example,  the  fleets  of 
Denmark,  Greece,  and  the  Netherlands  [Baker  2002]).  The  naval-combat  inventory 
model  is  more  complex  than  CCGIM,  however,  and  cannot  be  solved  with  KPP’s  simple 
decomposition  procedure. 

A  key  difference  between  ground  combat  and  naval  surface  combat  is  the 
flexibility  in  mission  assignments  within  a  specific  scenario.  A  high  degree  of 
interchangeability  arises  in  meeting  demands  because  ASCM  ranges  are  often  long 
compared  with  the  distances  between  battling  combat  ships  [Hughes  1995]. 
Consequently,  once  a  set  of  targets  becomes  evident,  the  combat  force  has  significant 
freedom  in  assigning  targets  to  shooters.  Targets  are  assigned  among  the  available 
shooters  based  on  available  supplies  and  tactical  positions.  In  some  cases,  each  target  is 
within  range  of  all  potential  shooters;  we  refer  to  this  situation  diS  fully  flexible.  When 
the  tactical  situation  does  not  allow  every  ship  to  engage  any  target,  we  refer  to  the 
targets  in  that  scenario  as  being  semi-flexibly  assigned.  In  any  case,  we  assume  each 
target  is  assigned  to  one  specific  shooter,  and  its  demand  is  not  shared.  This  dissertation 
focuses  on  the  fully-flexible  case,  which  provides  a  lower  bound  (optimistic  value)  on  the 
required  number  of  missiles  for  any  other  case. 

Other  conditions  necessary  for  KPP’s  simple  decomposition  approach, 
appropriate  for  ground  combat,  do  not  hold  for  naval  combat.  We  assume  the  existence 


7 


of  lower  and  upper  bounds  on  the  number  of  missiles  that  will  be  carried  aboard  a  ship 
because  of  fixed  capacity  constraints,  operational  considerations  and  doctrine.  The 
assumption  of  a  safety  stock  is  unreasonable  for  combat  ships,  which  typically  have 
onboard  inventories  of  at  most  eight  ASCMs  [Baker  2002].  Therefore,  each  ship  can 
expend  no  more  than  its  onboard  inventory,  and  “backorders”  on  missiles  do  not  arise. 
We  assume  that  if  the  required  number  of  missiles  to  cover  an  assigned  mission  exceeds 
the  number  available  onboard  the  ship,  the  ship  will  “stay  and  fight”  and  expend  its  entire 
inventory.  However,  scenarios  in  which  any  ship  cannot  satisfy  its  assigned  mission  are 
considered  unsatisfied.  A  small  adjustment,  which  we  do  not  pursue  here,  can  model 
situations  in  which  an  insufficiently  armed  shooter  avoids  engaging  its  assigned  target. 

Once  a  missile  is  procured,  we  assume  that  the  monetary  costs  of  loading  it  onto  a 
ship  prior  to  period  I  or  storing  it  at  a  shore -based  depot  are  essentially  equivalent,  and 
not  considered  in  the  decision  process.  However,  we  still  define  a  cost  ratio  between  the 
cost  of  storing  a  missile  in  the  depot  and  the  cost  of  initially  allocating  a  missile  to  a 
combat  ship.  This  ratio  represents  our  operational  preferences  regarding  allocation 
strategies,  by  defining  the  relative  value  of  missiles  initially  allocated  to  ships  compared 
with  those  stored  in  the  depot  for  period  II.  Essentially,  whenever  the  ratio  is  not  1,  we 
are  willing  to  buy  more  than  the  minimum  total  number  required  in  order  to  obtain  an 
allocation  plan  that  suits  us.  If  we  wish  to  reduce  the  operational  burden  of  carrying  out 
wartime  replenishment  operations  in  port,  we  set  the  depot  cost  higher  than  the  on-ship 
cost,  so  the  ratio  exceeds  one.  If  we  wish  to  avoid  the  risk  of  losing  missiles  due  to  own- 
force  casualties  or  expenditure  on  low-value  targets,  we  can  set  the  ratio  below  1.  For 
example,  setting  the  cost  ratio  to  1.25  means  that  we  view  a  solution  that  procures  five 
missiles  and  allocates  them  to  the  ships,  to  be  equally  desirable  as  a  solution  that  requires 
us  to  procure  only  four  missiles  that  must  be  stored  at  the  depot.  Similarly,  setting  the 
cost  ratio  to  l  +  £  ,  for  a  small  £■>0,  corresponds  to  a  solution  that  acquires  the 
minimum  total  number  of  missiles,  but  allocates  as  many  of  them  as  possible  to  the  ships. 

Combat  ships  cannot  normally  transfer  missiles  directly  between  themselves,  as 

ground-combat  units  can.  That  is,  no  direct  transshipment  occurs.  We  also  assume  that 

no  ship  will  be  required  to  make  a  port  call  for  the  sole  purpose  of  offloading  missiles  for 

the  use  of  some  other  ship.  Therefore,  the  model  assumes  that  no  transshipment 
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capability  exists  whatsoever,  and  any  requirement  for  missiles  that  a  ship  might  have 
following  the  first  period  of  combat  must  be  met  by  missiles  stored  in  onshore  depots. 
This  assumption  is  conservative  because,  in  practice,  some  ships  will  return  to  port  for 
other  reasons,  and  will  be  able  to  offload  missiles.  Because  of  the  flexible  assignments, 
lack  of  safety  stocks,  no  transshipment,  bounded  inventories,  and  arbitrary  cost  ratio,  we 
see  that  new  techniques  are  required  to  solve  a  CCIM  in  the  context  of  naval  combat. 

D.  DISSERTATION  OUTLINE 

This  dissertation  presents  two-period,  chance-constrained  models  for  the 
procurement  and  allocation  of  ASCMs  for  naval  surface  warfare.  Chapter  II  develops 
CCNIM  and  its  single -period  variant  as  stochastic  programs.  It  also  presents  a 
deterministic  equivalent  IP  of  CCNIM,  which  can  be  implemented  in  standard 
optimization  software.  Because  the  IP’s  solution  times  are  often  inadequate  for  practical 
purposes,  the  rest  of  this  dissertation  develops  faster  solution  methods.  Chapter  III 
describes  the  properties  of  the  “Monotonic  Assignment  Policy,”  by  which  a  commander 
assigns  targets  with  larger  demands  to  ships  with  larger  inventories.  Because  of  these 
properties,  we  argue  that  CCNIM  should  be  modified  to  include  such  a  policy,  and 
obtain  CCNIM-MAP  (CCNIM  with  constraints  enforcing  the  Monotonic  Assignment 
Policy).  Chapter  IV  explores  the  feasible  region  of  CCNIM-MAP,  and  identifies  points 
within  the  feasible  region  that  will  be  optimal  if  specific  conditions  on  the  objective 
function  hold.  In  Chapter  V,  we  use  these  observations  to  develop  CCNIM-dc,  a 
specialized  decomposition  algorithm,  which  can  be  implemented  in  a  standard  computer¬ 
programming  language.  CCNIM-dc  provides  an  optimal  solution  for  a  wide  range  of 
cost  ratios,  and  bounds  the  optimal  cost  otherwise.  It  requires  several  orders  of 
magnitude  less  time  than  the  IP  to  solve  problems  of  practical  size. 
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II.  MATHEMATICAL  MODELS 


This  chapter  describes  a  chance-constrained  inventory  model  for  naval  surface 
warfare,  CCNIM.  Initially,  we  develop  the  single-period  model  CCNIM-sp,  and  show 
how  it  relates  to  the  theory  of  probabilistic  integer  programs  (PIPs).  We  then  develop  the 
full  model,  CCNIM,  which  spans  two  periods  of  combat.  CCNIM  allows  full  flexibility 
in  mission  assignment,  i.e.,  any  ship  may  prosecute  any  mission,  although  shooter-to- 
mission  assignments  are  one-to-one.  It  assumes  no  transshipment  capability,  direct  or 
indirect,  between  the  ships,  and  explicitly  accounts  for  the  limited  onboard  inventories. 
We  then  develop  the  reduced  fully-flexible  assignment  model  (RFFAM),  a  deterministic 
equivalent  integer-programming  model  of  CCNIM  with  improved  performance 
compared  to  FFAM,  developed  in  Avital  [2004].  RFFAM  is  substantially  faster  than 
FFAM,  and  solves  some  practical  instances  in  a  reasonable  length  of  time.  However,  for 
other  instances,  RFFAM  may  fail  to  produce  even  feasible  solutions  in  one  hour  of 
computation.  For  this  reason,  we  dedicate  the  later  chapters  of  this  dissertation  to  the 
development  of  specialized  solution  techniques  for  CCNIM. 

A.  CHANCE-CONSTRAINED  MODELS  FOR  NAVAL  WARFARE 

I.  CCNIM  for  a  Single  Period 

We  first  develop  a  single-period  model  of  the  inventory  problem  we  consider.  Let 
n  denote  the  number  of  combat  ships  (shooters)  on  the  planner’s  side.  Let  i'e  S  denote 
an  index  set  of  combat  scenarios,  such  that  each  scenario  has  probability  of  occurrence 
(p,>Q  and  ^  =  1 .  Each  scenario  comprises  a  vector  of  n  demands  associated  with 

.V€S 

n  missions,  which  the  shooters  should  satisfy  in  scenario  s.  We  let  d  denote  the  n- 
dimensional,  integer,  random  vector  of  demands  drawn  from  the  discrete  set  of  scenarios 

In  practice,  the  vector  of  n  mission  demands  is  constructed  from  a  conceptualized 
set  T  of  potential  targets.  We  associate  a  specific  demand  with  target  t  in  scenario 
5 ,  which  represents  the  number  of  missiles  required  to  successfully  engage  the 
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corresponding  target  in  that  scenario.  Because  that  number  is  generally  stochastic,  we 
may  generate  several  scenarios  with  different  mission  demands  to  represent  the  same 
target  set.  We  set  =  0  if  target  t  does  not  actually  appear  in  scenario  5  (it  may  not 

even  be  present).  Let  denote  the  set  of  engaged  targets  in  scenario  5,  defined  as 
=  {t  G  T I  >  0}  ,  and  assume  |tJ  <  n  for  all  s. 

Without  loss  of  generality,  we  construct  by  arranging  components  in  such 
a  way  that  d,^>0  for  t  =  l,...,|Tj  and  d,^=Q  for  t  =  |Tj  +  l,...,max|n,  |t||  .  Then, 

dj  =  (djj,...,d„J .  Because  each  mission  comprises  at  most  one  target,  we  often  use  the 
term  “target”  in  place  of  “mission”  in  the  following  text. 

For  example,  suppose  that  there  are  three  shooters,  and  that  |t|  =  4 .  We  assume 
that  in  no  scenario  would  the  three  shooters  be  required  to  engage  all  of  the  enemy  targets 
concurrently.  In  scenario  ,  we  encounter  two  of  the  targets,  and  require  2  missiles  to 

successfully  engage  each  one,  and  thus  dj  =  (2,2,0) .  In  scenario  s^,  we  encounter  the 
same  targets,  but  require  two  salvos  to  successfully  engage  the  second  target,  setting 
d,=  (2,4,0). 

We  minimize  the  procurement  cost  of  missiles  required  to  satisfy  demands  in 
projected  scenarios,  while  satisfying  a  user-specified  minimum  probability  of  success  p  . 
A  scenario  is  considered  satisfied  if  all  ships  have  enough  missiles  to  satisfy  the  demands 
associated  with  their  assigned  missions  in  that  scenario.  The  model  chooses  an  allocation 
plan  V  that  allocates  y  missiles  to  ship  i  and  an  assignment  plan  that  assigns  one 
mission  to  each  ship  in  each  scenario.  Ship  i’s  pre-combat  inventory  v^  is  maintained 
between  a  discretionary  lower  bound  y_-  and  a  physical  upper  bound  v- . 

Because  mission  assignments  are  fully  flexible,  two  scenarios  are  essentially 

equivalent  if  one  scenario’s  demand  vector  is  a  permutation  of  the  other’s.  For  the  sake 

of  solution  efficiency,  if  equivalent  scenarios  are  presented  as  input  data,  they  should  be 

consolidated  and  appropriate  adjustments  made  to  the  probability  data.  (Or,  if  an 

automatic  scenario  generator  is  used  to  create  scenarios,  it  should  be  adjusted  to  avoid 
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producing  equivalent  scenarios.)  We  also  assume  that  each  scenario  is  potentially 
feasible,  meaning  that  the  ships  can  carry  the  required  inventories  to  satisfy  that  demand 
vector.  If  some  hypothesized  scenario  cannot  be  satisfied  because  some  of  its  demands 
exceed  the  capacity  upper  bounds,  it  should  be  excluded  from  the  set  of  scenarios,  and 
the  probability  threshold  requirement  adjusted  accordingly. 

The  flexibility  in  mission  assignments  is  expressed  through  the  nxn  permutation 
matrix  T(v,d) ,  belonging  to  the  set  Y  of  all  binary  matrices  with  a  single  1  in  each  row 

and  each  column.  In  our  formulation,  I^^(v,d)  =  l  if  we  assign  mission  m  to  ship  i. 

Because  targets  are  assigned  flexibly,  we  may  choose  an  assignment  plan  after  the 
demands  are  observed  and  the  existing  allocations  of  missiles  are  considered.  Therefore, 
the  single-period  problem  is  a  two-stage  stochastic  program  with  recourse,  in  which  the 
ship  inventories  v  constitute  first-stage  decision  variables,  and  the  elements  of  T(v,d) 

constitute  second-stage  recourse  variables.  In  practice,  target  assignments  may  not  be 
one-to-one.  However,  sharing  demands  typically  reduces  the  effectiveness  of  a  missile 
salvo,  and  handling  more  than  a  single  target  by  some  ship  may  increase  the  time  it 
requires  to  plan  the  engagement.  Consequently,  we  require  sufficient  inventories  to 
satisfy  the  scenarios  under  the  one-to-one  assumption. 

There  are  several  ways  to  explicitly  formulate  the  optimal  assignment  plan 
T(v,d).  The  expression  for  T(v,d)  in  (2.4)  minimizes  the  maximum  single-target 

shortage,  and  the  demand  vector  d  is  satisfied  if  T(v,d)  can  be  found  which  reduces 
that  shortage  to  zero.  The  Chance-Constrained  Naval  (Surface  Warfare)  Inventory  Model 


for  a  Single  Period  (CCNIM-sp)  can  be  expressed  as: 

(CCNIM-sp)  minify  (2.1) 

V 

S.t. 

Prjv  >  T(v,d)d|  >  p  (2.2) 

v<v<v  (2.3) 

where  T(v,d)G  argmin|max|(Td). -V;|  |  (2.4) 
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For  future  reference,  note  that  if  we  could  not  explicitly  express  the  form  of  the 
optimal  permutation  F  (v,d) ,  we  would  treat  that  functional  form  as  a  decision  variable, 

and  write  CCNIM-sp  as:  minify ,  s.t.  constraints  (2.2)  and  (2.3). 

v.r(-) 

Because  F(v,d)  represents  “recourse”  for  the  shooters,  CCNIM-sp  is  not  a  PIP, 

as  defined  by  Beraldi  and  Ruszczynski  [2001]  (see  section  A  of  Chapter  I),  except  for 
some  important  special  cases.  In  general,  the  concept  of  a  PEP  is  not  well  defined  for 
this  problem  because  the  optimal  permutation  of  d  depends  on  the  number  of  missiles 
allocated  to  each  ship  and  cannot  be  known  a  priori.  An  alternative  formulation,  where 
the  permutations  are  applied  to  the  vector  of  allocations  v  instead  of  to  d ,  is  also  not  a 
PIP,  in  general.  Although  PEPs  can  now  be  found  for  the  distribution  of  d ,  the  optimal 
solution  need  not  lie  in  the  union  of  cones  emanating  from  the  PEPs,  a  characteristic  of  a 
PIP.  Eor  example,  suppose  that  the  six  permutations  of  the  vector  (3,2,1)  are  the  possible 
scenarios,  with  equal  probability.  If  we  set  p  =  5/6,  the  only  PEP  is  (3,3,3),  but  the 
optimal  allocation  plan  is  clearly  (3,2,1),  because  we  can  assign  the  first  ship  to  the  target 
with  demand  1,  the  second  ship  to  the  target  with  demand  2  and  the  third  ship  to  the 
target  with  demand  3  in  every  scenario. 

There  are  some  important  special  cases  where  PEPs  do  apply.  Under  some 
assignment  policies  (perhaps  coupled  with  other  constraints),  the  matrix  T(v,d)  depends 

only  on  the  problem  data,  i.e.,  T(v,d J  =  T(d^ )  .  In  that  case,  we  can  define 

equivalent  demand  vectors  and  eliminate  the  matrix  Y  altogether,  obtaining  a  PIP.  Eor 
example,  in  CCGIM,  where  targets  are  rigidly  assigned,  we  have  Y  =  I ,  and  no 
manipulation  of  the  demand  vectors  is  needed.  If  Y  is  eliminated.  The  PEPs  of  the 
distribution  of  the  demand  vectors  are  useful,  and  can  be  found,  for  example,  by  the 
enumeration  schemes  described  by  Beraldi  and  Ruszczynski  [2002] . 

2.  CCNIM 

In  the  actual  problem  we  wish  to  solve,  we  assume  two-periods  of  combat, 
between  which  ships  may  return  to  port  to  replenish  their  supply  from  a  central  depot. 
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The  index  set  of  period-I  scenarios,  s\  is  defined  as  in  CCNIM-sp.  Let  s'g  s“  denote 
the  index  set  of  period- II  combat  scenarios.  Each  scenario  s'e  s“  has  an  associated 
conditional  probability  ,  such  that  ^  =1  V5  e  Therefore,  the  demand  vector 

i'es” 

for  period  II,  d“,  is  an  integer  vector  whose  distribution  may  depend  on  the  scenario 
occurring  in  period  1.  We  denote  a  realization  of  this  conditional  random  vector  by  d“ . 
We  refer  to  a  realization  of  (d’,d”)  as  a  “compound  scenario”,  as  defined  in  Chapter  1. 
For  each  period-I  scenario  5  e  s' ,  we  define  the  subset  s“  (5)  =  s“  I  (p^y  >  0| ,  which 
indexes  all  possible  demand  vectors  that  might  follow  in  period  11.  The  subsets  s“  (5') 
may  not  be  disjoint  or  even  different  for  the  various  5  e 

Let  v’  denote  the  total  number  of  missiles  initially  allocated  to  the  ships,  and  let 
denote  the  number  of  missiles  initially  allocated  to  the  depot  for  potential  use  in 
period- II.  Let  v/  and  vf  denote  the  number  of  missiles  allocated  to  ship  i  in  period  I  and 
in  period  II,  respectively.  We  assume  that  the  cost  of  allocating  missiles  to  the 
shooters  may  be  different  from  the  cost  of  storing  missiles  at  the  depot.  Initial 
allocation  and  assignment  plans  are  chosen  in  period  I  so  that  the  period-I  scenarios  are 
satisfied  with  a  user-specified  probability  .  Following  period  I,  CCNIM  calculates 
each  ship’s  remaining  inventory  .  Because  a  ship  may  not  expend  more  missiles  than  it 

carries,  if  its  assigned  demand  in  period  I  exceeds  its  inventory,  then  the  remaining 
inventory  is  zero.  Then,  for  each  ship  i,  the  model  supplements  its  remaining  inventory 
with  w,  >  0  missiles  drawn  from  the  depot,  thereby  setting  vf .  Missions  in  period  II  are 

assigned  so  that  period-II  scenarios  are  satisfied  with  a  user-specified  probability  pf , 
which  may  depend  on  the  realized  demand  scenario  in  period-I.  CCNIM  also  ensures 
that  each  ship’s  inventory  is  maintained  within  required  limits  at  the  outset  of  each 
combat  period.  We  illustrate  the  “flow”  of  missiles  in  Figure  1. 
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Figure  1.  Schematic  Missile  Flow  in  a  Compound  Scenario  for  Three  Ships. 
We  initially  procure  +  missiles  at  a  cost  of  c^x^  .  We  allocate 
v/  missiles  to  ship  i ,  and  store  in  the  depot.  The  ships  face  a  random 
set  of  targets  with  a  known  distribution,  which  are  represented  by  the 
demand  vector  d’ .  The  commander  assigns  one  demand  to  each  ship, 

represented  by  the  permutation  matrix  7’ ,  and,  following  the  engagement, 
r.  missiles  remain  on  ship  i.  We  supplement  ship  /’s  inventory  by  w. 

missiles,  for  a  total  of  vf .  In  period  II,  the  ships  face  the  random  demands 
represented  by  d“ ,  and  the  commander  again  may  choose  which  demand 
to  assign  to  each  ship,  represented  by  T" . 


Below,  we  formulate  CCNIM  as  a  chance-constrained  program,  which  may  be 
viewed  as  a  three-stage  stochastic  program.  The  first  stage  determines  the  total  number 
of  missiles  to  procure  and  how  many  to  allocate  to  each  ship  in  period  I.  Once  the 
period-I  scenario  is  revealed,  the  model’s  second  stage  assigns  missions,  and,  following 
the  battle  in  period  I,  the  number  of  missiles  to  add  to  each  ship’s  inventory  for  period  II. 
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We  denote  optimal  second-stage  decisions  by  F’ (v,d,Fjj  (*1(1))  and  w 
respectively,  where  Fjj(*ld)  denotes  the  conditional  distribution  of  the  period- II 
demands.  After  the  period  II  scenario  is  revealed,  the  model’s  third  stage  assigns 
missions  for  that  period;  we  denote  the  optimal  third  stage  decisions  by  F”(v,d). 

However,  the  only  “real”  decision  variables  are  the  first-stage  variables,  after  which  the 
system  is  set  in  motion,  and  we  simply  observe  whether  the  random  demands  can  be 
covered  with  sufficiently  high  probability.  The  optimal  second-stage  and  third-stage 
decisions  are  any  combination  of  assignment  plans  and  replenishment  vectors  that  allow 
us  to  procure  the  minimum-cost  package.  For  the  reader’s  convenience,  we  summarize 
the  notation  used  in  CCNIM,  followed  by  the  model’s  statement: 


Parameters  [units] 

d’,  d”  random  demand  vector  associated  with  period  I  and  period  II, 
respectively  [missiles] 

Fjj  ('Id)  conditional  distribution  of  period- II  demands 

p"  probability  threshold  for  period  I,  and  for  period  II,  if  scenario  s  occurs 
in  period  I 

Cj,  Cj  unit  cost  of  procuring  a  missile  and  allocating  it  to  a  ship  or  to  the 

depot,  respectively  [$/missile] 

y,  V  discretionary  lower  bound  and  physical  upper  bound  on  the  number  of 

missiles  that  may  be  allocated  to  the  ships  [missiles] 


Decision  Variables  [units] 

x\  total  number  of  missiles  allocated  initially  to  the  ships  and  to 

the  central  depot,  respectively  [missiles] 

v',  v“  numbers  of  missiles  allocated  to  the  ships  at  the  outset  of 

period  I  and  period  II,  respectively  [missiles] 

r  numbers  of  missiles  left  on  the  ships  following  period-I 

[missiles] 

w  (r,  F„  ('Id))  numbers  of  missiles  receives  from  the  depot  following  period-I 

[missiles] 
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y'  (v,d,Fjj  ('Id))  assignment  matrix  of  demands  to  ships 

II ,  respectively 

F"(v,d) 

in  period  I  and  period 

(CCNIM)  min  c.v'+Cjv" 

(2.5) 

S.t. 

Pr{v'>F'(v',d',Fi,(.|d‘))d‘}>p' 

(2.6) 

1 

II 

O 

(2.7) 

y<  v’  <  v 

(2.8) 

v'  e  Z; 

(2.9) 

where  r  =|v^  -F’  (v’,d’,Fj;  (‘Id'ljd' j 

(2.10) 

v"=r  +  w(r,F4.|d')) 

(2.11) 

and  where  w(»)e  Z"  is  chosen  such  that 

Pr{v">F"(v",d")d“}>pf 

(2.12) 

o 

VI 

1 

• 

u 

(2.13) 

y-r  <  (•!  d’))  <  v-r 

(2.14) 

where  F“(v,d)G  argmin<max|(Fd). -vj  | 

YeY  1  '  ‘  1 

(2.15) 

Note  that  CCGIM  is  a  variant  of  this  model.  Because  target  assignments  are 
rigid,  F’ (v,d,Fjj  (•!  d))  =  (v,d)  =  7  .  Because  backorders  and  transshipments  are 

allowed,  r.  and  w.  may  be  negative  in  some  cases. 

CCNIM  is  not  a  PIP  for  several  reasons.  Not  only  does  each  system  of  chance 

constraints  allow  permutations  of  the  demand  vector,  but  there  are  |s’|  +  l  such  systems. 

Even  if  the  permutation  matrix  were  fixed  for  a  given  missile  allocation,  defining  relevant 

PEPs  on  the  joint  distribution  of  d'  and  d“  would  be  impossible.  However,  we  could 
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still  define  PEPs  separately  for  and  for  d“  following  each  period-I  scenario  5,  a  fact 
we  exploit  in  the  algorithms  we  develop  later. 

B.  INTEGER  PROGRAMING  MODELS 

Stochastic  programs  are  often  converted  into  mixed  integer  programs  and  solved 
using  standard  optimization  software,  e.g.  Ruszczyhski  [2002].  An  IP  representation  of 
CCNIM,  denoted  the  Fully  Flexible  Assignment  Model  (FFAM),  is  presented  in  Avital 
[2004].  FFAM’s  computational  performance  is  inadequate,  and  solution  times  for 
problems  of  modest  size  often  exceed  one  hour.  For  the  reader’s  convenience,  a  full 
description  of  FFAM  is  provided  in  Appendix  A. 

The  following  provides  a  formal,  non-mathematical  description  of  the  Reduced 
Fully  Flexible  Assignment  Model  (RFFAM),  a  new  IP  representation  of  CCNIM.  This 
model’s  formulation  requires  significantly  fewer  constraints  and  variables  than  FFAM, 
and  initial  results  indicate  it  solves  dramatically  faster  than  FFAM.  We  differentiate 
between  “regular”  constraints,  which  are  the  minimum  necessary  to  model  CCNIM,  and 
“specialized”  constraints,  which  enforce  the  Monotonic  Assignment  Policy  (MAP)  and 
add  valid  inequalities  that  reduce  the  feasible  region  and  improve  solution  times.  The 
specialized  constraints  in  RFFAM  are  adapted  from  those  originally  developed  to 
improve  FFAM’s  performance.  Some  of  the  specialized  constraints  require  data 
obtained  from  auxiliary  models;  we  denote  such  data  by  “auxiliary  parameters.”  A 
detailed  explanation  of  the  specialized  constraints  and  the  auxiliary  models  follows  the 
mathematical  description. 

I.  RFFAM  Specification 

The  following  non-mathematical  description  of  RFFAM  itemizes  a  list  of 
problem  requirements  followed  by  the  constraint  keys  where  they  are  represented  in  the 
mathematical  model  presented  in  the  next  section. 

Regular  constraints 

•  Minimize  the  weighted  cost  of  allocating  missiles  to  the  combat  ships  in  period  I 
and  storing  extra  missiles  at  a  depot  for  possible  use  in  period  II,  (2.16). 
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Subject  to: 

•  Each  ship’s  binary  allocation  variable  agrees  with  its  inventory,  (2.17)  and  (2.31). 

•  Each  ship  in  each  scenario  is  successful  only  if  it  has  enough  missiles  to  satisfy 
the  demand  of  its  assigned  mission,  (2.18)  and  (2.32). 

•  Each  scenario  in  each  period  is  satisfied  only  if  every  ship  has  enough  missiles  to 
satisfy  the  demand  of  its  assigned  mission,  (2.19)  and  (2.33). 

•  In  each  period,  the  probability  of  successfully  covering  the  scenarios  exceeds  a 
user-specified  threshold,  (2.20).  In  period  II,  the  cumulative  probability  must  be 
achieved  for  every  possible,  preceding,  period-I  scenario,  (2.34). 

•  Each  ship  is  allocated  a  specific  number  of  missiles  in  each  scenario,  (2.21)  and 
(2.35). 

•  Each  ship  in  each  scenario  is  assigned  exactly  one  mission,  (2.22)  and  (2.36). 

•  Each  mission  in  each  scenario  is  assigned  to  exactly  one  ship,  (2.23)  and  (2.37). 

•  Each  ship’s  inventory  following  the  prosecution  of  a  successful  mission  equals  its 
initial  level  less  the  demand  associated  with  its  assigned  mission,  (2.24)  and 
(2.25).  Otherwise,  the  remaining  inventory  is  zero,  (2.26).  (All  ships  “stay  and 
fight.”) 

•  Eor  each  ship,  its  remaining  inventory  following  the  prosecution  of  its  period-I 
mission  is  a  non-negative  integer,  bounded  by  the  physical  capacity  limit  (2.27). 

•  Eor  each  ship,  the  number  of  missiles  it  may  carry  at  the  outset  of  each  combat 
period  is  bounded  from  below  (discretionary  operational  constraint)  and  from 
above  (physical  capacity  limit),  (2.28),  (2.29),  (2.39)  and  (2.40). 

•  Eor  each  period-I  scenario,  each  ship’s  total  inventory  of  missiles,  before 
prosecuting  any  period- II  mission,  equals  the  post-mission  inventory  after  period  I 
plus  any  missiles  that  are  replenished,  (2.30). 

•  Eollowing  each  period-I  scenario,  the  total  number  of  missiles  distributed  to  the 
ships  between  periods  of  combat  may  not  exceed  the  number  kept  in  the  depot, 
(2.38). 

Specialized  constraints 

•  In  each  period-I  scenario,  missions  with  greater  demands  are  assigned  to  ships 
with  greater  inventories  (MAP  constraints),  (2.41)  and  (2.42). 

•  The  number  of  missiles  allocated  to  ship  i+\  does  not  exceed  the  number 
allocated  to  ship  i,  (symmetry-breaking  constraints)  (2.43). 

•  The  total  number  of  missiles  allocated  in  each  period  must  exceed  some  minimum 
(total-allocation  valid  inequalities),  (2.44),  and  (2.45). 

•  Each  ship  is  allocated  at  least  some  minimum  number  of  missiles  in  period  I 
(single-ship  valid  inequalities),  (2.46). 
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2.  RFFAM  Mathematical  Description 

ships 

level  (number)  of  missiles  ( ^  =  {O, . . . ,  max  vj  }  ) 
missions 

scenario  5  in  period  I 
scenario  s'  in  period  II 

Subset  of  period-II  scenarios  that  may  occur  (with 
positive  probability)  following  period-I  scenario  5 

Initial  Parameters  [units] 

d  demand  associated  with  mission  m  in  scenario  s  [missiles] 

(p^  probability  that  period-I  scenario  5  occurs 

conditional  probability  that  period-II  scenario  s'  occurs,  given  that 
period-I  scenario  5  occurs 

probability  threshold  for  period  I  (probability  that  the  realized 
scenario  must  be  satisfied) 

probability  threshold  for  period  II,  if  scenario  s  occurs  in  period  I 

Cj  unit  cost  of  procuring  a  missile  and  allocating  it  to  a  ship 

[$/missile] 

C2  unit  cost  of  procuring  a  missile  and  storing  it  in  the  depot 

[$/missile] 

V;  discretionary  lower  bound  on  the  number  of  missiles  that  may  be 

allocated  to  ship  i  [missiles] 

V.  physical  upper  bound  on  the  number  of  missiles  that  may  be 

allocated  to  ship  i  [missiles] 


Indices 

ie  7 
K 

me  M 
seS^ 
s'e  s“ 

S“(5) 


Auxiliary  Parameters  [units] 

minimum  aggregate  ship  load-out  required  to  satisfy  the  period-I 
scenarios  [missiles] 

minimum  aggregate  ship  load-out  required  to  satisfy  the  period-II 
scenarios  following  scenario  s  in  period  I  [missiles] 
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period-I  data- specific  lower  bound  on  the  number  of  missiles  that 
may  be  allocated  to  ship  i  [missiles] 


V: 


X:. 


iks 


Decision  Variables  [units] 

(Note:  As  discussed  later,  variables  vl ,  and  w.^  can  be  substituted  out.) 

v’  number  of  missiles  allocated  to  ship  i  in  period  I  [missiles] 

1  if  ship  i  has  k  missiles  before  the  first  combat  period,  and  0 
otherwise 

number  of  missiles  allocated  to  ship  i  in  period  II  following  period-I 
scenario  5  [missiles] 

1  if  ship  i  is  replenished  to  level  k  missiles  following  period-I 
scenario  5,  and  0  otherwise 

number  of  missiles  allocated  initially  to  the  central  depot  for 
potential  use  in  period  II  [missiles] 

r.^  number  of  missiles  left  on  ship  i  following  period-I  scenario  s 
[missiles] 

number  of  missiles  i  receives  from  the  depot  following  period-I 
scenario  5  [missiles] 

1  if  period-I  scenario  5  is  satisfied,  and  0  otherwise 

I  if  ship  i  successfully  covers  its  assigned  mission  in  period-I 
scenario  5  ,  and  0  otherwise 

zf,^  I  if  period-II  scenario  s'  is  successful  following  period-I  scenario  5, 
and  0  otherwise 

I  if  ship  i  successfully  covers  its  assigned  mission  in  period-II 
scenario  s'  following  period-I  scenario  5,  and  0  otherwise 

I  if  ship  i  is  assigned  mission  m  in  period-I  scenario  5,  and  0 
otherwise 

1  if  ship  i  is  assigned  mission  m  in  period-II  scenario  s'  following 
period-I  scenario  5,  and  0  otherwise 


Formulation 

Due  to  the  specialized  constraints  (2.41)  and  (2.42),  constraints  (2.22)  and  (2.23) 
are  redundant,  and  remain  in  RFFAM  for  exposition  purposes  only. 
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(RFFAM)  min  c.Yv^+cy  (2.16) 

r,u,x,v,w,z 

i 

S.  t. 

Period  I: 

v;=Z*^4  Vi  (2.17) 

k 

1  (  ^ 

zl<-  Yj  +1  m,  5e  s'  (2.18) 

) 

z]<zl  Vi,seS^  (2.19) 

Y^sZl>P  (2.20) 

Z4=l  Vi  (2.21) 

k 

X«L  =  1  Vi.sES'  (2.22) 

m 

XmL=1  Vm,5GS'  (2.23) 

i 

nM-'L‘^...>‘L+^p-i)  Vi.^ss'  (2.24) 

m 

ris  ^  v;  -  X  ^  ^  s'  (2.25) 

m 

>^,<v,zl  V/,.GS'  (2.26) 

/^,e{0,l,...,i^},  V/,  ^es'  (2.27) 

Vi^{vi’Vi  +  '^’---Y}  (2.28) 

4=0  V/,  fcl(/t<  v,.)v(fc>ii)  (2.29) 

zl,  zl,  ul^  binary. 

Period  II: 

Vi,s£S‘  (2.30) 

Z*4=v“  Vi,  SES'  (2.31) 

k 


Z^'s  ^  ^  Z  ^^“4  -  4/wL-,  +1  V/,  m,  5  e  s' ,  5  e  S" (5)  (2.32) 

\k>d,„'  j 
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yi,seS\s'eS\s) 

(2.33) 

5'gS^^  (5') 

(2.34) 

Z4.=l  V/,  5eS' 

k 

(2.35) 

V/,  5es\  5'es“(5) 

m 

(2.36) 

Z^.w.  =1  Vm,  5eS‘, /eS“(5) 

i 

(2.37) 

^w.  V^eS' 

i 

(2.38) 

V/,  5eS' 

(2.39) 

=  0  Vi,  A:|(A:<V;)v(A:>\^.),  5eS* 

(2.40) 

4>  wL',  binary. 

II  n? 

X  ^  £  . 

Specialized  constraints: 

u]i^=\  Vi,  5es‘ 

(2.41) 

“L=0  yi,m^i,sGS^ 

(2.42) 

Z4'^Z4i,r  V/<n-l, 

(2.43) 

Al 

"p^~ 

(2.44) 

Z4^4  V^eS' 

i 

(2.45) 

x\  =  0  V/,  k  <  y, 

(2.46) 

with  these  restrictions  on  the  data: 

1 

s 

VI 

> 

+ 

ls^~ 

Al 

ls^~ 

(2.47) 

Vm<n-1,  V5e  s‘us". 

(2.48) 
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This  model  contains  some  variables  that  can  be  substituted  out,  namely  , 

and  w-^.  Defining  these  variables  in  constraints  (2.17),  (2.31),  and  (2.38)  generates 
“branching  constraints,”  however,  and  branching  on  these  variables  accelerates  the 
branch- and-bound  solution  process  for  the  integer  model;  see  Appleget  and  Wood  [2000]. 

3.  Explanation  of  the  Specialized  Constraints 

It  often  happens  that  RFFAM  finds  an  optimal  integer  solution  fairly  quickly,  but 
the  solver  spends  a  great  deal  of  time  pruning  the  branch-and-bound  tree  before  it  can 
declare  that  solution  optimal.  The  specialized  constraints  included  in  RFFAM  help 
reduce  solution  times  by  substantially  tightening  the  linear-programming  (LP)  lower 
bound,  and  reducing  the  number  of  nodes  in  the  branch-and-bound  tree  that  must  be 
explored. 

Three  sets  of  constraints  restrict  target  assignments  to  reflect  the  MAP  (described 
shortly).  While  this  restricts  the  problem  (with  respect  to  CCNIM),  the  optimal  solution 
value,  in  our  experience,  seldom  increases  at  all  compared  with  that  of  the  unrestricted 
model.  Other  constraints  tighten  the  LP  relaxation  by  enforcing  lower  bounds  on  the 
number  of  missiles  that  each  ship  requires  individually  and  that  the  combat  force  requires 
as  a  whole  to  meet  the  demands  of  each  combat  period.  These  lower  bounds  are  found  by 
solving  instances  of  single-period  inventory  models. 

a.  Operational-Assignment  and  Symmetry -Breaking  Constraints 
Identifying  optimal  assignment  plans  is  difficult  for  RFFAM,  and 
obviously  would  also  be  a  very  difficult  task  for  a  force  commander  to  determine  in  real¬ 
time.  (We  assume  a  human  commander  assigns  the  targets,  although  an  automated 
decision  aid  could  also  be  used.)  In  all  likelihood,  the  commander  attempts  to  satisfy  the 
current  demands  by  following  some  intuitive  assignment  rules,  and  cannot  assess  the 
implications  of  each  possible  outcome  on  the  post-battle  distribution  of  missile 
inventories.  If  an  optimal  plan  requires  fewer  missiles  than  one  prescribed  by  any  simple 
rules,  then  it  is  unlikely  that  the  commander  will  be  able  to  identify  it  in  real-time.  By 
using  some  alternative  plan,  he  may  use  more  missiles  than  the  scenario  prescribed,  and 
the  number  of  missiles  in  the  depot  may  no  longer  be  sufficient  for  the  period-II 
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requirements.  By  incorporating  an  assignment  policy  in  RFFAM,  we  guarantee  that  the 
procurement  levels  can  cover  scenarios  with  the  required  probability  if  that  policy  is 
indeed  followed. 

A  plausible  assignment  rule,  which  we  call  the  Monotonic  Assignment 
Policy  (MAP),  assigns  targets  with  larger  demands  to  ships  with  larger  inventories.  As 
well  as  being  computationally  convenient,  the  MAP  is  desirable  from  an  operational 
standpoint.  Assignment  plans  that  follow  the  MAP  minimize  the  maximum  shortage 
over  all  targets.  As  a  result,  the  MAP  guarantees  that  all  the  missions  are  satisfied,  if  that 
is  at  all  possible.  Due  to  these  properties,  which  we  formally  prove  in  Chapter  III,  we 
believe  that  a  commander  would  have  little  motivation  to  choose  a  different  assignment 
plan. 

In  principle,  there  can  be  many  assignment  plans  (combinations  of 
variables)  that  conform  to  the  requirements  of  the  MAP.  If  we  assume  that  ships  are 
listed  in  order  of  non-increasing  upper  bounds  i  =  I,..., n,  we  can  express  the  MAP  in  an 
efficient  manner.  We  first  reorder  the  demands  in  each  scenario  in  non-increasing  order, 
so  that 

Vm<n-1,  5es‘us“.  (2.48) 

We  then  assign  the  missions  to  the  ships  in  that  order  for  period  I,  obtaining  the 
constraints 

ul=l  Vi,  5  e  s'  (2.41) 

“L  =  0  V/,  m^i,  seS\  (2.42) 

Period-I  allocations  are  now  guaranteed  to  be  non-increasing  in  the  ships’ 
index.  Experience  indicates  that  including  explicit  symmetry-breaking  constraints, 

Vi  <«-!,*,  (2.43) 

k  '>k  k  '>k 

which  force  an  allocation  plan  in  which  ships  with  a  lower  index  are  allocated  more 
missiles  than  ships  with  a  higher  index,  improves  the  model’s  solution  time. 

Similar  constraints  cannot  be  easily  applied  to  period  II  because  different 
permutations  of  the  period-II  allocation  plan  have  different  costs,  which  depend  on  the 
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number  of  missiles  that  are  added  to  the  remaining  inventory  from  period  1.  Hence,  only 
after  observing  the  period-I  remainders,  can  we  know  specifically  what  inventory  level  to 
allocate  to  each  ship. 


We  require  complex  constraints  to  ensure  that  the  inventory  of  ship  i  is  no 
less  than  that  of  ship  i'  if  ship  i'  is  assigned  mission  m  and  ship  i  is  assigned  any 
mission  m  <m  .  Otherwise,  there  is  no  practical  restriction  on  the  difference  between  the 
two  inventories.  Such  constraints  can  be  formulated  as  follows,  where  v  =  maxjy  }  : 


^  +  +  Vi,  /V/,  m>2,  s',  /e  s“(5).  (2.49) 


These  constraints  are  not  included  in  RFFAM  because  testing  indicates  that  RFFAM 
solves  faster  without  them.  In  any  case,  the  period-II  allocation  in  any  feasible  solution 
of  RFFAM  will  successfully  cover  the  period-II  demands  with  probability  pf  if  the 

MAP  is  followed  (see  Corollary  1  in  Chapter  III).  Therefore,  the  solution  of  RFFAM 
guarantees  a  sufficient  number  of  missiles  is  procured  if  the  MAP  is  followed  in  both 
periods. 

b.  Aggregate-Allocation  Valid  Inequalities 

The  aggregate-allocation  valid  inequalities  (integer  cutting  planes)  exploit 
the  solutions  of  single -period  problems.  Recall  that  is  the  minimum  aggregate  ship 
load-out  required  to  satisfy  the  period-I  scenarios  and  is  the  minimum  aggregate  ship 
load-out  required  to  satisfy  the  period-II  scenarios  following  scenario  s  in  period  I.  We 
can  use  b^  and  Z?"  to  generate  lower  bounds  on  the  number  of  missiles  that  must  be 

allocated  to  combat  ships  in  each  period.  Since  any  feasible  solution  to  the  two-period 
problem  must  satisfy  the  period-I  constraints,  we  obtain 

Y^v]>b\  (2.44) 


Similarly,  the  total  number  of  missiles  allocated  to  the  ships  in  period  II 
must  reach  at  least  following  period-I  scenario  5,  and  from  this  fact  we  obtain 


(2.45) 
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We  obtain  by  solving  a  single-period  model,  defined  by  constraints 
(2.16)-(2.21),  (2.29),  and  (2.41)-(2.43),  but  where  we  set  q  =1  and  Cj  =  0.  We  refer  to 

this  model  as  RFFAM-sp  (“sp”  stands  for  “single  period.”)  We  obtain  by  solving 
RFFAM-spII  (“spll”  stands  for  “single  period  on  period  II”),  a  similar  model  defined  on 
the  period- II  scenarios  s' ^  s"(5)  with  their  respective  demands,  probabilities  and 
specified  probability  thresholds.  RFFAM-spII  comprises  the  following  constraints  from 
RFFAM:  (2.31)-(2.37)  and  (2.39)-(2.40),  but  requires  the  modified  objective  function: 


i'GS  I 


We  also  include  a  version  of  the  MAP,  to  accelerate  solution  of  the  model: 

m“,=1  Vi,  5eS\ /eS“(5).  (2.51) 

The  explicit  formulation  of  both  models  appears  in  Appendix  B. 

Because  of  the  MAP,  RFFAM-sp  and  RFFAM-spII  are  both  PIPs 

(RFFAM-spII  must  first  be  decomposed  into  |s’|  problems,  each  with  a  single 

probabilistic  constraint)  and  could  be  solved  by  enumerating  PEPs,  by  a  specialized 
method  utilizing  PEPs,  or  with  the  MSP  algorithm  of  KPP.  However,  for  problems  of 
our  size,  they  can  be  solved  as  IPs  by  branch  and  bound.  The  computation  time  for  both 
models  is  negligible  compared  to  that  of  RFFAM  for  problems  requiring  more  than  a  few 
seconds  to  solve. 

c.  Single-Ship  Valid  Inequalities 

The  discretionary  lower  bound  on  a  single  ship’s  inventory,  set  in 
constraints  (2.28),  (2.29),  (2.39)  and  (2.40),  derives  from  generic  operational 
considerations.  Tighter  constraints  can  be  derived  from  specific  problem  data. 

We  define  a  modified  single-period  model  denoted  RFFAM-lb  (“lb” 
stands  for  “lower  bound”)  in  which  we  require  the  set  of  n  |s*  |  missions  to  be  assigned  to 

the  ships  as  before.  This  model  uses  the  same  constraints  as  RFFAM-sp,  except  that  we 
drop  constraints  (2.19)  and  replace  by  in  constraint  (2.20),  so  that  each  ship  selects 

its  set  of  successful  scenarios  independently  of  the  other  ships  (although  mission 
assignments  are  not  independent.)  This  relaxes  RFFAM-sp  because  each  ship  can  reach 
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the  specified  probability  threshold  by  satisfying  demands  from  a  set  of  scenarios 
different  from  those  of  another  ship. 


We  solve  RFFAM-lb  and  obtain  an  optimal  allocation  plan 
y  =  {yi,---,y^}  that  defines  the  number  of  missiles  k  allocated  to  each  ship  (as  specified 
for  RFFAM).  We  then  modify  the  lower  bounds  on  the  capacities  by  setting 

4  =  0  V/,  k  <  y,.  (2.46) 

before  solving  RFFAM. 

These  inequalities  are  valid  because  constraints  (2.41)  and  (2.42)  together 
completely  control  the  period-I  target  assignment  plans.  Therefore,  the  target-to-shooter 
assignment  plans  generated  by  RFFAM-lb  are  identical  to  those  generated  by  RFFAM 
in  period  I.  The  only  open  decision  actually  left  to  RFFAM  in  period  I  and  to  RFFAM- 
lb  is  the  choice  of  an  allocation  plan  with  a  sufficient  number  of  missiles  to  meet  the 
probability  threshold. 

Because  the  target  assignments  are  forced  through  the  assignment  rules, 
RFFAM-lb  is  a  PIP.  In  fact,  the  solution  of  RFFAM-lb  provides  the  same  lower  bound 
on  the  level  of  PEPs  for  RFFAM-sp  as  would  have  been  obtained  by  using  the  method 
described  in  Dentcheva  et  al.  [2000].  There,  the  lower  bounding  vector  on  any  PEP  is 
obtained  by  using  the  marginal  distribution  for  each  component  of  the  demand  vector. 

Equivalent  constraints  cannot  be  formulated  for  period  II,  because  the 
remainders  are  not  known  in  advance.  Eor  example,  consider  the  two  period-II  demand 
vectors  (6,3,1)  and  (5,3,3),  and  assume  that  only  one  of  them  must  be  satisfied.  Then 
solving  a  single  period  model  with  this  data  produces  the  optimal  allocation  (6,3,1). 
However,  if  the  three  ships  have  3  missiles  remaining  each,  following  the  period-I 
scenario,  then  it  would  be  less  costly  to  choose  to  satisfy  the  scenario  (5,3,3),  and 
replenish  2  missiles  from  the  depot,  rather  than  satisfying  the  scenario  (6,3,1),  which 
requires  3  additional  missiles.  If  (6,3,1)  were  used  as  single  ship  lower  bounds  in  period 
II,  the  minimum  allocation  would  be  (6,3,3),  which  is  excessive. 
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C.  SUMMARY  OF  COMPUTATIONAL  BEHAVIOR 

We  compare  the  solution  times  of  RFFAM  to  those  of  FFAM  using  the  CPLEX 
solver  version  9.0  [ILOG  2003]  on  a  Pentium  IV,  2  GHz  personal  computer  with  1  Gbyte 
of  RAM.  We  set  all  cost  ratios  using  the  smallest  integers  possible,  and  set  the  relative 
optimality  criterion  to  0.0%,  but  the  absolute  optimality  criterion  to  0.99. 

We  randomly  generate  scenarios  by  drawing  demands  from  a  discrete  uniform 
distribution  between  0  and  8  (details  on  test  cases  appear  in  Appendix  C).  All  of  the 
scenarios  are  equi-probable,  and  the  probability  threshold  in  each  period  requires  us  to 
satisfy  all  but  one  scenario  in  each  case.  We  refer  to  a  complete  specification  of  the 
parameters  required  to  solve  FFAM  as  an  instance.  We  generate  several  instances  from 
each  case  specification  by  varying  the  cost  ratio.  Not  only  does  this  change  the  objective 
function  and  result  in  different  optimal  solutions,  but  the  computational  effort  required  to 
solve  the  different  instances  changes  as  well.  We  obtain  the  parameters  required  for  the 
specialized  constraints  by  solving  one  instance  each  of  RFFAM-sp,  RFFAM-spII,  and 
RFFAM-lb.  The  auxiliary  models  solve  very  quickly,  so  we  ignore  this  computational 
effort  when  presenting  results. 

Table  1  summarizes  the  computation  time  of  both  FFAM  and  RFFAM  for  the 
four  cases  examined.  For  each  case,  we  list  the  number  of  ships  involved,  and  number  of 
period-I  and  period-II  scenarios.  For  each  instance,  we  report  the  solution  time,  in 
seconds,  for  FFAM  and  RFFAM. 

The  table  shows  that  FFAM’s  run-time  performance  is  too  slow  for  many 
practical  instances,  which  may  have  as  many  as  fifteen  ships,  and  often  more  than  ten 
scenarios  in  each  period.  RFFAM  outperforms  FFAM  in  every  one  of  the  instances 
presented,  usually  by  several  orders  of  magnitude.  In  fact,  RFFAM  solves  every 
instance  but  one  in  less  than  3  seconds.  In  that  instance  (case  2c,  with  cost  ratio  0.5)  a 
solution  within  3%  optimal  is  found  in  less  than  one  second. 
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Case  Name 

(2,1) 

(1,1) 

(1,2) 

FFAM 

RFFAM 

FFAM 

RFFAM 

FFAM 

RFFAM 

Case  2a 

(3,3,3) 

4.437 

0.359 

21.42 

0.249 

11.47 

0.249 

Case  2b 

(4,4,4) 

393.6 

0.656 

111.6 

1.453 

146.3 

0.328 

Case  2c 

(5,5,5) 

>3600" 

98.59 

>3600*’ 

2.812 

730.8 

2.046 

Case  2d 

(6,6,6) 

2511.2 

0.874 

>3600^^ 

1.374 

2839.1 

1.473 

Table  1.  Comparison  of  FFAM  and  RFFAM  Computation  Times. 

The  table  lists  the  amount  of  time,  in  seconds,  required  by  CPLEX, 
version  9.0,  to  solve  various  instances  of  FFAM  and  RFFAM.  The  top 
row  lists  the  cost  coefficients  used.  Each  subsequent  row  lists  the  times 
required  for  that  instance.  Case  size  is  specified  by  three  parameters, 

[n,  |s'|,  |s”|) .  Each  scenario  is  equi-probable,  and  one  scenario  in  each 

period  may  be  unsatisfied.  Instances  that  do  not  solve  in  an  hour  are 
indicated  by  >3600,  and  the  remaining  relative  gaps  are  as  follows: 
a)  4.23%.  b)  2.22%.  c)  6.28%. 

Although  RFFAM  is  much  faster  than  FFAM,  its  run  times  are  less  predictable 
than  FFAM’s  and  may  become  excessive.  As  is  apparent  from  the  table,  RFFAM 
requires  more  time  to  solve  the  instances  of  case  2c  than  those  of  case  2d  for  equivalent 
cost  ratios,  although  case  2c  has  fewer  ships  and  scenarios.  In  fact,  instances  of 
comparable  size  can  be  found,  for  which  the  solver  is  unable  to  report  an  optimal  solution 
following  one  hour  of  computation  (case  2e  with  costs  (2,1)).  In  a  larger  instance, 
involving  8  ships  and  8  scenario  in  each  period  (case  2f  with  costs  (1,  1)),  the  solver  fails 
to  even  identify  any  feasible  integer  solution  in  the  allotted  time.  Ultimately,  we  believe 
RFFAM  is  not  a  reliable  option  for  solving  practical  instances  of  CCNIM. 
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D.  CONCLUSIONS 

This  chapter  has  developed  CCNIM,  the  two-period  chance-constrained 
inventory  model  for  naval  surface  warfare.  The  model  is  a  relatively  complex  three-stage 
stochastic  program.  We  also  formulate  the  deterministic-equivalent  RFFAM,  an 
improved  formulation  over  FFAM,  which  was  developed  in  Avital  [2004].  Solution 
times  for  RFFAM  are  substantially  better  than  for  FFAM,  but  neither  model  can  reliably 
solve  problems  of  practical  size.  In  the  following  chapters,  we  develop  alternative,  faster 
solution  techniques. 
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III.  THE  MONOTONIC  ASSIGNMENT  POLICY 


We  introduced  the  MAP  in  Chapter  II,  which  is  the  mission  allocation  policy  that 
assigns  targets  with  larger  demands  to  ships  with  larger  inventories,  in  part,  to  improve 
the  solution  times  for  the  IP  models.  More  important  than  improving  solution  times,  we 
believe  the  MAP  is  a  reasonable  approximation  of  how  actual  assignments  should  and 
would  be  carried  out  because  assigning  targets  according  to  the  MAP  guarantees  that  all 
the  missions  are  satisfied,  if  that  is  at  all  possible.  Furthermore,  the  MAP  guarantees  the 
lowest  cost  solution  to  any  single-period  problem,  so  it  is  also  logistically  efficient.  For 
these  reasons,  we  believe  that  the  MAP  should  be  implemented  in  any  algorithm  that 
calculates  a  procurement  and  allocation  plan  for  naval  surface  warfare,  and  formulate 
CCNIM-MAP,  a  restriction  of  CCNIM  that  incorporates  MAP  constraints. 

A.  OPERATIONAL  PROPERTIES  OF  THE  MAP 

We  identify  three  properties  of  a  MAP,  which  hold  for  any  given  allocation. 

Property  3.1:  A  MAP  minimizes  the  maximum  shortage  over  all  targets  in  any 
scenario. 

Property  3.2:  A  MAP  satisfies  any  single-period  scenario  that  can  be  satisfied  by 
a  given  set  of  inventories. 

Property  3.3:  A  MAP  expends  the  greatest  number  of  missiles  in  any  scenario  (a 
desirable  characteristic  when  planning  procurement,  as  we  explain  later). 

We  proceed  to  prove  these  properties  in  several  lemmas,  below. 

Because  each  property  depends  only  on  the  relation  of  a  set  of  ships’  inventories 
to  the  ships’  assigned  missions,  we  assume,  without  loss  of  generality,  that  the 
inventories  are  non-increasing,  i.e.,  Vj  >  V2  >  •••  >  v,, ,  for  ships  i  =  Let 

d  =  (dp...,d„)G  Z"  represent  the  demands  of  an  arbitrary  scenario.  There  are  n\ 
permutations  of  the  elements  of  d .  Let  K  =  {l,...,n !}  be  the  index  set  of  all  permutations, 
and  let  d*" ,  G  K ,  denote  the  permutation  of  d .  We  denote  the  monotonically  ordered 
demand  vector  by  d* ,  and  let  k*  denote  some  corresponding  permutation  index.  Note 
that  d*  is  unique,  although,  due  to  ties  in  the  data,  there  may  be  several  corresponding 
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permutation  indices.  We  represent  alternative  assignment  plans  by  different  permutations 
of  the  demand  vector,  indexed  by  k,  and  assign  demand  to  ship  i . 

Lemma  1. 

Let  ve  Z"  be  a  non-increasing  vector.  Let  de  Z"  be  a  vector  of  demands,  and 
let  d*  be  a  permutation  of  d  whose  components  are  non-increasing.  Let 
=  maxjdf -vj  denote  the  greatest  single-target  shortage  associated  with  the 

assignment  plan  represented  by  d^ .  Then  S'"  =  min  S'‘ . 

teK 

Proof: 

Let  k  e  argmin  S'‘ ,  and  suppose  that  d^  d* .  Then,  there  exist  demands  I  and  m 

teK 

such  that  I  <m  and  df  <  .  Consider  d^  ,  an  alternative  permutation  of  d ,  in  which 

<'  =  <  (3.1) 

df  =  di  (3.2) 

df  =  dfyi^l,m.  (3.3) 

Let  (/,m)  =  maxjdf -vj  .  Then,  =  maxj^^  (/,m),  df  - v,,  d^ -v^ j  and 

d^-v,,  df  -v^j  .  By  our  assumptions,  we  have  d^  -v,  <  d^  -v„ 
and  df  -  v„  <  d^-v  ,  so  . 

If  d^  ^  d* ,  set  k  =  k' ,  and  find  a  new  permutation  d^  that  satisfies  conditions 
(3.1)-(3.3).  Because  |k|  is  finite  and  each  k'  indexes  a  different  permutation  of  d,  the 

sequence  of  permutation  indices  k,k',...,k''’^^  must  reach  d*  without  increasing  S'‘ .  | 

The  proof  of  the  following  corollary  is  then  obvious. 
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Corollary  1. 

Let  ve  Z"  be  a  non-increasing  vector.  Let  de  Z"  be  a  vector  of  demands,  and 

let  d*  be  a  permutation  of  d  whose  components  are  non-increasing.  If  v  >  d^  for  some 
permutation  index  k,  then  v  >  d* .  | 

Lemma  1  formalizes  property  3.1,  i.e.,  that  a  monotonic  plan  always  minimizes 
the  maximum  single-target  shortage.  In  particular,  that  shortage  is  zero  in  satisfied 
scenarios.  This  leads  directly  to  Corollary  1,  which  implies  that  if  any  successful 
assignment  plan  exists,  then  a  monotonic  plan  is  successful  as  well. 

As  an  example  of  Lemma  1,  consider  a  situation  with  three  targets,  corresponding 
demand  vector  d  =  (4,4,2),  and  an  inventory  vector  for  three  shooters  of  v  =  (4,3,1). 
Under  a  MAP,  the  target  is  assigned  to  the  ship,  z  =  1,  2,  3 .  The  maximum  single¬ 
target  deficiency  is  one  missile.  If  the  commander  were  to  switch  the  last  two 
assignments,  the  second  ship  would  not  use  its  entire  inventory,  but  the  second  target 
would  be  engaged  with  only  a  single  missile,  increasing  the  maximum  deficiency  to  three 
missiles. 

Finally,  Lemma  2  proves  property  3.3,  i.e.,  that  adhering  to  the  MAP  guarantees 
expending  at  least  as  many  missiles  in  any  scenario  as  any  other  assignment  plan,  without 
engaging  any  target  with  more  missiles  than  it  requires,  of  course.  As  a  result,  if  a 
different  assignment  plan  is  chosen  for  any  reason  (for  example,  when  the  assumption  of 
full-flexibility  does  not  hold),  the  total  number  of  missiles  left  over  from  period  I  does 
not  decrease.  However,  the  number  of  remaining  missiles  on  each  ship  may  be  different 
following  an  alternative  assignment.  Therefore,  some  ships  may  now  carry  more  missiles 
than  they  require  in  period  II,  while  other  ships  need  to  replenish  more  missiles  than 
originally  planned.  Because  of  the  no-transshipment  policy,  the  depot  inventory  may 
now  be  insufficient  for  period  II.  However,  by  recalling  some  “overstocked”  ships  to 
port  and  violating  the  no-transshipment  policy,  we  are  guaranteed  to  have  enough 
missiles  to  meet  the  period-II  required  allocations.  This  violation  may  occur  naturally, 
because  some  overstocked  ships  may  be  recalled  to  port  for  other  purposes.  Although 
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property  3.3  does  not  directly  affect  the  efficiency  of  any  algorithm,  it  implies  additional 
robustness  for  a  procurement  solution  provided  by  CCNIM,  and  further  motivates  the 
use  of  a  MAP. 

Lemma  2. 

Let  ve  Z"  be  a  non-increasing  vector.  Let  de  Z"  be  a  vector  of  demands,  and 
let  d*  be  a  permutation  of  d  whose  components  are  non-increasing.  Let 
//f  =  minjvpdf  I  be  the  number  of  missiles  fired  by  ship  i  associated  with  the 

assignment  plan  represented  by  d^ .  Then  ^  ?jf  =  max 

i  i 

Proof: 

Let  k  G  argmax  ^  Jjf  ,  and  suppose  that  d*  d* .  Then,  there  exist  demands  I  and 

ieK  i 

m  such  that  I  <m  and  df  <  .  Consider  d^  ,  an  alternative  permutation  of  d  that 

satisfies  conditions  (3.1)-(3.3).  We  obtain 

Z -  Z b’l  ’  b’>n ’ }]  -  [min  {v, ,  df  }  +  min  {v„ ,  d^ ]  .  (3.4) 

i  i 

Because  we  assume  v,  >  and  df  <  d^ ,  there  are  six  cases  to  examine: 

•  If  >  d* ,  then  equation  (3.4)  reduces  to  [d^  +  df  ]  -  [df  -I-  d^ ]  =  0 . 

•  If  V,  <  df  ,  then  we  are  left  with  [v,  -l-  ]  -  [v,  -l-  ]  =  0  . 

•  If  V,  >  d^  >  >  df  ,  we  have  fd^  -I-  df  1  -  Tdf  -I-  ]  =  d^  -  >  0 . 

•  If  v,>d^^>d^>v^,wehawe  [d;^ +v„]-[df +v„]  =  d;^ -df  >0. 

•  If  dt^Vi>d^  >v„,wehave  [v, +v^]-[df +v^]  =  v, -df  > 0. 

•  If  dt>vi>v^>d^  ,wehawe  [v, +df  ]-[df +v„]  =  v, -v„  >  0  . 
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If  ^  d* ,  set  k  =  k' ,  and  find  a  new  permutation  d^  that  satisfies  conditions 
(3.1)-(3.3).  Because  |k|  is  finite  and  each  k'  indexes  a  different  permutation  of  d,  the 

sequence  of  permutation  indices  k,k',...,k^'^^  must  reach  d*  without  decreasing  .  | 


B.  OPTIMAL  SOLUTION  TO  CCNIM-sp 

This  section  proves  that,  assuming  the  shooters  can  be  listed  such  that  v  and  v 
are  both  non-increasing,  then  assigning  targets  according  to  the  MAP  guarantees  an 
optimal  solution  to  CCNIM-sp.  We  believe  this  assumption  is  reasonable  in  practical 
situations,  as  will  be  discussed  below.  The  optimality  of  the  MAP  in  single-period 
problems  is  significant  because,  by  applying  the  MAP,  CCNIM-sp  is  reduced  to  a  PIP. 
The  specialized  algorithm  for  solving  CCNIM,  which  we  develop  later  in  this 
dissertation,  can  then  take  advantage  of  the  properties  of  PIPs  to  solve  problems  quickly. 

To  prove  that  assigning  targets  according  to  the  MAP  is  optimal  for  instances  of 
CCNIM-sp,  we  first  prove  that  the  MAP  is  optimal  for  the  Flexible  Minmax  Subset 
Problem  (FMSP).  FMSP  is  a  combinatorial  problem  that  is  based  on  the  Minmax  Subset 
Problem  (MSP)  defined  by  KPP  [2004].  We  obtain  FMSP  by  eliminating  the  capacity 
constraints  (2.3)  in  CCNIM-sp. 

Let  d^  =  Z"  denote  a  demand  vector  for  scenario  5GS  with 

associated  probability  (p^>0  .  Let  ={l,...,n!}  be  the  index  set  of  all  permutations  of 
dj ,  and  let  d^’ ,  g  ,  denote  the  vector  induced  by  the  k*  permutation  of  d^ .  Let 
k  =  ^A:j,...,k|g|  j  denote  a  vector  of  permutation  indices,  one  for  each  scenario,  which  we 
refer  to  as  an  arrangement  indicator.  Let  D*'  =  uii;-}  ,  3.n(i  let  A  —  k  .  Note  thst 

\A\  =  {nf. 


A  subset  of  scenario  indices  c  S  is  said  to  be  p-feasible  for  a  given  probability 


parameter  p  if  '^(p^>  p  .  Let  Sp 

seSf 


denote  the  set  of  p-feasible  subsets. 
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Let  )  =  max  ;  the  vector  d*^  exactly  prescribes  the  largest  demand  that  each 

ship  will  need  to  satisfy  in  the  scenarios  that  comprise  the  p-feasible  subset.  Also,  define 

n  _ 

D  ,  k )  =  ^  ) .  We  now  define  FMSP  precisely: 

/-I 

Definition:  FMSP  {Flexible  Minmax  Subset  Problem):  Find  a  p-feasible  subset  S  and 
associated  arrangement  indicator  k  such  that  D(s,k)  =  min^D(s^,k) ;  we  refer  to  S  as 

a  “minimal  p-feasible  subset.”  | 

We  may  also  formulate  FMSP  as: 

n 

D(s,k)=  min  V maxd,^’  (3.5) 

S.eSp.k'^  seS, 

■'  l=l  ’ 

s.t.  'Z(p,>p 

seSf 

Let  A]y[  =|k  I  d.;  >  V/  <  j,  s|  be  the  set  of  arrangement  indicators  where 
each  vector  d^'  is  non-increasing.  Because  monotonic  orderings  are  unique  to  within  ties 
in  the  data,  the  vectors  d^’  and  their  upper-bounding  vector  d*‘  are  unique,  regardless 
which  ke  A]y[  is  used.  We  denote  these  vectors  by  d*  and  d*  respectively,  by  D*  the 
arrangement  of  d* ,  and  by  k*  some  kGAjy;.  (If  the  data  exhibit  no  ties,  then 
A„=|k*j.)  Furthermore,  let  A^j^j  =  |k  I  V/ <  j|  be  the  set  of  arrangement 

indicators  for  which  the  elements  of  the  resulting  vector  of  maxima  d‘‘  form  a  non¬ 
increasing  sequence  (“MM”  stands  for  “monotonic  maxima”).  For  notational 
convenience,  we  drop  the  set  argument  if  d^  (s)  is  calculated  over  the  entire  set  S . 
Clearly,  A„  e  A^^^  e  A . 

We  proceed  to  prove  that  solving  FMSP  using  a  fixed  arrangement  D*, 
corresponding  to  a  monotonic  plan,  yields  the  same  optimal  value  as  an  unrestricted 

solution  of  FMSP,  i.e.,  minD(s.,k*)=  min  D(s.,k). 

S^eSp  '  ''  SpSSp.k  '  '' 
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Lemma  3. 


Let  D  =  UK}  be  a  set  of  demand  vectors,  and  let  D*  be  an  arrangement  of  those 
vectors  in  which  d*  is  non-increasing  for  all  s .  Recall  that  =  max  ,  and  let 

seS 

d*  =  maxd* .  Then  d  *  =  min  d^  V/ . 

seS  keA„M 

Proof: 

We  prove  the  lemma  by  showing  that  there  is  a  sequence  of  arrangement 
indicators  leading  from  any  k  e  to  some  g  ,  where  each  successive 
arrangement,  obtained  by  interchanging  elements  in  a  permutation  of  a  single  scenario, 
maintains  the  relation  d*^  >  •  •  •  >  * . 

Let  k  G  argmin  ,  and  suppose  that  k  ^  Aj^j .  Then,  there  exist  scenario  t  and 

demands  I  and  m  with  I  <m  and  d,*'  <  d*', .  Consider  D*'  ,  an  alternative  arrangement  of 
D*^ ,  in  which 

(3.6) 

4=4  (3.7) 

d-'is  -  I,  in  if  s  =  t  and  Vz  if  s^t .  (3.8) 

Let  S,  =S-{t},  and  note  that  d;*^  =max|d,.‘‘ (s,),  d^'j  Vz  .  By  our  assumptions, 

^  and  d,"'  <  d‘; .  Because  d^  >  d^; ,  we  obtain  d,‘‘  =  d,‘‘  (S, )  >  d*'  and  d,‘‘  >  d;^; . 
Therefore,  d,*^  =max|d,*‘ (s,),  d^;|=d,*‘ (3,)  =  ^^,*^ .  Because  d*(<d^;,  we  get 
d^  <  d^  ,  and,  considering  that  d,)^  =  d,*^  Vz  ,  we  obtain  d,.‘‘  <  d,.‘‘  Vz . 

Because  d^  ^d^ ,  d*^  may  no  longer  be  monotonic  and  it  may  be  true  that 

k'^A]y[„.  Let  7  =  argmaxjd^^  <  d,.*^  I .  If  j  0  ,  let  the  arrangement  indicator  k‘  be 

i>m 

such  that  d,^'  =  df}. ,  Vzzi  <  z  <  j  - 1,  V5  ,  d =  d^l  \/s ,  and  d,^”  =  df;  \/i<mvi>  j,  \/s  . 

IS  «Ti,j  o  >  ’  jS  fits  ^  IS  IS  'f  ’ 
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Thus,  we  are  effectively  permuting  the  entries  of  without  changing  their  values. 
Because  Vi  >  m  and  >d^,we  have  k‘  e  ,  and  set  k'  =  k‘ . 

If  k'^  A„,  set  k  =k',  and  find  a  new  arrangement  D*'  that  satisfies  conditions 
(3.6)-(3.8).  Because  |a„^|  is  finite  and  each  k'  leads  to  a  different  member  of  A^j^j 
(possibly  through  k' )  the  sequence  of  arrangement  indicators  k,k',...,k^'^'  must  reach 
some  k^'"*  g  A„  ,  and  d*^*  *  =  d*  by  definition.  | 

The  proof  of  the  following  corollary  is  then  obvious. 

Corollary  2. 

n  _ 

Let  D(s,k)  =  .  Then  D(s,k*)  =  min  D(s,k) .  | 


Lemma  4. 

For  a  fixed  set  of  vectors  D  and  corresponding  index  set  S,  the  arrangement 
D*  =  UK}  minimizes  the  sum  of  component-wise  maxima  over  all  arrangements,  i.e. 

D(s,k*)  =  minD(s,k) . 

Proof: 

Let  kG  argminD(s,k) ,  so  D(s,k)  <  D(s,k*) .  If  k^Aj^jj^j,  there  must  be 

keA 

indices  I  and  m  with  I  <m  and  .  Consider  D*'  ,  an  alternative  arrangement  of 

D ,  in  which 

<'=<;  V.  (3.9) 

<;  =  <'  V.  (3.10) 

V/9^/,mV5.  (3.11) 
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In  effect,  the  arrangement  D*^  switches  two  components  of  d*^  so  that  =d^, 
and  Obviously,  the  equality  D(s,k')  =  D(s,k)  holds. 

If  k'^  ^mm’  k  =  k' ,  and  find  a  new  arrangement  D*^  that  satisfies  conditions 
(3.9)-(3.11).  Because  |a|  is  finite  and  each  k'  leads  to  a  different  member  of  A  the 
sequence  of  arrangement  indicators  k,k',...,k^'^'  must  reach  some  k^*^' g  Aj^j^  such  that 
D^S,k^'"*  j  =  D(s,k).  By  Corollary  2,  D(s,k*)<  D^S,k^'"*  j,  and  so,  by  our  initial 
assumption,  D(s,k)  =  D(s,k*) .  | 

Theorem  1. 

Let  S  and  keA  be  an  optimal  solution  of  FMSP.  Let  k*  g  A„  and  let 
G  argminD(s^,k*) .  Then  D(s*  ,k*)  =  D(s,k) . 

S^SSp 

Proof: 

By  assumption,  D(s,k)  <  D(s^,k*) .  Consider  the  arrangement  D*  and  the 

n  _ 

resulting  D(s,k*)  =  (s) .  Lemma  4  implies  that  D(s,k)  >  D(s,k*) ,  and 

D(s,k*)  >  D(s^,k*)  follows  from  g  argminD(s^,k*) ,  so  D(s,k)  =  D(s^,k*) .  | 

Sy-GSp 

Note  that  the  results  of  Theorem  1  are  invariant  to  row  permutations  on  D* .  More 
accurately,  let  7r{»)  denote  a  specific  permutation.  If  for  every  dj'  gd‘‘,  dj’  =.;r(d*), 

then  D(s^,k)  =  D(s*  ,k*). 

We  now  obtain  the  conditions  required  to  extend  Theorem  1  to  CCNIM-sp.  The 
theorem  states  that  solving  FMSP  on  a  set  of  ordered  scenarios  yields  d*,  where 
Prjd*  >  d*|  >  p  ,  and  l^d*  is  minimal.  Obviously,  if  there  is  some  ordering  of  the  ships 

for  which  y  <  d*  <  v ,  then  we  could  set  v*  =  d* ,  and  v*  would  be  an  optimal  solution  to 
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CCNIM-sp.  If  that  condition  does  not  hold,  we  require  some  transformation  from  an 
optimal  solution  of  FMSP  to  that  of  CCNIM-sp. 

We  consider  the  effect  of  the  lower  bounds  first.  It  is  easy  to  see  that  if  y  ^  d*, 
then  setting  v  =  max|d*,y|  may  not  solve  CCNIM-sp  optimally.  For  example,  consider 

the  case  where  only  one  of  the  two  scenarios  (5,1)  and  (4,3)  must  be  successfully 
covered,  and  y  =  (3,3).  Then  d*=(5,l)  leads  to  the  allocation  (5,3),  but  the  optimal 
solution  to  CCNIM-sp  is  (4,3). 

We  transform  the  data  pertaining  to  FMSP  to  obtain  an  optimal  solution  that 
satisfies  the  lower  bounds.  If  D ,  (p,  and  p  define  an  instance  of  FMSP,  then  we  label  by 
FMSP-lb  an  instance  of  FMSP  on  the  following  data  (marked  by  the  prime  symbol): 
D'  =  DUy,  ^'=0.5,  q)[=0.5q)^  V5GS,and  p  =0.5{\+  p) .  By  Theorem  1,  we  arrange 

d'  monotonically  and  obtain  an  optimal  solution  d'’*  and  associated  minimal  p'  -feasible 
subset  s'f.  Because  '^cp'^Kp',  we  must  have  ves'j,  so  y<d'*.  Because 

^'=0.5^^  V5GS  and  =0.5p ,  we  have  ^(p^>p  so  s'^*-{y}  is  a  p -feasible 

seS'^ 

subset  in  FMSP.  Therefore,  we  may  set  v  =  d'*  as  an  optimal  solution  of  CCNIM-sp  as 
long  as  we  list  the  ships  monotonically  by  lower  bound,  such  that  Vj  >,...,>  y„,  and 

d'*  <  V .  But,  we  already  assume  that  every  scenario  is  feasible,  i.e.  there  is  some 
permutation  of  the  demands  Td^  that  satisfies  the  inequality  v  >  Td^.  By  Corollary  1,  if 

V  is  also  non-increasing,  then  v  >  d*  \/s . 

We  conclude  that  if  an  ordering  of  the  ships  exists  where  y  and  v  can  be  non¬ 
increasing  simultaneously,  we  can  optimally  solve  CCNIM-sp  as  a  simpler  single-stage 
problem  by  ordering  the  demand  vectors  monotonically  a  priori.  If  such  an  ordering 
does  not  exist,  then,  by  ordering  the  ships  such  that  v  is  non-increasing,  we  guarantee 
that  the  solution  of  FMSP-lb  is  a  feasible  solution  of  CCNIM-sp.  It  seems  reasonable  to 
assume  that  both  bounds  be  non-increasing  simultaneously,  because  the  lower  bounds  are 
generally  fixed  for  all  ships,  or  are  set  proportionally  to  the  capacities.  Note  that  the 
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resulting  allocation  plan  may  not  be  unique,  as  there  may  be  alternative  plans  of  equal 
cost  that  satisfy  a  different  subset  of  scenarios.  From  this  point  forward,  we  assume  that 
ships  are  listed  monotonically  by  capacities,  that  the  lower  bounds  vector  is  monotonic, 
and  that  all  demand  vectors  are  monotonically  ordered  unless  otherwise  specified. 

Of  course,  the  MAP  constitutes  a  restriction  of  CCNIM,  and  may  increase  the 
cost  of  an  optimal  solution  through  an  increased  depot  inventory.  When  calculating  the 
procurement  requirements  by  CCNIM,  an  unrestricted  assignment  plan  may  assign  ships 
with  low  inventories  to  high-demand  missions  in  period  I  in  order  to  spare  missiles  for 
period  II  in  scenarios  we  do  not  intend  to  satisfy.  Furthermore,  a  different  assignment 
plan  alters  the  numbers  of  missiles  remaining  on  individual  ships  in  satisfied  scenarios, 
and  may  reduce  the  number  of  missiles  left  on  ships  where  they  are  not  required  in  period 
II.  For  these  two  reasons,  assignments  that  lead  to  a  smaller  depot  requirement  may  be 
found  if  the  MAP  is  not  enforced. 

Returning  to  the  example  that  follows  Lemma  1,  consider  the  situation  with  three 
targets,  corresponding  demand  vector  d  =  (4,4,2),  and  an  inventory  vector  for  three 

shooters  of  v  =  (4,3,1).  Under  the  MAP,  no  ship  has  any  missiles  left  following  the 

engagement,  but  the  scenario  is  not  satisfied.  If  the  commander  were  to  assign  the  target 
requiring  two  missiles  to  the  ship  carrying  four  missiles,  then  there  would  be  two  missiles 
left  following  the  engagement.  These  missiles  could  reduce  the  depot  requirement  for 
period  II.  However,  as  mentioned  in  Chapter  II,  it  appears  that  in  many  cases  a 
monotonic  plan  is  optimal  for  CCNIM. 


C.  CCNIM-MAP 

Because  enforcing  the  MAP  is  desirable  from  an  operational  perspective,  and  the 
cost  of  its  enforcement  appears  to  be  small,  we  modify  CCNIM  to  incorporate  the  MAP 
in  “CCNIM-MAP”: 
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(CCNIM-MAP)  min_Ci;c'+C2x“  (3.12) 

s.t. 

Pr{v’ >d'}>;?'  (3.13) 

rv‘-.x:‘=0  (3.14) 

v<v'<v  (3.15) 

(3.16) 

where  r=(v‘-d’)  (3.17) 

v“  =r  +  w(r,F„(.|d‘))  (3.18) 

and  where  w(»)g  Z"  is  chosen  such  that 

Pr{v“  >  Cp  >  p:  g  s'  (3.19) 

l^w(r,F„(.|d'))-x“<0  (3.20) 

y-r<w(r,F„(.|d'))<v-r  (3.21) 

where  (v)g  argjfy^'v)  >0  V/<n-l|  (3.22) 

and  with  these  restrictions  on  the  data: 

<  >  ,  Vi  <  n  - 1,  V5  G  s'  (3.23) 

df.  >  V/  <  n  -1,  V/g  s"  .  (3.24) 

y>y.^jV/<n-l  (3.25) 


CCNIM-MAP  is  only  a  two-stage  model,  and  is,  therefore,  simpler  than 
CCNIM.  The  flexibility  in  target  assignment  allows  us  to  preprocess  the  demand 
vectors,  and  order  each  one  monotonic  ally.  We  assign  demand  d]^  to  ship  i  in  every 
period-I  scenario,  and  calculate  the  remaining  inventory  r.  on  each  ship.  As  in  CCNIM, 
we  then  decide  the  number  of  missiles  wf  to  add  to  each  ship’s  inventory  for  period  II. 
This  determines  v" ,  the  period-II  inventories.  To  enforce  the  MAP,  some  permutation 
T^5['^p(v),  which  assigns  larger  demands  to  ships  with  larger  inventories,  must  be 
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determined.  Because  demands  are  ordered  monotonically,  this  permutation  is  simply  the 
inverse  of  the  permutation,  which  arranges  the  ships  inventories  in  non-increasing  order. 
An  inverse  of  a  permutation  matrix  is  well  defined,  because  its  determinant  is  either  1  or 
-1.  The  remainder  of  this  dissertation  is  dedicated  to  finding  practical  solution  methods 

to  CCNIM-MAP. 
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IV. 


IDENTIFYING  OPTIMAL  SOLUTIONS 


The  IP  model  RFFAM  may  fail  to  reach  an  optimal  solution  in  a  reasonable 
length  of  time  for  problem  instances  that  involve  more  than  a  few  ships  and/or  scenarios, 
depending  in  the  case  specifications  and  the  cost  ratio  used.  Therefore,  this  chapter 
explores  the  properties  of  feasible  solutions  of  CCNIM-MAP,  and  identifies  several 
specific  “points  of  interest,”  i.e.,  a  small  set  of  potentially  optimal  solutions  that 
correspond  to  optimal  solutions  for  certain  cost  ratios.  These  points  may  be  found  by 
solving  instances  of  RFFAM-sp  and  restrictions  of  RFFAM.  However,  these  methods 
prove  unreliable,  so  in  the  next  chapter  we  develop  a  specialized  algorithm  to  identify 
these  points  based  on  the  properties  of  PEPs. 

A.  STRUCTURE 

The  objective  function  (3.12)  is  the  weighted  sum  of  two  variables,  and  v”, 
where  is  the  total  number  of  missiles  assigned  to  the  ships  in  the  period  I,  and  is 
the  number  of  missiles  stored  at  the  depot.  We  define  an  integer  pair  x  =  to  be 

feasible  if  there  exist  missile  allocations  and  target  assignments  that  satisfy  the 
probability  requirements  in  CCNIM-MAP,  with  corresponding  values  of  x^  and  . 
Because  the  objective  function  is  entirely  determined  by  the  value  of  x ,  we  go  as  far  as 
to  refer  to  x  as  a  solution  to  CCNIM-MAP  in  the  following  text  (although  knowing  x 
does  not  reveal  all  of  the  other  variables  involved).  All  other  variables  are  treated  as 
secondary,  and  are  needed  only  to  assess  the  feasibility  of  x .  We  refer  to  the  value  of  the 
point  X  as  cx  =  c^x^  +  Cjv" .  In  much  of  the  following  discussion,  we  refer  to  xl  as  a 

first-period  optimal  solution  and  to  the  integer  pair  x*  =  (v’,  v")  as  an  optimal  solution  of 
the  full  problem. 

We  define  the  quasi-feasible  region  for  CCNIM-MAP,  F  cz  ,  as  the  set  of  all 
feasible  integer  pairs  .  We  bound  F  by  aggregating  the  individual  ship-capacity 

constraints.  Combining  these  with  the  lower  bounds  obtained  for  x^  and  v”,  and  the 
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non-negativity  restriction  on  ,  we  conclude  that  F  is  contained  within  the  rectangle 
defined  by  <\^ y  and  0 <  <  maxZ^f  .  Note  that  is  not  explicitly  bounded  in 

CCNIM-MAP,  but  will  not  exceed  maxZ?“,  the  greatest  number  of  missiles  we  may 
need  in  period  II,  in  any  actual  solution. 

We  also  define  the  following  subset  of  F , 

0(f)  =  {(x’,x“)|(x',x")g  F,  (x’-l,x“)^  F,  (x’,x“-l)^  f}.  (4.1) 

Only  the  points  in  o(f)  are  potentially  optimal,  because  both  points  and 

have  a  lower  objective  value  (assuming  positive  costs)  than  If 

either  of  those  solutions  is  feasible,  then  it  is  obviously  preferable  to  .  In  other 

words,  any  optimal  solution  of  CCNIM-MAP  must  be  Pareto  optimal  [Rardin  1998,  p. 
379]  with  respect  to  the  two  variables  v’  and  .  To  solve  CCNIM-MAP,  we  only  need 
to  identify  which  xg  o(f)  is  optimal,  and  in  many  cases  we  can  do  so  without  solving 
RFFAM. 

We  explore  the  properties  of  o(f)  to  find  features  that  will  lead  us  to  specialized 
algorithms.  KPP  use  a  specialized  decomposition  algorithm  to  solve  CCGIM,  so  we 
begin  our  exploration  by  examining  the  structure  of  o(f)  for  that  model,  which  is 

simpler  than  CCNIM-MAP.  In  CCGIM,  safety  stocks  exist  that  must  be  replenished,  if 
used,  and  units  may  transship  between  themselves  at  the  end  of  period  I.  Because  of  the 
safety  stocks,  for  any  given  period-I  scenario,  all  allocation  plans  expend  the  same 
amount  of  ordnance.  Furthermore,  KPP  prove  that  every  two-period  optimal  solution 
requires  the  same  amount  of  ordnance  in  total.  It  is  clear  that  o(f)  ,  in  the  KPP  setting, 

is  a  line  segment  in  ,  maintaining  the  relation  =  const . 

The  geometric  view  provides  further  intuition  into  the  correctness  of  the  KPP 
decomposition  procedure.  KPP  seek  the  point  =  argmin  v’ ,  which  is  the  potentially 

(x'.F)eo(F) 

optimal  point  that  minimizes  the  number  of  missiles  in  the  initial  (period  I)  ship 
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allocation.  They  obtain  it  by  decomposing  the  two-period  model  into  two  separate  one- 
period  problems,  which  they  solve  sequentially.  In  the  first  problem,  they  minimize  , 
the  amount  of  ordnance  allocated  to  the  units  in  order  to  satisfy  demands  in  the  first 
combat  period  with  a  specified  probability.  The  specific  period-I  allocations  are  used  as 
parameters  in  the  period-II  problems,  where  KPP  minimize  the  amount  of  ordnance  to  be 
added  to  the  units’  inventories  from  the  depot  following  each  possible  period-I 

scenario  5.  By  following  this  procedure,  they  obtain  v“=max|l^wl.  By  the 
assumption  that  Cj  >  C2 ,  must  be  globally  optimal  for  CCGIM. 

The  structure  of  o(f)  may  be  more  complex  in  the  case  of  CCNIM-MAP.  We 
define  several  specific  integer  pairs,  and  show  which  of  these  are  optimal  under  which 
circumstances.  As  above,  let  =  argmin  be  the  potentially  optimal  point  that 

(A:‘,x")eO(F) 

minimizes  the  number  of  missiles  in  the  initial  (period  I)  ship  allocation,  and  let 
x^  =  argmin  be  the  potentially  optimal  point  that  minimizes  the  number  of  missiles 

(x‘,x")eo(F) 

Stored  in  the  depot.  Let  x(f)=  argmin  be  the  set  of  points  that  minimize  the 

(x‘,x“  jeo(F) 

total  number  of  missiles. 

By  these  definitions,  it  is  obvious  that  in  cases  where  there  is  a  clear  preference 
for  allocation  strategy,  that  is,  when  the  ratio  Cj/Cj  is  very  different  from  1,  we  seek 

either  the  point  x^  or  the  point  x^  as  the  optimal  solution.  (It  is  difficult  to  determine 
how  large  or  small  “different  from  1”  must  be  without  solving  CCNIM-MAP.)  When 
Cj  /Cj  =  1 ,  any  point  xg  x(f)  is  an  optimal  solution.  When  the  cost  ratio  is  set  such  that 

Cj  /Cj  ~  1 ,  then  identifying  the  optimal  points  in  o(f)  may  be  more  difficult. 

The  first  point  of  note  is  that  x^  may  not  belong  to  x(f)  ,  in  which  case  it  is  not 
optimal  when  Cj  =  Cj .  As  an  example,  consider  case  4a  in  table  2,  which  has  five  combat 
ships,  s'  ={5p52’'^3}  s"  ={54}  .  In  this  case,  x^  =(7,2)  because  seven  missiles  are 
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enough  to  satisfy  scenario  5j  and  meet  the  probability  requirement.  If  that  scenario 

actually  occurs,  we  will  have  no  inventory  carried  over  into  period  II,  so  two  missiles  are 
needed  in  the  depot.  Increasing  the  period-I  allocation  total  to  eight  allows  the  period-I 
allocation  v'  =  (2,2,2, 1,1) ,  which  covers  scenarios  and  53 .  If  either  of  these  scenarios 

occurs,  we  obtain  the  (reordered)  remainder  vectors  (1, 1, 0,0,0)  or  (1, 1, 1, 0,0), 
respectively,  both  of  which  satisfy  the  period-II  requirement  without  requiring  any 
additional  missiles.  If  scenario  occurs,  the  remainder  is  (2,2, 1, 1,0).  Hence  we  have 

the  optimal  solution  (8,0) ,  which  requires  one  missile  less,  in  total,  than  . 

The  reason  that  may  not  belong  to  x(f)  is  that  a  ship  carries  no  safety  stock. 

Therefore,  the  number  of  missiles  expended  in  unsatisfied  scenarios  depends  on  the 
allocation.  Furthermore,  adding  missiles  to  the  ships’  inventories  in  period  I  relaxes  the 
limitations  on  the  scenarios  we  may  choose  to  satisfy  in  that  period.  We  may  now  choose 
a  different  allocation  plan,  requiring  more  than  x\  missiles,  that  uses  fewer  missiles  in 

the  unsuccessful  period-I  scenarios,  thereby  saving  more  missiles  for  period-II  and 
lowering  the  requirement  for  depot  stockpile.  The  difference  between  the  depot 
requirements  may  exceed  the  difference  between  the  total  period-I  allocations.  We  refer 
to  this  effect  as  the  no-safety-stock  effect  in  future  discussion. 

The  reader  may  correctly  assume  that  an  enumeration  procedure,  which  generates 
every  pkfeasible  subset,  can  overcome  these  difficulties.  Such  an  algorithm,  however,  is 
inefficient  compared  with  the  algorithm  we  provide  in  the  next  chapter,  and  solves 
CCNIM-MAP  only  for  cost  ratios  where  Cj/Cj  <1.  We  detail  such  an  algorithm  in 
Appendix  D. 
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Vi 
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Vi 
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Table  2.  Case  4a  -  Example  Where  ^  • 


The  table  lists  the  probabilities  of  each  scenario  {(p^,  the  success 
thresholds  specified  for  each  period  (p\  p"),  the  demands  associated 
with  the  various  missions  in  each  scenario  and  the  capacity  limits 
imposed  on  each  ship  ( y,,  y ).  A  minimum  period-I  allocation  leads  to  the 
solution  =  (7,2) .  If  we  allow  8  missiles  to  be  allocated  to  the  ships  in 
the  first  period,  we  can  choose  an  allocation  that  satisfies  the  scenarios 
and  53 ,  instead  of  scenario  .  The  remaining  inventory  following  any  of 
the  three  scenarios  will  satisfy  scenario  s\ . 


The  no-safety-stock  effect  causes  a  second  difficulty  in  the  analysis  of  o(f). 

CCNIM-MAP  assumes  that  ships  will  expend  their  entire  inventory  when  the  assigned 
demand  exceeds  that  inventory.  Increasing  the  allocation  to  any  such  ship  will  increase 
the  number  of  missiles  it  expends  in  such  a  scenario.  However,  if  the  original  allocation 
plan  was  already  p* -feasible,  then  increasing  it  constitutes  a  “waste”  of  missiles  (from 
the  logistics  point  of  view),  because  those  missiles  will  not  necessarily  reduce  the  depot 
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requirements  for  period  11.  Thus  we  may  have  cases  where  but 

(v*  +l,v"  -l)d  F . 

The  point  may  also  not  belong  to  the  set  x(f),  due  to  the  no-safety-stock 

effect,  as  well  as  the  no-trans shipment  effect.  The  no-trans shipment  effect  arises 
because,  following  period  I,  the  number  of  missiles  remaining  on  some  ships  may  exceed 
the  number  of  missiles  actually  required  by  them  to  satisfy  period-II  scenarios.  By 
assumption,  these  extra  missiles  are  essentially  wasted  because  they  cannot  be 
transshipped  to  other  ships  that  may  require  them. 

We  illustrate  these  effects  by  examining  o(f)  for  case  4b  in  Table  3.  There, 
x^=(15,15)g  x(f)  ,  and  the  initial  allocation  of  (5,4,4,2)  satisfies  every  scenario  but  s^, 
whose  associated  demand  vector  is  (5,4,4,3)  when  ordered  monotonically.  We  also  know 
that  (18,12)g  x(f)  ,  and  the  reader  may  verify  that  (19,11)  is  not.  That  is  because  there 

is  no  p’ -feasible  allocation  of  19  missiles  that  does  not  expend  at  least  16  missiles  in  any 
monotonic  assignment  plan  if  scenario  54  occurs.  In  fact,  there  are  exactly  five 
monotonic  19-missile  allocations  that  would  expend  only  16  missiles,  e.g.,  (7,4,4,4). 
There  are  three  total  missiles  remaining  following  54 ,  just  as  when  we  optimally  allocate 

18  missiles.  Hence,  the  depot  requirement,  which  is  driven  by  s^,  is  still  12,  and  the 
point  (19,12)  is  feasible.  Note  that  (19,12)^  o(^)  because  (18,12)g  o(f)  ,  so  the  next 
point  (in  order  of  decreasing  x" )  that  belongs  to  0  (f)  is  (20,1 1),  and  the  sum  of  missiles 
expended  has  increased  by  one.  The  reader  may  verify  that  a  similar  increase  in  total 
expenditure  occurs  from  the  point  (26,5)g  0(f)  to  (28,4)g  o(f),  because  27  missiles 

cannot  be  allocated  in  period  I  without  expending  17  missiles  when  scenario  s^  occurs. 
Thus,  we  require  the  same  depot  level  as  when  we  allocate  26  missiles  optimally,  and 
(27,4)^  0(f). 
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Table  3.  Parameter  Specifications  for  Case  4b. 

All  labels  defined  as  in  Table  2. 

The  increases  in  the  number  of  missiles  expended,  shown  above,  result  from  the 
no-safety-stock  effect.  Had  there  been  some  safety  stock,  which  must  be  replenished,  the 
ships  would  have  expended  17  missiles  whenever  scenario  occurred,  regardless  of  the 

initial  allocation.  It  is  perhaps  counter-intuitive  that  optimal  allocations  do  not  satisfy 
every  demand,  even  when  enough  missiles  are  allocated  in  total  for  that  period. 

When  minimizing  it  happens  that  (32,0)^  o(f)  and  x^=(32,l).  This 
further  increase  in  total  missile  requirement  results  from  the  no-transshipment  effect,  and 
occurs  when  increases  from  29  missiles  to  30,  requiring  a  depot  inventory  of  3 
missiles  in  both  cases.  The  period-II  scenario  specifications  are  independent  of  the 
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period-I  scenarios,  and  Z7“=15  for  all  s,  corresponding  to  the  period- II  minimal 
allocation  of  (5, 4, 4, 2).  There  are  only  two  (monotonic)  period-I  allocations  using  30 
missiles:  (8, 8, 8, 6)  and  (8, 8, 7, 7).  If  scenario  occurs,  the  remainders  are  (3,3,4,3)  and 

(3, 3,3,4),  respectively,  which  are  essentially  equivalent.  Three  more  missiles  are 
required,  in  total,  to  increase  the  inventories  of  three  ships  to  the  requirement  of  the 
period- II  minimal  allocation,  but  the  fourth  ship,  which  requires  an  inventory  level  of  two 
missiles,  is  also  carrying  three  and  has  one  in  excess.  29  missiles  can  be  initially 
allocated  according  to  (8, 8, 8,5),  leaving  (3, 3,4,2)  missiles  if  scenario  occurs,  and  still 
requiring  a  depot  inventory  of  three  missiles. 

Figure  2  gives  a  full  mapping  of  o(f)  for  case  4b  of  CCNIM-MAP,  based  on 

the  discussion  above.  Based  on  the  above  analysis,  we  find  all  of  the  optimal  solutions 
for  any  cost  ratio,  and  summarize  the  results  in  Table  4.  Some  solutions  cover  a  span  of 
ratios,  while  others  are  optimal  for  only  a  specific  cost  ratio;  these  latter  solutions  are 
never  uniquely  optimal.  The  solutions  obtained  by  solving  different  instances  of  case  4b 
in  RFFAM,  differentiated  by  cost-ratios,  correspond  to  some  of  the  points  of  interest  in 
o(f)  ,  as  can  be  seen  in  Table  4. 
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Figure  2.  A  Map  of  0  (f)  for  Case  4b  in  Table  3. 

The  figure  shows  every  point  belonging  to  o(f)  when  solving  case  4b. 
Only  these  points  could  lead  to  optimal  solutions  of  CCNIM-MAP. 
There  are  three  occasions  where  o(f)  breaks  from  linearity.  Two  breaks 

occur  due  to  the  no-safety-stock  effect,  which  results  in  no  feasible 
solution  to  CCNIM-MAP  if  we  use  the  total  period-I  allocation  and  depot 
inventory  pairs  (19,11)  or  (27,4).  The  pair  (30,2)  also  leads  to  no  feasible 
solution  of  CCNIM-MAP,  due  to  the  no-transshipment  effect. 
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C2/C 

Optimal  Solutions,  x* 

RFFAM  Initial  Allocation 

<1 

(15,15) 

(5,4,4,2) 

1.0 

(15,15),  (16,14),  (17,13),  (18,12) 

(8,4,4,2) 

A 

A 

00 

(18,12) 

(8,4,4,2) 

8/7 

(18,12),  (26,5) 

(8,4,4,2) 

8/7<---<3/2 

(26,5) 

(8,8,8,2) 

3/2 

(26,5),  (29,3),  (32,1) 

(8,8,8,8) 

>3/2 

(32,1) 

(8,8,8,8) 

Table  4.  Optimal  Solutions  for  Case  4b  in  Table  3. 


The  “Cj/cj  ”  column  gives  the  cost  ratio,  thereby  completing  the  definition 
of  the  instance.  The  “Optimal  Solution”  column  gives  the  optimal  number 
of  missiles  to  be  loaded  for  period-I  combat  and  the  number  of  missiles  to 
be  stored  for  later  use.  “RFFAM  Initial  Allocation”  gives  initial  load-outs 
for  the  ships,  as  obtained  by  solving  RFFAM  for  the  corresponding 
instance  of  case  4b.  Notice  there  are  several  specific  ratios  that  have 
multiple  optimal  solutions. 


Incidentally,  if  the  cost  ratio  is  greater  than  1.5,  the  optimal  solution  (32,1)  is  a 
direct  result  of  the  MAP,  and  is  of  strictly  higher  cost  than  the  optimal  solution  in  the 
unrestricted  case.  Consider  the  initial  allocation  (8, 8, 8,7),  requiring  a  total  of  31  missiles. 
If  the  targets  of  scenario  are  assigned  according  to  (5,4,3, 5),  then  the  remainders  are 

(3, 4,5, 2),  and  only  a  single  extra  missile  must  be  stored  at  the  depot.  For  any  other 
period-I  scenario,  no  more  missiles  would  be  required.  Hence,  the  solution  (31,1)  is 
feasible,  and  costs  less  than  =  (32,1) . 

Finally,  note  that  loading  the  ships  to  full  capacity  guarantees  that  the  resulting 
depot  inventory  is  minimal.  However,  it  may  be  that  v- A:, F  for  some  integer 

A:  >  0 ,  so  the  point  (l^v,v“)  may  not  belong  to  o(f)  .  This  can  be  caused  by  either  the 
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no-transshipment  effect  or  by  the  no-safety-stock  effect  because  both  effects  cause 
portions  of  the  convex  hull  of  F  to  parallel  the  v’  axis. 

B.  EXAMPLE  OF  REF  AM-BASED  TECHNIQUES 

As  observed  in  Chapter  II,  case  2f  is  difficult  to  solve  by  RFFAM.  In  this 
section,  we  attempt  to  solve  it  for  various  cost  ratios  by  exploiting  our  newly  gained 
insight  into  the  feasible  region  and  by  exploiting  other  techniques  involving  variations  of 

RFFAM. 

The  points  and  correspond  to  “preemptive”  allocation  policies  that 
respectively  minimize  the  total  number  of  missiles  allocated  to  the  ships  in  period  I,  or 
minimize  the  number  of  missiles  that  are  kept  in  the  depot.  We  can  try  to  use  RFFAM 
to  identify  x^  and  x^  directly  by  setting  extreme  cost  ratios.  However,  even  setting  a 

cost  ratio  of  100  to  1,  in  either  direction,  does  not  yield  an  integer  solution  in  one  hour  of 
computation. 

We  may  be  able  to  identify  x^  and  x^  significantly  faster  by  taking  advantage  of 
their  preemptive  nature.  The  point  x^  corresponds  to  a  solution  that  minimizes  the 
number  of  missiles  allocated  to  ships  in  period  I,  and  then  minimizes  the  number  of 
missiles  held  in  inventory  for  that  allocation.  We  know  that  x\=b^ ,  where  ,  the 

minimum  feasible  v’ ,  is  obtained  by  solving  RFFAM-sp.  Let  (v')  denote  the  value 

of  in  an  optimal  solution  of  RFFAM  where  is  fixed.  Because  is  unbounded,  a 
solution  [b^ )  must  exist.  We  expect  that  solving  RFFAM  with  the  added  constraint 

')Jy  =  b^  will  yield  the  point  x^  =  faster  than  by  setting  Cj/Cj^cl  in 

RFFAM.  However,  for  case  2f,  an  integer  solution  is  not  obtained  in  one  hour  even 
when  RFFAM  is  restricted  in  that  manner. 

If  we  prefer  to  reduce  the  operational  burden  of  replenishment  operations,  we 
seek  the  solution  x^  ,  which  first  minimizes  and  then  minimizes  x^ .  To  do  so,  we  can 

solve  RFFAM  twice.  First,  we  fix  an  assignment  plan  that  sets  x^  to  its  upper  bound  j’ 
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by  fixing  v,.  =  v.  V/ ,  so  that  each  ship  is  loaded  to  capacity.  Solving  RFFAM,  we  obtain 
=  v"  ,  the  minimum  depot  inventory.  We  now  fix  and  solve  RFFAM  to 

obtain  j ,  where  denotes  the  value  of  in  the  solution  of 

RFFAM  for  a  fixed  value  of  v” .  Obtaining  a  solution  by  this  method  should  be 
significantly  faster  than  simply  solving  RFFAM  with  Cj/cj  » 1  because  the  number  of 
potential  allocations  is  reduced.  When  we  fix  the  allocations  of  the  8  ships  in  case  2f  to 
their  upper  bounds,  we  obtain  the  solution  x“=19  in  5.3  seconds  of  computation. 
However,  when  we  solve  RFFAM  again,  with  the  depot  inventory  fixed,  no  integer 
solution  is  found  in  one  hour. 

Although  the  restrictions  designed  to  yield  and  x^  do  not  improve  the 

performance  of  RFFAM  significantly,  we  find  other  restrictions  that  do.  In  particular, 
one  restriction,  denoted  RFFAM-mll  (“mil”  stands  for  monotonic  in  period  II),  solves 
case  2f  in  less  that  one  second  for  every  cost  ratio.  In  RFFAM-mll,  we  restrict  the 
period- II  assignment  variables  by  adding  the  constraints 

=  0  Vz,  m  z,  5  e  s\  s'e  s“(5) .  (4.2) 

Recall  that  is  set  to  1  if  ship  z  is  assigned  mission  m  in  period-II  scenario  s', 

following  scenario  s  in  period  I,  and  0  otherwise.  Because  the  targets  have  been 
monotonically  ordered  by  demands,  ship  z  is  now  assigned  the  z^**  largest  demand  in  each 
period.  This  restriction  may  not  be  optimal  because  it  ignores  the  numbers  of  missiles 
remaining  on  each  ship  following  period  I. 

Lower  bounds  on  the  minimum  number  of  missiles  that  must  be  procured  to 
satisfy  RFFAM  can  be  found  by  relaxing  the  integrality  constraints  on  some  of  the 
variables  in  RFFAM.  We  suggest  removing  the  binary  requirement  from  the  variables 
controlling  period-II  assignments.  An  initial  attempt  to  solve  case  2f  with  zz“ g  [0,1] 

yields  no  optimal  solutions  in  one  hour  of  computation,  but  when  zz“ is  unrestricted,  an 

optimal  solution  is  obtained  in  about  four  seconds  for  every  cost  ratio.  We  refer  to  this 
model  as  RFFAM-rx  (“rx”  stands  for  relaxation).  RFFAM-rx  solves  significantly  faster 
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than  RFFAM  itself,  provides  an  integer  allocation  plan,  and  the  optimal  relaxed  objective 
value  provides  a  lower  bound  on  the  number  of  required  missiles. 

For  cost  ratios  that  satisfy  Cj  <  Cj ,  RFFAM-rx  and  RFFAM-mll  yield  the  same 
initial  allocation  plan  (8,7,7, 6,5,4,3,2)  and  depot  inventory  (40),  which  can  therefore  be 
declared  optimal.  When  Cj  =  Cj ,  the  plans  are  not  identical,  but  are  of  equal  cost. 
Therefore,  we  can  declare  the  solution  given  by  RFFAM-mll  optimal.  When  Cj/cj  >  2  , 

however,  the  solutions  of  the  two  models  do  not  have  the  same  cost.  RFFAM-rx  yields  a 
period-I  allocation  with  all  ships  loaded  to  their  capacities  (64  missiles  in  all),  and  19 
missiles  are  required  at  the  depot.  The  initial  allocation  prescribed  by  RFFAM-mll  is 
v^  =  (8, 8, 8, 8, 8, 8,6,3) ,  with  26  missiles  required  at  the  depot.  Although  both  solutions 
require  the  same  number  of  missiles,  the  optimality  gap  is  about  7%.  We  already  know 
that  =  19 ,  so  the  solution  obtained  from  RFFAM-rx  is  feasible,  but  in  general  that 
may  not  be  the  case. 

It  is  possible  that  investing  further  effort  can  lead  to  dependable  and  fast  solutions 
of  CCNIM-MAP  based  on  RFFAM,  for  problems  with  eight  ships  and  eight  scenarios  in 
each  period.  But  we  require  optimal  solutions  for  significantly  larger  cases:  this  prompts 
the  development  of  a  specialized  algorithm  in  the  next  chapter.  The  specialized 
algorithm  can  identify  x^,  x^  ,  or  some  xg  x(f),  depending  on  the  cost  ratio,  without 
requiring  an  LP-based  solution  at  all. 
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V.  SPECIALIZED  ALGORITHMS  EOR  CCNIM-MAP 


This  chapter  develops  a  specialized  enumerative  algorithm  that  calculates  the 
values  of  the  points  discussed  in  the  previous  chapter,  and  therefore  solves  CCNIM- 
MAP  for  a  wide  range  of  cost  ratios.  The  specialized  algorithm  calculates  the  value  of 
either  or  if  and  at  least  one  point  from  x(f)  if  q  =C2.  It  does  so  by 

decomposing  the  two-period  problem  into  single-period  problems  and  using  the 
remainders  following  period  I  as  data  for  the  period- II  problem.  However,  as  discussed 
in  the  previous  chapter,  the  optimal  period-I  allocation  and  the  resulting  optimal  depot 
inventory  need  not  be  an  optimal  solution  of  CCNIM.  Instead,  the  algorithm  we  propose 
enumerates  a  set  of  candidate  period-I  allocations,  calculates  the  depot  inventory 
requirement  resulting  from  each,  and  identifies  the  minimum-cost  solution  from  these 
candidates. 

A.  DECOMPOSITION 

The  period-I  problem  comprises  constraints  (3.13),  (3.15),  and  (3.16)  from 
CCNIM-MAP.  To  simplify  constraints  (3.15),  we  replace  the  demand  vector  in  (3.13) 

with  the  scenario  requirements  vectors  dd  =  max|d\y| ,  which  represent  the  minimum 

allocations  required  in  satisfied  scenarios.  We  obtain  the  following  formulation  of  the 


period-I  problem,  CCNIM-pI. 

(CCNIM-pI)  minqLv’  (5.1) 

s.t. 

Pr|v‘>dd’|>p'  (5.2) 

v'<v  (5.3) 

v'  e  (5.4) 


As  discussed  in  Chapter  IV,  an  optimal  solution  of  CCNIM-MAP  must  specify  a 
period-I  allocation  that  is  feasible  in  CCNIM-pI.  CCNIM-pI’ s  feasible  region  is 
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characterized  by  the  PEPs  of  the  distribution  of  dd  ,  v’  ,  as  seen  by  reformulating 
constraint  (5.2)  as  v'  g  U(  v'  +  Z" )  (see  Chapter  I  for  details). 

j^j 


For  any  given  period-I  allocation  ,  we  derive  5*  period- II  problems,  each  one 


being  a  PIP.  In  each,  we  calculate  ,  the  depot  requirement  for  each  combination  of 

period-I  allocation  and  scenario.  Because  we  require  success  in  period  II  for  every 
period-I  scenario,  the  depot  inventory  requirement  induced  by  each  allocation  plan  is 

v”  (v' )  =  maxv”.  The  remainders  vector,  =  (v'  -d')  ,  is  used  as  data  in  this  model, 

which  we  call  CCNIM-pII.  (“pll”  stand  for  “period  II”). 

(CCNIM-pII)  minv“  (5.5) 

"5.4' 


s.  t. 

Pr{v^^Cp(v:;))d"}>pf 

l^w,-v^=0 

w,  e  z: 

where  v"=r^.  +  w,. 

and  arg|(y  ‘v)  -(y  ‘v)  >0  V/<n-l| 


(5.6) 

(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.11) 


Because  the  vectors  f„.  need  not  be  monotonic,  we  may  wish  to  allocate  missiles 
to  the  ships  in  some  order  other  than  lexicographically,  and  assign  targets  accordingly. 
However,  because  of  the  MAP,  the  demand  in  the  component  of  d”  is  always  assigned 
to  the  same  ship.  Therefore,  we  may  determine  the  mission  assignments  simultaneously 
with  the  allocations  themselves.  The  assignments  are  determined  by  any  permutation 
matrix  that  orders  the  chosen  allocations  monotonically  to  match  large  demands  with 
large  inventories.  In  the  formulation  of  CCNIM-pII,  we  apply  the  inverse  permutation, 
denoted  yMAp(’'')>  to  the  demand  vector.  For  example,  consider  the  very  simple  case 
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where  the  two  demand  vectors  are  (4,4)  and  (6,3),  and  only  one  must  be  satisfied.  If  the 
remaining  inventory  is  (3,6),  we  can  satisfy  the  second  scenario  without  allocating  any 
more  missiles  to  the  ships. 

As  in  the  period-I  problem,  the  period-II  missile  requirements  take  the  inventory 
lower  bounds  into  account.  Because  ships  will  not  return  to  port  if  they  do  not  need  to 
load  more  missiles  to  satisfy  their  projected  demands  in  the  scenarios  we  plan  to  satisfy, 
the  number  of  missiles  required  by  each  ship  from  the  depot  is  the  positive  part  of  the 
difference  between  the  demands  (adjusted  for  lower  bounds)  and  the  remainders.  If  the 
difference  is  negative,  then  the  excess  missiles  are  “wasted”  due  to  the  no-transshipment 
policy  (although  they  may  be  used  by  the  ship  in  scenarios  we  are  not  planning  to 
satisfy).  Hence,  the  number  of  missiles  each  ship  requires  from  the  depot  to  satisfy 
period-II  scenario  s' ,  following  period-I  scenario  5 ,  is  given  by 

dd“j^.  =  (max{yLL(’v“)d“^,v|-r^^.j  ,  if  period-I  allocation  j  is  chosen.  We  let  ddj; 

denote  the  n-dimensional,  integer,  random  vector  of  demands.  Obviously,  to  determine 
the  requirement  prior  to  solving  CCNIM-pII,  we  must  be  able  to  determine  the  optimum 
permutation  matrix  )  before  choosing  v“. 

As  a  first  step  in  solving  CCNIM-pII,  we  claim  that  if  ship  capacities  are  equal, 
constraint  (5.8)  can  be  dropped,  and  the  optimal  permutation  of  demands  lists  them  in  the 
same  order  as  the  remainders.  For  notational  convenience,  we  may  view  this  as  a 
permutation  of  the  remainders  in  non-increasing  order.  If  ships  have  unequal  capacities, 
a  more  elaborate  procedure  would  be  required  to  ensure  that  the  specified  allocations  can 
actually  be  loaded  onto  the  ships,  but  we  ignore  this  complication  in  this  dissertation,  and 
assume  that  v.=v\/i.  The  validity  of  the  sorting  operation  on  the  remainders  is 
provided  by  the  following  Lemma  5,  below. 

Let  r  e  Z"  represent  the  remainders  following  an  arbitrary  period-I  scenario. 
There  are  n  \  permutations  of  the  elements  of  r .  Let  K  =  {l,...,n!}  be  the  index  set  of  all 
permutations,  and  let  ,  k  e  K ,  denote  the  k*  permutation  of  r .  We  denote  the 
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monotonically  ordered  remainders  vector  by  r* ,  and  let  k*  denote  some  corresponding 
permutation  index. 


Lemma  5. 

Let  D  =  Uid.}  be  a  set  of  scenarios,  with  corresponding  index  set  S ,  where  each 

s 

column  s  is  ordered  in  non-increasing  order.  Let  g  K  be  a  permutation  of  the 
remainders,  and  let  k*  be  a  permutation  index  such  that  the  components  of  r*  are  non¬ 
increasing.  Letting  be  the  required  depot  inventory  if  scenario  5 

i 

occurs,  for  remainders  permutation  ,  we  have  =  nun  . 

Proof: 


Let  k  G  argmin  ,  and  suppose  that  r’'  ^r* .  Then,  there  exist  indices  I  and  m 


A:eK 


such  that  I  <m  and  r/  <  .  Consider  ,  an  alternative  permutation  of  r  ,  in  which 


(5.12) 

(5.13) 

(5.14) 


We  obtain 


■^sk'  -^sk 


(5.15) 


Because  we  assume  and  ry  <  ,  there  are  six  cases  to  examine: 

If  <  r/  ,  then  equation  (5.15)  reduces  to  [O  -I-  O]  -  [O  -I-  O]  =  0 . 


•  If  >  rj; ,  we  have 


ms  m 


[dis  -rt)  +  -  r/ )]  -  -  r/ )  +  [d^  -  )]  =  0 . 


•  If  d,^>r^>d^^>r,\  we  have  [(d,, )  +  (d„, -r/ )]-[(4 -r/ )  +  0 


d  -C<0. 

ms  m 


•  If 


d,s  ^  ^  ^  d^^ ,  we  have  [(d,,  -  )  +  o]  - [(d,,  -  r/ )  +  0 


=  <-<<0. 
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•  4  >  d,„ ,  we  have  [0  +  0]  -  [(d„  -  ^‘ )  +  O]  =  -  d„  <  0 . 

•  If  '■„*  i,..  ^  '■,* .  »e  have  [o  +  (d„  -  4 )]  -[(d„  -  ) +o]  =  d.,  -  d„  <  0 . 

We  conclude  that  <0  in  every  case.  If  ^r*,  set  k  =  k'  and  find  a 

new  permutation  that  satisfies  conditions  (5. 12)-(5. 14).  Because  |k|  is  finite  and  each 
k'  indexes  a  different  permutation  of  r ,  the  sequence  of  permutation  indices  k,k',...,k^'^^ 
must  reach  r*  without  increasing  the  value  of  xf^. .  | 

It  follows  that  the  permutation  r*  minimizes  the  depot  inventory  required  to 
satisfy  any  given  set  of  scenarios  D ,  and  in  particular  any  -feasible  subset.  Therefore, 
assume  that  the  remainder  vectors  are  listed  in  non-increasing  order.  We  arrive  at  this 
simplified  formulation  CCNIM-pIIm  (“pllm”  stands  for  “period  II  modified”): 


(CCNIM-pIIm)  minv“ 

^ sj 

(5.16) 

s.  t. 

Pr{w,.  >dd“}>pf 

(5.17) 

l^w,-<=0 

(5.18) 

(5.19) 

W.,  e  z: 

(5.20) 

where  dd"  =  (max  {d" ,  y}  -  Cp  )  Ij )" 

(5.21) 

and  Cp(r)Garg{(Tr).-(Tr).^j>0  V/<n-l} 

(5.22) 

YeY 

CCNIM-pIIm  is  a  simple  PIP,  whose  solution  is  given  by  finding  the  PEP  w*^. , 

such  that  is  minimal.  There  are  two  simple  ways  to  find  using  the 

“forward  enumeration”  described  by  Beraldi  and  Ruszczynski  [2001]  (hereafter  referred 
to  as  PEP  enumeration),  and  the  MSP  algorithm  (hereafter  referred  to  as  MSP  A), 
described  in  Kress  et  al.  [2004]. 
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The  forward  enumeration  begins  by  obtaining  a  vector,  which  bounds  the  PEPs 
using  the  marginal  distributions  of  the  requirements  from  below.  This  lower  bound  is  set 
as  the  root  node  in  a  PEP  search  tree.  Prom  each  node  in  the  tree,  we  generate  the  next 
level  of  nodes  by  generating  all  candidate  allocation  vectors  that  have  one  additional 
missile  in  one  component.  If  a  candidate  allocation  satisfies  the  demand  vectors  with  the 
required  probability,  it  is  a  PEP,  and  no  more  levels  are  generated  from  it. 

MSPA  begins  by  ordering  all  scenarios  by  their  aggregate  demand  in  non¬ 
decreasing  order.  An  upper  bound  on  the  value  of  is  obtained  from  the  allocation 

required  to  satisfy  the  p-feasible  subset  defined  by  the  first  scenarios  (in  the  aggregate 
ordering)  whose  aggregate  probability  of  occurrence  satisfies  p.  Any  scenario  which 
requires  more  missiles  to  be  satisfied  cannot  be  part  of  the  optimal  p-feasible  subset,  and 
is  eliminated.  If  the  aggregate  probability  of  the  remaining  scenarios  exceeds  the 
threshold,  MSPA  generates,  based  on  single-ship  criteria,  a  list  of  scenarios  that  might 
not  be  included  in  an  optimal  p-feasible  subset.  MSPA  then  enumerates  the  p-feasible 
subsets  resulting  from  elimination  of  appropriately  weighted  sets  of  scenarios  from  that 
list.  Because  MSPA  typically  enumerates  only  a  small  fraction  of  all  few  p-feasible 
subsets,  it  may  be  viewed  as  an  advanced  version  of  CCNIM-e. 

Clearly,  CCNIM-pIIm  can  be  formulated  as  an  IP  and  solved  by  an  LP-based 
branch- and-bound  solver.  KPP  indicate  that,  for  a  very  large  number  of  scenarios, 
branch- and-bound  is  more  efficient  than  MSPA,  but  as  we  show  later  in  this  chapter, 
under  those  circumstances,  PEP  enumeration  is  also  much  more  efficient,  so  it  is  unlikely 
that  an  IP  would  be  needed.  If  the  problem  becomes  very  large,  then  implementing  some 
of  the  specialized  decomposition  procedures  described  in  Chapter  I  may  be  worthwhile. 
This  dissertation  focuses  on  the  decomposition  properties  of  CCNIM-MAP,  which  are 
independent  of  the  solution  method  for  CCNIM-pIIm.  We  believe  that  the  simpler 
algorithms  are  more  efficient  for  problems  of  practical  size,  and  so  restrict  our 
computational  analyses  to  those  algorithms.  Another  advantage  of  using  MSPA  or  PEP 
enumeration  is  that  we  can  solve  CCNIM-MAP  without  requiring  specialized  solver  like 
CPLEX. 
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B.  ALGORITHM 

Because  of  the  differences  between  CCNIM-MAP  and  CCGIM,  and  because  we 
seek  solutions  for  a  greater  range  of  cost  ratios,  a  simple  decomposition  procedure,  as 
used  by  KPP,  does  not  guarantee  we  obtain  an  optimal  solution.  As  shown  in  section  A 
of  Chapter  IV,  the  point  exposed  by  following  their  decomposition  procedure,  ,  may 

not  be  optimal  for  some  Cj  <  Cj ,  and  is  certainly  not  optimal  when  >  Cj  (except  for  the 

degenerate  case  where  there  is  only  one  feasible  solution.)  We  therefore  need  a  more 
extensive  list  of  period-I  candidate  allocations. 

We  refer  to  the  algorithm  that  solves  CCNIM-MAP  by  CCNIM-dc  (“dc”  stands 
for  “decomposition”).  CCNIM-dc  has  two  distinct  parts.  The  first  part  optimally  solves 
CCNIM  in  the  case  that  Cj  >  Cj .  If  Cj  <  Cj ,  we  invoke  the  second  part  of  CCNIM-dc  to 

calculate  ,  the  optimal  solution  for  cases  where  Cj  <sc  Cj .  This  solution  may  not  be 
optimal,  but  by  setting  and  invoking  the  first  part  of  CCNIM-dc,  we  can  verify  its 

optimality  or  obtain  a  lower  bound  on  the  optimal  solution  value.  If  Cj  ^  Cj ,  we  obtain 
every  feasible  period-I  allocation  plan  that  corresponds  to  the  identified  solutions. 

The  difference  between  the  two  parts  of  CCNIM-dc  lies  solely  in  the  way  in 
which  period-I  candidate  allocations  are  generated.  Theorem  2  proves  that  if  c^>c^, 
some  period-I  allocation  which  is  a  PEP  of  the  distribution  of  the  scenario  requirements 
leads  to  an  optimal  solution  of  CCNIM.  Furthermore,  if  Cy>  then  every  optimal 

allocation  must  be  a  PEP  of  that  distribution.  Therefore,  we  need  only  consider  PEPs  of 
the  distribution  of  the  scenario  requirements  as  potential  period-I  allocations.  The  set  of 
PEPs  is  a  small  subset  of  all  potentially  legal  allocations,  so  this  enumeration  scheme  is 
relatively  efficient. 


Theorem  2. 

Let  Pj  =|v|,...,V|p||  be  the  set  of  all  PEPs  on  the  distribution  of  dd  ,  and  let  P  be 
the  set  of  period-I  allocations  feasible  in  CCNIM-MAP;  Pj  c  P  by  definition.  Let 
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^(v)  =  Cj1^v  +  C2a:“(v)  be  the  total  cost  of  a  procurement  plan  as  a  function  of  the 
period-I  allocation,  and  let  v'e  P  be  an  arbitrary  feasible  allocation  plan.  Then, 


(1) 

If  Cj  >  C2  , 

minz(v)  = 

minz(v') . 

VGPi  ^  ' 

(2) 

If  Cj  >  C2  , 

minz(v)< 

min  z(v') 

VGPi  ^  ' 

V€P-E^  ^ 

Proof: 

Let  v'e  argmin  z(v) .  If  v'e  Pj  then  result  (1)  is  trivial.  Otherwise,  there  must  be 

VGP 

some  v’  G  Pj  such  that  v'  >  v' ;  denote  s'  =  g  s’  I  v’.  >  d’ |  as  the  index  set  of  scenarios 
satisfied  by  v’  . 

For  every  5  g  s'  ,  the  vector  of  remainders  maintains  L  ( ^  ~  ' 

Because  v'>v’  ,  it  is  also  true  that  r^(v')  =  v'-d' ,  and  r^(v')-L(v2)  =  v'- >  0 . 
Thus,  every  missile  above  the  level  v’  is  carried  over  into  period  II  for  every  scenario 
s  G  s' ,  and  could  have  been  placed  in  the  depot  inventory,  rather  than  on  a  ship.  This 
would  have  reduced  the  total  cost  by  -Cj  per  missile. 

Let  s]  denote  the  complement  of  the  set  s' ,  and  suppose  that  for  some  scenario 

/g  s],  v.j-dl  <0  for  some  ship  i.  In  that  case,  if  v\>v].,  then  up  to  d\-v]j  more 

missiles  would  be  expended  in  period  I  and  not  be  carried  over  into  period  II.  These 
“wasted”  missiles  would  still  have  to  be  allocated  out  of  the  depot  in  period  II,  if  scenario 
s'  were  to  occur. 

We  conclude  that  (v')  >  v" (v’ )  “  X  )  ’  which  leads  to: 

i 

z{V)-z[y))  =  c,Y,[v'i-v]j)  +  c^[x\V)-x^\y))) 

(5-23) 
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By  assumption,  z(v')<z(v'),  which,  combined  with  (5.23),  proves  claim  (1). 
Claim  (2)  follows  from  the  fact  that  if  v'  g  P  -  Pj ,  then  ^  -  v]j )  >  0 ,  and  so  if  Cj  >  Cj , 

i 

inequality  (5.23)  is  strict.  | 

We  enumerate  the  entire  set  of  PEPs  using  the  forward  enumeration  scheme 
detailed  in  Beraldi  and  Ruszczyhski  [2002],  and  calculate  the  objective  function 

Zj  =Cjl^v’ +C2v"(v' )  associated  with  each  one  (recall  that  )  denotes  the  depot 

inventory  resulting  from  period-I  allocation  v' ).  The  optimal  solution  is  given  by 

obtaining  g  argmin  Zj ,  and  retrieving  the  appropriate  values  for  the  period-I 

j 

allocation  and  depot  inventory. 

Note  that  an  alternative  approach  would  enumerate  all  -feasible  subsets  of  the 
period-I  scenarios  and  calculate  the  minimal  allocations  to  satisfy  those  subsets.  Perhaps 
some  candidate  subsets  could  be  eliminated,  but  we  have  not  explored  this  possibility. 
This  scheme  may  be  more  efficient  if  the  demand  values  are  themselves  large,  if  the 
number  of  p-feasible  subsets  is  small,  or  if  the  probability  threshold  is  high.  We  do  not 
expect  to  encounter  these  conditions  in  most  practical  applications  of  CCNIM-MAP, 
because  the  demands  are  usually  bounded,  and  there  are  more  than  a  handful  of  scenarios 
in  each  period. 

If  Cj  <  Cj ,  then,  as  discussed  in  Chapter  IV,  we  seek  the  point  ,  which 
minimizes  the  depot  inventory.  To  find  ,  we  initialize  the  set  of  potentially  optimal 
allocations  A  =  |vj|,  where  v|  =  v,  and  calculate  x“=x"(v[)  by  solving  the 
appropriate  instance  of  CCNIM-pIIm.  We  then  begin  an  iterative  trial-and-error  method 
to  find  the  minimum  period-I  allocation  that  yields  ,  thus  obtaining  x^ .  From  the 
allocation  in  A ,  we  create  a  list  of  predecessors,  which  are  the  allocations  that  have  one 
missile  less  in  exactly  one  ship’s  inventory.  Because  the  demands  in  every  scenario  are 
ordered,  we  only  need  to  generate  predecessors  that  maintain  the  ordering  u.  >  u.  for 
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i  <  i' .  Therefore,  from  v|  we  generate  only  one  predecessor,  in  which  Vj,  =  v  for  z  <  n 
and  Vj„=v-1.  We  refer  to  predecessors  that  are  p’ -feasible  as  candidates,  and 
eliminate  those  that  are  not. 

We  use  the  following  enumeration  scheme,  adapted  from  the  backwards 
enumeration  scheme  proposed  by  Beraldi  and  Ruszczyhski  [2001].  In  order  to  avoid 
creating  duplicate  predecessors,  this  scheme  relies  on  the  fact  that  u.  >  u.  i  <  i'  in  any 

candidate  allocation.  With  each  potential  allocation  in  A,  we  maintain  b- ,  the  index  of 
the  ship  where  allocation  j  differs  from  its  successor.  For  each  g  A  ,  we  generate  only 
predecessors  that  reduce  the  inventory  to  either  the  ship  indexed  by  5^  -1  or  by  5j . 

To  check  p’-fgasibility  efficiently,  we  first  calculate  p’ ,  the  minimum  allocation 
plan  that  satisfies  every  period-I  scenario.  Clearly,  if  u^  >p’,  then  it  must  be  pC 
feasible.  If  this  test  is  not  passed,  a  rigorous  check  for  feasibility  is  performed.  Note  that 
it  is  possible  that  J  p\  because  some  scenario  may  have,  for  example,  dl  =v.  Vz ,  and 
never  be  satisfied  by  an  optimal  plan. 

We  calculate  the  induced  depot  inventory  level  )  for  all  of  the  candidate 

allocations.  Any  candidate  inducing  a  depot  level  of  missiles  is  an  improvement  over 
the  allocations  of  the  previous  iteration,  and  is  potentially  optimal.  These  potentially 
optimal  allocations  are  used  to  create  the  next  generation  of  candidate  allocations.  The 
algorithm  repeats  until  no  more  legal  allocations  can  be  found,  or  none  of  the  tested 
allocations  yields  a  two-period  solution  requiring  only  missiles  in  the  depot.  The 
incumbent  period-I  allocations  all  require  missiles,  and  are  equally  favorable. 

The  solution  is  optimal  if  Cj  is  sufficiently  large,  but  calculating  the  minimum 
value  of  Cj  at  which  is  optimal  is  difficult.  We  can,  however,  obtain  an  upper  bound 
on  Cj ,  beyond  which  x^  is  guaranteed  optimal,  by  solving  CCNIM-dc  again,  this  time 
setting  Cj=C2,  and  obtaining  We  then  calculate  A  =  -I- -I- ) .  If 
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A  =  0 ,  then  is  optimal  for  any  Cj  <  Cj .  Otherwise,  a  lower  bound  on  the  optimal 
solution  is  given  by 

i  =  c,P;,-l-A)  +  c,(.i,;+l).  (5.24) 

The  threshold  cost  ratio  for  which  is  guaranteed  optimal  is  calculated  such  that 
z  =  Zj,  and  is  given  by  Cj  /cj  =  1  +  A  .  We  define  the  relative  gap  in  CCNIM-MAP  by 

gap  =  {za-z)lz.  (5.25) 

A  skeletal,  pseudo-code  description  of  CCNIM-dc  is  given  here  using  vector 
notation,  with  vectors  always  being  n-dimensional  columns.  We  assume  that  the  ships 
are  listed  in  non-increasing  order  of  their  inventory  lower  bounds,  and  that  their 
capacities  (upper  bounds)  are  equal.  If  the  capacities  are  not  equal,  a  small  modification 
is  required  to  ensure  a  feasible  solution,  and  that  solution  is  not  guaranteed  optimal  for 
CCNIM-MAP.  Appendix  E  provides  a  more  detailed  description  of  CCNIM-dc. 


Data 

ships 

period-I  scenarios. 

1  period- II  scenarios. 

d'  vector  of  demands  for  period  I  scenario  5 
d“  vector  of  demands  for  period  II  scenario  s' 

(Each  demand  vector  is  preordered  in  non-increasing  order.) 
p’  probability  threshold  for  period  I 

p“  probability  threshold  for  period  II,  following  period-I  scenario  s 
(p^  probability  of  period-I  scenario  s 

conditional  probability  of  period- II  scenario  s'  on  period-I  scenario  5 
V  missiles  capacity  (upper  bound,  for  each  ship,  on  the  number  of  missiles 
it  can  carry.)  Assume  v-=v  for  all  i. 
y  discretionary  lower  bounds  on  missile  load-outs 
Assume  bound  vectors  are  non-increasing. 

Cj  cost  of  allocating  a  missile  to  a  ship 

Cj  cost  of  allocating  a  missile  to  the  depot 
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(1) 

Algorithm  CCNIM-dc 

(2) 

begin 

(3) 

if  Cj  >  Cj  then 

//  Part  1 

(4) 

for  5  :=  1  to  |s’  1  do 

(5) 

dd’  :=max{d‘,y| ; 

//  number  of  missiles  required  in  successful  scenarios 

(6) 

end; 

(7) 

GENERATE  A:={vp. 

v|^|| ,  all  PEPs  of  dd’ ;  //  see  theorem  2 

(8) 

for  j  :=  1  to  |a|  do 

//  allocation  index 

(9) 

for  5  :=  1  to  |s'  1  do 

(10) 

//  remaining  inventories  following  scenario  s 

(11) 

:=  SORT  (r^^. , '  descend ;  //by  inventory  size 

(12) 

dd’v;  :=(max{d“ 

,y|-r^^.)  ;  //  period- II  requirement  vector 

(13) 

:=  minjl^wj 

J  w  ^  J 

(14) 

s.t.  w  is  a  -efficient  point  of  dd“^  j ; 

(15) 

end; 

(16) 

//  we  refer  to  lines  (9)-(16)  as  calculating  ( v’ ) 

(17) 

end; 

(18) 

]:=argmin{c,f  \]+C2 

j 

x“j;  v’:=v(; 

1  i  J  J 

(19) 

else 

H  Part  2 

(20) 

C:={v}; 

//  initialize  current  set  of  candidates 

(21) 

6  =  {n} ; 

//  initialize  predecessor  index 

(22) 

CN  :=  0  ; 

//  initialize  set  of  candidates  for  next  generation 

(23) 

Jl. 

//  see  Part  1 

(24) 

while  C^0 

(25) 

for  j  :=  1  to  |c|  do 

(26) 

for  /:=max|5^. - 

l,l|  to  5^.  do 

(27) 

u:=C,; 

(28) 

u.,  :=  u.  -1; 

//  reduce  one  component  only 

(29) 

if  (u  is  -feasible)  and  >  u.  >  m^+i)  then 

(30) 

CN:=CNU{u};  5N:=5nU{/}; 

(31) 

end; 

(32) 

end; 

(33) 

end; 
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(34)  for  j  :=  |cn|  to  1  by  -1  do 

(35)  if  xj  (u’ )  >  Xj  then 

(36)  CN:=CN\{u‘};  5N:=5N\{5N.} ; 

(37)  end; 

(38)  end; 

(39)  ifCN==0then 

(40)  A  :=  C ; 

(41)  end; 

(42)  C:=CN;  5:=5N; 

(43)  end; 

(44)  v’:=v‘eA; 

(45)  end; 

(46)  end; 

C.  COMPUTATIONAL  RESULTS 
1.  Comparison  with  RFFAM 

To  test  the  performance  of  CCNIM-dc,  we  repeat  each  case  reported  in  Chapter 
11.  The  algorithm  is  implemented  in  Matlab  7.1,  and  run  on  a  2.8  GHz  personal  computer 
under  the  Microsoft  XP  operating  system.  Because  the  run  times  are  very  short 
compared  with  the  timekeeping  resolution  (1/64  seconds),  we  measure  the  time  required 
to  solve  1000  replications  of  each  instance,  and  divide  the  result  by  1000.  The  algorithm 
expends  the  same  amount  of  effort  for  any  ratio  c^jc^  <  1,  and  the  same  amount  of  effort 

when  c^jcy  >  1 .  The  results  are  summarized  in  Table  5.  For  cases  with  c^jc^  <  1,  we 
also  report  the  number  of  dominating  period-I  allocations  that  are  actually  examined  by 
the  algorithm.  For  cases  with  C2/cj  >  1 ,  we  also  report  the  number  of  times  the  algorithm 
solved  a  period- II  problem. 

The  ability  of  the  Matlab  code  to  solve  all  of  the  specified  cases  in  a  few 
milliseconds  is  very  encouraging.  We  perform  more  extensive  testing  of  the 
computational  behavior  by  randomly  generating  instances  for  cases  of  various  sizes. 
Solving  CCNIM-dc  when  Cj/q^l  is  a  two-step  procedure,  and  we  investigate  the 

behavior  of  each  step  separately.  In  particular,  we  compare  the  performance  of  the  two 
methods  available  to  solve  CCNIM-pIIm,  namely  PEP  enumeration  (until  the  first  PEP 

is  found)  and  MSPA.  This  analysis  will  enable  us  to  choose  the  appropriate  algorithm 
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based  on  a  problem’s  characteristics.  When  C2/cj  >  1 ,  we  are  only  solving  instances  of 
the  period- II  problem  while  enumerating  period-I  allocations  of  interest.  In  this  case,  the 
algorithm’s  computation  time  is,  therefore,  directly  proportional  to  the  time  required  to 
solve  a  single  period-II  instance  of  that  size. 


cjc,<\ 

A 

(N 

Case  Name 

Number  of 

Period-I  PEPs 

Average  CPU- 

Time  (msec.) 

Number  of 

Period-II 

Problems 

Average  CPU- 

Time  (msec.) 

case  2a 

2 

3.9 

2 

2.2 

case  2b 

2 

3.8 

2 

2.4 

case  2c 

2 

5.8 

2 

4.6 

case  2d 

I 

2.3 

2 

4.1 

case  2f 

2 

6.7 

2 

6.4 

Table  5.  CCNIM-dc  Solution  Results  for  Specified  Cases. 


“Case  Name”  references  the  cases  defined  for  RFFAM;  see  Appendix  C 
for  details.  “Average  CPU-time”  is  the  average  time  over  1000  runs 
required  to  optimally  solve  the  instance,  in  milliseconds.  “Number  of 
Period-I  PEPs”  is  the  number  of  period-I  allocations  for  which  the  two- 
period  solution  must  be  calculated  to  determine  the  optimal  solution  when 
C2/cj  =  1 .  “Number  of  period-II  problems”  specifies  the  number  of  times 

a  period-II  problem  was  solved.  This  includes  the  initial  solution  of  v" . 
Solution  times  are  extremely  short,  even  for  case  2f,  for  which  RFFAM 
cannot  obtain  an  integer  solution  in  one  hour  of  computation. 


2.  Period-I  Solution  Times 

Here,  we  randomly  generate  100  instances  for  each  case  size;  all  demands  for 
missiles  are  uniformly  distributed  integers  on  [0,8].  Allocation  lower  and  upper  bounds 
are  set  at  2  and  8  missiles,  respectively,  for  all  ships.  We  assume  every  scenario  is  equi- 
probable.  We  observe  the  performance  of  the  PEP  enumeration  as  each  of  the  three 
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parameters,  n  ,  |s'| ,  and  varies  separately.  Because  solution  times  may  be  very  short 

compared  with  the  resolution  of  the  timekeeping  mechanism,  we  repeatedly  solve  each 
instance  until  a  discretization  error  of  less  than  1%  is  achieved  (a  minimum  runtime  of 
1.56  seconds),  and  report  the  observed  time  divided  by  the  number  of  times  the 
enumeration  algorithm  solved. 

It  is  important  to  stress  that  the  specific  results  are  strongly  affected  by  the 
distribution  of  the  generated  scenarios.  The  purpose  of  these  trials  is  only  to  observe 
general  trends,  and  verify  that  problems  of  useful  size  can  be  solved  in  a  reasonable 
length  of  time. 

We  first  observe  the  effect  of  the  number  of  ships  on  the  enumeration  time  of  the 
entire  set  of  PEPs  of  the  period-I  requirements  vectors.  The  number  of  period-I  scenarios 

is  set  at  |s'|  =  12 ,  and  the  probability  threshold  is  set  at  =  0.75 .  The  distribution  of  the 

computation  times  is  summarized  in  the  left  graph  of  Figure  3,  which  depicts  the  median, 
the  quartiles  and  the  minimum  and  maximum  observed  times.  The  graph  on  the  right  of 
Figure  3  displays  the  same  statistics  for  the  enumeration  gap,  i.e.,  the  difference  between 
the  maximum  level  of  any  PEP  of  the  requirements  vector  and  the  level  of  the  PEP  lower 
bounding  vector.  (Recall  that  the  “level”  of  an  integer  vector  is  defined  as  the  sum  of  its 
components.) 

There  are  two  of  observations  we  wish  to  emphasize.  First,  note  that  every 
statistic  except  for  the  minimum  time  increases  exponentially  as  the  number  of  ships 
increases.  The  minimum  value  appears  to  increase  only  linearly,  if  we  ignore  the  sample 
for  n  =  24 .  Second,  it  is  clear  that  there  is  great  variability  in  the  computation  times  for 
cases  of  equal  size,  and  that  this  variability  increases  as  n  increases.  For  cases  of  20  ships 
or  more,  the  computation  time  of  the  entire  set  of  PEPs  ranges  over  three  orders  of 
magnitude.  This  variability  makes  it  difficult  to  predict,  based  on  the  case  size  alone, 
how  difficult  solving  a  specific  case  will  be.  Clearly,  the  easiest  cases  involving  28  ships 
require  less  effort  than  the  hardest  cases  involving  as  few  as  8  ships. 
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PERIOD-1  PEPS:  COMPUTATION  TIME  VERSUS  #  OF  SHIPS.  PERIOD-1  PEPS:  ENUMERATION  GAP  VERSUS  #  OF  SHIPS. 


n  n 


Figure  3.  Period-I  PEP  enumeration  Versus  Number  of  Ships. 

The  graph  on  the  left  presents  the  median  (diamonds),  first  and  third 
quartiles  (circles),  and  minimum  and  maximum  values  (crosses)  of  the 
computation  times  of  the  PEPs.  The  computation  times  are  drawn  on  a 
logarithmic  scale,  as  n,  the  number  of  ships,  varies  linearly.  We  set 
|s'|  =  12  and  p*  =  0.75 .  The  graph  on  the  right  presents  the  same  statistics 

for  the  enumeration  gap,  but  on  a  linear  scale.  The  PEP  enumeration  time 
is  strongly  correlated  with  the  enumeration  gap. 


The  Spearman’s  rank  correlation  [e.g.,  Conover  1999,  pg.  314]  between  the 
computation  time  and  the  enumeration  gap,  taken  over  600  observations,  is  0.8127  with  a 
p-value  of  practically  zero.  It  is  no  surprise  that  the  PEP  enumeration  time  is  strongly 
correlated  with  the  enumeration  gap  because  the  enumeration  procedure  generates  a 
search  tree,  whose  depth  equals  that  gap.  Any  node  in  the  tree,  which  does  not  represent 
a  PEP,  contributes  up  to  n  successors,  depending  on  the  number  of  monotonic  vectors, 
which  conform  to  the  capacity  bounds,  that  can  be  generated  by  adding  1  to  any  single 
component  of  the  vector  represented  by  that  node.  Therefore,  if  we  denote  the 
enumeration  gap  by  g,  the  number  of  nodes  in  the  tree  is  O  (n* ) . 


76 


The  Spearman’s  correlation  is  fundamental  to  understanding  the  behavior  of  the 
computation  time  as  a  function  of  different  case  parameters,  as  we  show  in  the  following 
figures.  To  understand  how  changes  in  some  case-size  parameters  affect  the  enumeration 
time,  it  suffices  to  understand  how  they  affect  the  enumeration  gap.  For  some  parameters 
this  may  be  easier  than  for  others,  but  in  any  case,  this  behavior  is,  to  some  extent,  an 
artifact  of  the  scenario-generation  scheme. 

From  Figure  3,  it  seems  that  the  majority  of  enumeration  gaps  increase  linearly 
with  n.  Intuitively,  this  can  be  explained  as  follows.  The  PEP  lower  bound  is 
constructed  from  the  marginal  distribution  of  each  element  in  the  requirements  vector, 
based  on  all  of  the  scenarios.  There  is  a  positive  probability  that  the  value  of  the  PEP 
exceeds  the  lower  bound  for  each  element  (except  those  equaling  the  capacity  bound). 
As  the  number  of  ships  increases,  and  the  requirement  vectors  contain  more  elements,  the 
expected  number  of  PEP  elements  that  exceed  their  lower  bound  increases,  contributing 
to  the  observed  increase  in  the  enumeration  gap. 

We  next  observe  the  behavior  of  the  PEP  enumeration  algorithm  as  the  number  of 
scenarios  varies,  as  shown  in  Eigure  4.  There  are  12  ships  in  each  scenario,  and  the 
probability  threshold  set  at  p’  =  0.75 .  We  observe  a  slow  increase  in  the  enumeration 
gap,  which  is  responsible  for  the  increase  in  computation  time.  Again,  the  reason  for  this 
increase  is  an  artifact  of  the  scenario-generation  mechanism,  and  does  not  necessarily 
reflect  on  problem  size.  Intuitively,  the  PEP  lower-bounding  vector  has  some  asymptotic 
value,  and,  as  the  number  of  scenarios  increases,  the  computed  lower  bound  approaches 
this  asymptotic  vector.  However,  as  the  number  of  scenarios  increases,  the  number  of  p- 
feasible  subsets  increases  exponentially,  and  the  likelihood  of  there  being  more  than  a 
few  PEPs  increases.  Then,  the  probability  that  some  PEP  would  be  of  a  high  level 
increases,  increasing  the  gap.  These  effects  can  be  seen  in  Eigure  5. 
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PERIOD-1  PEPS:  COMPUTATION  TIME  VERSUS  NUMBER  OF  SCENARIOS. 


PERIOD-1  PEPS:  ENUMERATION  GAP  VERSUS  NUMBER  OF  SCENARIOS. 
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Figure  4.  Period-I  PEP  enumeration  Versus  Number  of  Scenarios. 

The  statistics  presented  are  the  same  as  in  Figure  3.  We  vary  the  number 
of  scenarios,  and  set  n  =  12  and  =  0.75 . 
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PERIOD-1  PEPS;  LOWER  BOUND  VERSUS  NUMBER  OF  SCENARIOS. 


PERIOD-1  PEPS:  MAXIMUM  LEVEL  VERSUS  NUMBER  OF  SCENARIOS. 


Figure  5.  PEP  Lower  Bound  Level  Compared  with  Maximum  Level. 

As  the  number  of  scenarios  increases,  the  PEP  lower  bound  level 
converges  to  some  asymptotic  value,  but  the  maximum  level  appears  to 
increase  slowly.  This  explains  the  increase  in  the  enumeration  gap, 
observable  in  Ligure  4. 


Linally,  we  observe  the  performance  of  the  enumeration  procedure  as  the 
probability  threshold  varies.  Lor  12  ships  and  10  scenarios,  we  vary  the  probability 
threshold  in  steps  of  0.1,  so  that  at  each  observed  probability,  the  size  of  a  p-feasible 
subset  changes  by  one.  As  shown  in  Ligure  6,  decreasing  the  probability  threshold  tends 
to  increase  the  enumeration  gap,  and  as  a  consequence,  the  computation  time  as  well. 
Lower  probability  thresholds  tend  to  have  a  larger  gap  because  the  PEP  lower-bound 
level  tends  to  decrease  faster  than  the  PEP  maximum  level.  At  the  low  extreme,  the  PEP 
lower  bounding  vector  is  the  10%  sample  quantile  marginal  distributions,  and  the 
maximum  PEP  is  the  largest  level  of  the  10  vectors.  At  the  other  extreme,  the  lower 
bounding  vector  is  the  100%  quantile,  there  is  only  one  PEP,  comprising  the  maximum 
over  all  10  vectors,  and  the  two  values  coincide. 
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PERIOD-1  PEPS:  COMPUTATION  TIME  VERSUS  PROBABILITY. 
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PERIOD-1  PEPS:  ENUMERATION  GAP  VERSUS  PROBABILITY. 


Figure  6.  Period-I  PEP  Enumeration  Versus  Success  Probability. 

The  statistics  presented  are  the  same  as  in  Eigure  3.  In  each  case,  we  set 
n  =  l2  and  |s’|  =  12  .  The  probability  threshold  varies  between  0.1  and  1. 

At  lower  success  probabilities,  the  enumeration  gap  tends  to  be  larger, 
resulting  in  longer  computation  times. 


In  summary,  the  enumeration  of  Period-I  PEPs  does  not  require  much 
computation  effort,  and  can  be  completed  in  about  one  second  of  computation  for  cases 
of  practical  size.  The  Period-I  enumeration  need  not  be,  therefore,  an  issue  of  concern  in 
practical  applications. 


3.  Period-II  Solution  Times 

We  now  consider  the  computation  time  for  the  depot  inventories,  required  for 
period  II.  In  order  to  solve  CCNIM-MAP,  CCNIM-dc  must  compute  the 

required  depot  inventory  for  a  single  period-I  allocation  times  several  times.  In  this 
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section  we  report  the  times  required  to  compute  ) ,  for  a  randomly  selected  period-I 
PEP  from  those  obtained  when  |s'|  =  12  and  p’=0.75.  Note  that  computing  ) 

itself  requires  C)(|s’|)  computation  time,  because  the  depot  inventory  is  calculated 

separately  for  the  remainders  following  each  period-I  scenario  s.  We  compare  the 
performance  of  PEP  enumeration  with  that  of  MSPA,  as  a  function  of  n  ,  |s“  | ,  and  pf  . 

This  may  help  us  choose  the  appropriate  algorithm  based  on  the  size  of  the  case  we  wish 
to  solve. 


DEPOT  INVENTORY:  COMPUTATION  TIME  VERSUS  #  OF  SHIRS. 
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Eigure  7.  Depot  Inventory  Calculation:  Time  Versus  Number  of  Ships. 
The  statistics  presented  are  the  same  as  in  Eigure  3.  In  each  case,  |s°|  =  12 

and  =  0.75 .  The  graph  on  the  left  presents  the  computation  times 
achieved  by  PEP  enumeration  and  the  graph  on  the  right  presents  the  same 
data  for  computation  by  the  MSPA.  As  expected,  the  MSPA,  whose 
computation  time  is  theoretically  linear  in  n,  is  much  more  efficient  for 
long  demand  vectors  than  is  PEP  enumeration. 
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We  set  |s“|  =  12  and  =0.75,  and  vary  the  number  of  ships.  The  graphs  in 
Figure  7  summarize  the  distribution  of  the  computation  time  of  the  optimal  depot 
inventory  for  a  single  allocation,  when  it  is  calculated  by  PEP  enumeration  (on 

the  left)  or  by  MSPA  (on  the  right).  As  expected,  the  computation  time  of  the  PEP 
enumeration  algorithm  is  exponential  in  n.  In  contrast,  the  MSPA  algorithm’s  behavior  is 
theoretically  linear  in  n,  and  this  is  corroborated  by  the  results  of  Eigure  7. 

It  is  interesting  to  compare  the  PEP  enumeration  to  the  performance  achieved  in 
period  I.  In  period  II,  the  PEP  enumeration  algorithm  halts  once  the  first  PEP,  of 
minimal  level,  is  found.  However,  the  median  computation  time  is  3.3  times  longer  than 
in  period  I,  averaged  over  n,  because  the  requirement  vectors  in  period  II  are  not 
monotonic.  The  gap  between  the  PEP  lower  bound  level  and  the  level  of  the  minimum 
PEP  tends  to  be  actually  larger  than  the  gap  between  the  PEP  lower  bound  level  and  that 
of  the  maximum  PEP  in  period  I.  Eurthermore,  each  node  generates  more  successors,  on 
average,  because  the  monotonicity  restriction  is  relaxed.  A  worst-case  time  of  750 
seconds  was  observed,  about  1000  times  longer  than  the  worst  case  Period-I  PEP 
enumeration  time  for  20  ships. 

We  next  compare  the  behavior  of  the  two  algorithms  as  the  number  of  scenarios 
varies.  We  assume  a  combat  force  of  12  ships,  and  set  the  probability  threshold  at 
p’  =  0.75 .  The  computation  times  for  both  algorithms  increase  exponentially  as  the 
number  of  scenarios  increases.  However,  the  rate  of  increase  for  MSPA  is  greater  than 
that  for  PEP  enumeration.  When  the  number  of  scenarios  is  small  MSPA,  is  more 

efficient;  for  |s"  |  >  24  it  appears  that  PEP  enumeration  is  faster. 

We  end  this  section  by  comparing  the  behavior  of  the  two  algorithms  as  the 
probability  threshold  varies.  Eor  cases  with  12  ships  and  |s”|  =  20,  we  vary  the 

probability  threshold  over  the  range  0.5-0.95  (we  already  observed  in  period-I  that  the 
result  is  trivial  when  p  =  1 ).  The  results  are  summarized  in  Eigure  9  in  similar  fashion  to 


82 


Figure  7.  As  can  be  seen,  there  is  little  difference  in  the  distribution  of  the  computation 
time  between  the  two  algorithms,  but  MSPA  seems  to  hold  a  slight  advantage  for  all 
thresholds  p”  <  0.9 . 

We  believe  that  in  most  applications  to  naval  forces,  the  number  of  ships  will  not 
exceed  20.  The  number  of  hand-generated  scenarios  is  probably  less  than  10,  but  may  be 
much  greater  if  automatic  scenario  generation  is  applied.  From  the  above  analysis,  it 
seems  that  the  MSPA  is  preferable  for  solving  the  period-II  depot  inventory  problem, 
with  the  caveat  that  the  number  of  scenarios  should  not  be  too  large. 


DEPOT  INVENTORY:  COMPUTATION  TIME  VERSUS  NUMBER  OF  SCENARIOS. 
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Figure  8.  Depot  Inventory  Calculation  Time  Versus  Number  of  Scenarios. 
The  format  of  theses  graphs  is  similar  to  that  of  Figure  7.  We  vary  the 
number  of  period-II  scenarios  and  set  n  =  I2  and  p“=0.75.  Both 
algorithms’  computation  time  increases  exponentially,  but  the  rate  of 
increase  for  MSPA  is  greater  than  for  PEP  enumeration.  Although  MSPA 
is  more  efficient  for  a  small  number  of  scenarios,  PEP  enumeration 
becomes  more  efficient  if  the  number  of  scenarios  exceeds  24. 
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DEPOT  INVENTORY:  COMPUTATION  TIME  VERSUS  PROBABILITY  THRESHOLD. 
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Figure  9.  Depot  Inventory  Calculation  Time  Versus  Probability  Threshold. 

The  format  of  theses  graphs  is  similar  to  that  of  Figure  7.  We  vary  the 

probability  threshold  for  period  II.  and  set  n  =  12  and  |s"|  =  20.  There  is 

not  much  difference  between  distribution  of  the  computation  times  of  the 
depot  inventory  by  PEP  enumeration  and  by  MSPA. 

4.  CCNIM-dc  Solution  Times 

We  next  solve  multiple  instances  of  CCNIM-dc,  for  both  cost  ratios,  to  get  a  feel 
for  solution  times  as  case  size  increases.  We  use  the  MSPA  to  obtain  the  depot  inventory 

in  each  case.  We  generate  100  instances  of  each  case  size,  and  maintain  n  =  |s’|  =  |s"| 
and  p’ =  p“=0.75.  Figure  10  summarizes  the  distribution  of  computation  times  for 
c^jc^<\  (left  graph)  and  Cj^lc^>\  (right  graph).  The  solution  times  increase 
exponentially  as  case  size  increases,  which  is  to  be  expected  from  previous  results. 
However,  even  for  cases  of  size  24,  most  instances  solve  in  less  than  ten  seconds,  and  the 
worst-case  observed  is  only  270  seconds.  A  case  of  that  size  is  probably  as  large  as  we 
would  need  to  solve,  and  the  probability  thresholds  are  not  likely  to  be  set  lower. 
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CCNIM-dc:  SOLUTION  TIMES  VERSUS  CASE  SIZE. 
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Figure  10.  Distribution  of  CCNIM-dc  Solution  Times. 

The  graphs  summarize  the  distribution  of  CCNIM-dc  solution  times  as 
the  case  size  varies,  for  any  cost  ratio;  p'  =  p“=0.75.  As  case  size 
increases,  solution  times  increase  exponentially,  and  the  relative 
dispersion  of  the  solution  times  increases  as  well. 


5.  CCNIM-dc  Optimality  Gaps 

CCNIM-dc  solves  optimally  for  most  cost  ratios,  but  if  Cj  >  Cj ,  an  optimality  gap 
(5.25)  usually  remains.  The  size  of  the  optimality  gap  is  data  dependent,  and  instances 
where  it  is  large  can  be  created.  For  any  given  case  parameters,  the  optimality  gap  is 
maximized  when  C2  =  Cj  +  £•  and  ^  0  .  For  the  random  scenario  generation  scheme  we 
adopt  for  these  tests,  however,  the  gap  tends  to  be  small.  The  mean  relative  optimality 
gap  is  roughly  5%,  and  the  maximum,  over  all  cases,  is  18%.  The  equilibrium  cost  ratio, 
however,  tends  to  increase  with  case  size,  and  in  some  instances  is  as  high  as  21.  It 
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appears,  then,  that  in  large  cases,  there  is  a  much  wider  range  of  cost  ratios  for  which 
some  optimality  gap  persists,  but  the  gap  is  not  necessarily  larger  than  with  smaller  cases. 
Table  6  provides  more  detail. 


Case 

Mean 

Max 

Mean  Threshold 

Max  Threshold 

Size 

Maximum  Gap 

Maximum  Gap 

Cost-  Ratio 

Cost-  Ratio 

8 

0.0583 

0.1795 

4.48 

15 

12 

0.0528 

0.1610 

6.13 

20 

16 

0.0449 

0.1234 

6.85 

20 

20 

0.0445 

0.0963 

8.5 

20 

24 

0.0409 

0.0868 

9.32 

21 

Table  6.  CCNIM-dc  Optimality  Gaps  for  Specified  Cases. 


“Case  Size”  references  the  number  of  ships  and  scenarios  in  each  period, 
which  are  all  equal.  The  maximum  optimality  gap  is  given  by  -  z)/z  , 

and  the  threshold  cost-ratio,  above  which  is  guaranteed  optimal,  is 
given  by  Cj  /cj  =  1  -I-  A  .  Notice  that  although  the  threshold  increases 
significantly  with  case  size,  the  optimality  gap  does  not.  In  fact,  both  the 
mean  and  the  maximum  observed  gaps  tends  to  decrease  with  case  size. 

D.  CONCLUSIONS 

This  chapter  has  developed  CCNIM-dc,  an  algorithm  to  directly  solve  CCNIM 
without  the  use  of  an  LP-based  solver.  CCNIM-dc  requires  that  all  ships  have  equal 
capacity,  and  yields  an  optimal  solution  for  a  wide  range  of  cost  ratios.  Although  the 
computation  time  increases  exponentially  with  case  size,  these  times  are  still  quite 
reasonable  even  for  large,  practical  problems.  All  cases  are  solved  optimally  except 
when  1  <  Cj/cj  <  1-1- A .  But,  when  Cj/cj  is  within  this  range,  the  average  maximum 
optimality  gap  is  only  about  5%. 
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VI. 


CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

We  seek  to  determine  minimum-cost  initial  ship  load-outs  plus  depot- level 
inventory,  while  ensuring,  with  sufficiently  high  probability,  that  all  ships  can  satisfy 
their  assigned  missions  for  each  combat  period.  Each  mission  represents  a  target  and  its 
associated  “demand,”  which  is  the  number  of  missiles  required  to  successfully  engage 
that  target.  A  “cost  ratio,”  the  cost  of  a  missile  stored  in  the  depot  for  potential  use  in 
period  II  divided  by  those  initially  allocated  to  the  ships,  reflects  operational  preferences 
rather  than  actual  monetary  costs. 

We  initially  formulate  the  problem  as  a  two-period,  which  we  denote  by  CCNIM 
(Chance-Constrained  Naval  Inventory  Model).  We  assume  that  all  targets  will  be  in 
range  of  each  combat  ship  in  the  fleet,  and  the  force  commander  can  assign  any  target  to 
any  ship;  we  denote  these  assignment  options  by  the  term  “flexible  assignments.”  The 
force  commander  assigns  at  most  one  target  to  each  ship,  and  does  so  based  on  the 
available  missile  inventories  and  the  target  demands.  We  further  assume  that  ships  have 
bounded  capacities,  carry  no  safety  stocks,  are  not  recalled  to  port  to  offload  excess 
inventories  prior  to  period  II,  and  there  is  no  direct  transshipment  of  missiles  between 
ships.  Due  to  the  flexibility  in  target  assignments  inherent  in  naval  combat,  the  natural 
model,  CCNIM,  involves  three  stages,  and  is  computationally  difficult  to  solve. 

We  show  that  by  following  a  monotonic  assignment  policy  (MAP),  under  which 
targets  with  larger  demands  are  assigned  to  ships  with  larger  inventories,  the  commander 
will  satisfy  every  scenario  possible  with  his  current  supply  of  missiles.  And,  because  this 
policy  is  reasonable  in  the  heat  of  battle,  we  assume  the  commander  will  follow  it.  Given 
the  MAP,  we  can  reduce  CCNIM  to  a  two-stage  stochastic  program,  CCNIM-MAP, 
which  is  easier  to  solve. 

We  go  on  to  show  that  the  MAP  optimally  solves  the  Flexible  Minmax  Subset 
Problem  (FMSP).  In  this  generic  combinatorial  problem,  a  subset  of  weighted  columns 
must  be  chosen,  such  that  the  aggregate  weight  exceeds  p,  and  the  sum  of  the  row 
maxima  is  minimized.  We  prove  that  if  we  may  permute  each  column  independently. 
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then  an  optimal  solution  can  be  found  by  ordering  every  column  monotonically,  and 
selecting  among  the  ordered  columns.  We  show  that  a  single  period  of  CCNIM-MAP 
can  be  transformed  to  an  instance  of  FMSP.  Algorithms  for  solving  FMSP  are  already 
known  in  the  literature,  and  include  the  MSP  algorithm  presented  in  Kress  et  al.  [2004] 
and  appropriately  tailored  enumeration  schemes  for  “p-efficient  points”  (PEPS),  e.g., 
Beraldi  and  Ruszczynski  [2002] . 

We  propose  a  practical  algorithm,  CCNIM-dc,  that  solves  CCNIM-MAP  for  a 
wide  range  of  cost  ratios,  assuming  that  the  ships’  missile  capacities  are  identical.  If  the 
cost  of  storing  missiles  at  the  depot  is  only  a  little  higher  than  the  cost  of  allocating 
missiles  to  the  ships,  the  resulting  solution  is  not,  generally,  guaranteed  optimal. 
However,  CCNIM-dc  solves  much  faster  than  the  deterministic  equivalent  integer 
program,  RFFAM,  which  often  cannot  solve  problems  of  practical  size  for  any  cost  ratio. 
CCNIM-dc,  which  we  implement  in  Matlab  version  7,  solves  most  test-problem 
instances  with  as  many  as  24  ships  and  24  scenarios  in  each  period,  in  less  than  one 
minute  of  CPU-time  on  a  Pentium  IV,  2.8  GHz  personal  computer.  If  the  result  is  not 
provably  optimal,  the  optimality  gap  remaining  is  usually  just  a  few  percent. 

A  second  advantage  of  CCNIM-dc  is  that,  except  when  the  cost  ratio  is  1,  it 
provides  all  of  the  initial  inventory  combinations  that  solve  the  problem  optimally. 
Planners  can  then  select  a  preferred  allocation  according  to  other  criteria  not  explicitly 
handled  by  the  algorithm.  Furthermore,  CCNIM-dc,  does  not  require  specialized 
optimization  software,  and  may  be  implemented  in  any  generic  computing  language. 

B.  FUTURE  WORK 

A  potentially  important  procurement  policy  would  minimize  the  total  number  of 
missiles  procured,  and  allocate  as  many  of  them  as  possible  to  the  combat  forces.  This 
policy  can  be  summarized  by  a  cost  ratio  slightly  greater  than  one.  CCNIM-dc  does  not 
solve  optimally  for  such  a  cost  ratio,  and  such  an  extension  may  prove  useful. 

A  key  assumption  of  CCNIM-MAP  is  that  the  number  of  targets  does  not  exceed 
the  number  of  shooters  in  any  scenario.  To  analyze  the  opposite  situation  with  our 
techniques,  the  analyst  must  first  predetermine  which  pairs  of  targets  (or  more)  will 
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comprise  a  single  mission,  so  that  the  total  number  of  missions  dos  not  exceed  the 
number  of  shooters.  We  are  unaware  of  any  optimal  method  for  doing  this  except  brute- 
force  enumeration.  Thus,  extending  CCNIM-MAP  to  handle  this  situation  may  prove  to 
be  an  interesting  challenge. 

Flexible  target  assignments  usually  arise  in  naval  combat,  but  in  most  situations, 
and  in  particular  as  force  sizes  grow,  the  assignments  are  not  fully  flexible  nor  are  they 
necessarily  one-to-one.  Semi-flexible  assignments  pose  a  challenge  to  the  current 
modeling  scheme,  particularly  if  legal  target  assignments  are  not  known  before  combat 
actually  occurs.  Methods  to  approximate  such  cases  should  be  developed,  as  solutions  of 
CCNIM-MAP  may  be  optimistic. 

Another  issue  that  merits  attention  arises  from  the  fact  that  assigning  targets  in  the 
real  world  is  a  stochastic  process.  In  practice,  not  all  of  the  targets  in  a  scenario  are 
detected  at  the  time  an  attack  begins,  and  some  shooters  will  have  to  be  committed  to 
targets  before  the  force  commander  can  see  the  full  extent  of  the  attack.  Consequently, 
targets  may  be  identified  or  at  least  assessed  incorrectly,  causing  the  expenditure  of  too 
many  or  too  few  missiles.  A  special  model  or  sub-model  will  be  required  to  handle  such 
situations  (for  example,  see  Washburn  [2001]).  A  practical  allocation  plan  should  be 
robust  to  the  effects  of  combat  uncertainty,  and  techniques  to  ensure  this  should  be 
explored. 

Finally,  we  note  the  need  for  handling  this  issue:  potential  enemy  interdiction  of 
our  assets.  Interdiction  may  occur  through  (a)  direct  attack  on  the  missile  depots,  or  (b) 
disruption  of  access  to  those  depots  by  submarines  or  offensive  mining  operations. 
Furthermore,  missiles  that  are  placed  on  ships  may  be  lost  in  combat  if  those  ships  come 
under  attack.  In  theory  it  is  not  too  difficult  to  add  scenarios  to  incorporate  potential 
interdictions,  but  actually  solving  such  a  model  might  require  substantial  research  effort. 
In  any  case,  it  seems  that  handling  larger  targets  sets  is  a  prerequisite  to  a  consistent 
assessment  of  the  effect  of  loss  of  shooters. 
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APPENDIX  A  -  FFAM 


For  comparison  purposes,  we  provide  the  formulation  of  FFAM,  the  original  IP 
used  to  solve  CCNIM,  which  appeared  in  Avital  [2004].  This  formulation  follows 
identical  assumptions  to  those  of  RFFAM,  but  has  many  more  binary  variables.  It  uses 
auxiliary  models,  similar  to  those  used  by  RFFAM,  to  generate  some  of  the  required 
parameters. 

A.  FFAM  SPECIFICATION 

The  following  non-mathematical  description  of  FFAM  itemizes  a  list  of  problem 
requirements  followed  by  the  constraint  keys  where  they  are  represented  in  the 
mathematical  model  presented  in  the  next  section. 

Regular  constraints 

•  Minimize  the  weighted  cost  of  allocating  missiles  to  the  combat  ships  in  period  I 
and  storing  extra  missiles  at  a  depot  for  possible  use  in  period  II,  (7.1). 

Subject  to: 

•  Each  scenario  in  each  period  is  successful  only  if  every  ship  has  enough  missiles 
to  satisfy  the  demand  of  its  assigned  mission,  (7.2)  and  (7.17). 

•  In  each  period,  the  probability  of  successfully  covering  the  scenarios  exceeds  a 
user-specified  threshold,  (7.3).  In  the  second  period,  the  cumulative  probability 
must  be  achieved  for  every  possible,  preceding,  period-I  scenario,  (7.18). 

•  Each  ship  is  allocated  a  specific  number  of  missiles  in  period  I,  (7.4). 

•  Each  ship  in  each  scenario  is  assigned  one  mission,  (7.5)  and  (7. 19). 

•  Each  mission  in  each  scenario  is  assigned  to  exactly  one  ship  (with  a  specific 
number  of  missiles),  (7.6)  and  (7.20). 

•  Eollowing  assignment  and  prosecution  of  a  mission,  each  ship  in  each  period-I 
scenario  maintains  an  inventory  equal  to  its  initial  level  less  the  demand 
associated  with  its  assigned  mission,  if  that  demand  can  be  met,  (7.7).  Otherwise, 
the  remaining  inventory  is  zero,  (7.8).  (All  ships  “stay  and  fight.”) 

•  Each  ship’s  post-scenario  inventory  equals  the  number  of  missiles  remaining  in  its 
inventory  after  prosecuting  a  mission,  (7.9)  and  (7.10). 

•  Eor  each  ship,  the  number  of  missiles  that  may  be  placed  on  it  is  bounded  from 
below  (discretionary  operational  constraint)  and  from  above  (physical  capacity 
limit),  (7.11)  and  (7.21). 


95 


•  For  each  period-I  scenario,  each  ship’s  total  inventory  of  missiles,  before 
prosecuting  any  period- II  mission,  equals  the  post-mission  inventory  after  period  I 
plus  any  missiles  that  are  replenished,  (7.12)-(7.14). 

•  For  each  ship  and  each  period-I  scenario,  only  one  level  of  replenishment  may 
take  place  following  period  I,  (7.15). 

•  Following  each  period-I  scenario,  the  total  number  of  missiles  distributed  to  the 
ships  between  periods  of  combat  may  not  exceed  the  number  kept  in  the 
depot,(7.16) . 

Specialized  constraints 

•  In  each  scenario,  missions  with  greater  demands  are  assigned  to  ships  with  greater 
inventories  (MAP  constraints),  (7.22)  and  (7.23). 

•  The  number  of  missiles  allocated  to  ship  i+\  does  not  exceed  the  number 
allocated  to  ship  /,  (symmetry-breaking  constraints)  (7.24). 

•  The  total  number  of  missiles  allocated  in  each  period  must  exceed  some  minimum 
(total-allocation  valid  inequalities),  (7.25)-(7.27). 

•  Each  ship  is  allocated  at  least  some  minimum  number  of  missiles  in  period  I 
(single-ship  valid  inequalities),  (7.28). 


B.  FFAM  MATHEMATICAL  DESCRIPTION 
Indices 

i  e  I  ships 

K  level  (number)  of  missiles  ( =  {0, . . . ,  max  y.}  ) 

niG  M  missions 

5  e  s'  scenario  5  in  period  I 

s'e  s“(5)  scenario  /  in  period  II 

Subset  of  period- II  scenarios  that  may  occur 
following  period-I  scenario  5 


Parameters  [units] 

(p^  probability  that  period-I  scenario  5  occurs 

conditional  probability  that  period- II  scenario  s'  occurs,  given 
period-I  scenario  5  occurs 

p’  probability  threshold  for  period  I  (probability  that  the  realized 
scenario  must  be  satisfied) 
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probability  threshold  for  period  II,  if  scenario  s  occurred  in  period  I 

d  demand  associated  with  mission  m  in  scenario  s  [missiles] 

Cj  unit  cost  of  procuring  a  missile  and  allocating  it  to  a  ship 

[$/missile] 

C2  unit  cost  of  procuring  a  missile  and  storing  it  in  the  depot 

[$/missile] 

y_.  discretionary  lower  bound  on  the  number  of  missiles  that  may  be 

allocated  to  ship  i  [missiles] 

V.  physical  upper  bound  on  the  number  of  missiles  that  may  be 

allocated  to  ship  i  [missiles] 


Auxiliary  Parameters  [units] 

minimum  aggregate  ship  load-out  required  to  satisfy  the  period-I 
scenarios  [missiles] 

minimum  aggregate  ship  load-out  required  to  satisfy  the  period-II 
scenarios  following  scenario  s  in  period  I  [missiles] 

V.  period-I  data- specific  lower  bound  on  the  number  of  missiles  that 

may  be  allocated  to  ship  i  [missiles] 


Decision  Variables  [units] 

1  if  ship  i  has  k  missiles  before  the  first  combat  period,  and  0 
otherwise 

V.  number  of  missiles  initially  allocated  to  ship  i  [missiles] 

number  of  missiles  allocated  initially  to  the  central  depot  [missiles] 
r'  1  if  ship  i  has  k  missiles  remaining  following  mission  m  in  period-I 

IKttlS 

scenario  5,  and  0  otherwise 

I  if  ship  i  has  k  missiles  left  following  period-I  scenario  5,  and  0 
otherwise 

I  if  ship  i  receives  k  missiles  following  period-I  scenario  5,  and  0 
otherwise 

X  _  I  if  ship  i  has  exactly  k'  missiles  left  after  scenario  5,  and  k" 
missiles  are  added  at  the  replenishment  opportunity,  and  0  otherwise 

x^^^  1  if  ship  i  is  replenished  to  level  k  missiles  following  period-I 

scenario  5,  and  0  otherwise 
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1  if  the  demand  vector  in  period-I  scenario  s  is  satisfied  by  the 
allocation  plan,  and  0  otherwise 

1  if  the  demand  vector  in  period-II  scenario  s'  is  satisfied  by  the 
allocation  plan  given  that  scenario  s  occurs  in  period  I,  and  0 
otherwise 

1  if  ship  i  has  k  missiles  and  is  assigned  mission  m  in  period-I 
scenario  5,  and  0  otherwise 

1  if  ship  i  has  k  missiles  and  is  assigned  mission  m  in  period-II 
scenario  s',  following  period-I  scenario  5,  and  0  otherwise 


Formulation 


u,x,y,z  .  ^ 

(7.1) 

s.t. 

Period  I: 

V/, 

m 

(7.2) 

Al 

In; 

(7.3) 

> 

II 

(7.4) 

yi,k,ses' 

m 

(7.5) 

Vm,  5GS' 

i  k 

(7.6) 

Kk^d^,.m  =  fL.s  m,  5gS‘ 

(7.7) 

Z  Kk',m  =  r'tns  yi,k  =  0,m,sG  s' 

(7.8) 

Y.^Ls=riks  Vi,  ^,5  g  s' 

m 

(7.9) 

k 

(7.10) 

xlk=0  Vi,  kl(^k  <  v.^U(k  >v.) 

(7.11) 

All  variables  binary. 

Period  II: 

V  <4^  yi,  k,  k'e  K,  seS^ 

ikks 
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(7.12) 


X  <  w,.,-  Vi,  k,  k'e  K,  se 

ikk's 

(7.13) 

Vi,  t,  iES' 

k  k  \k  +k  =k 

(7.14) 

Yariks  ='^  V/,  5GS' 

k 

(7.15) 

(■  k 

(7.16) 

Y  Y  ^ikm  ^  4u  5  G  s\  5'g  S"(5) 

m  k>d^^' 

(7.17) 

Y  <Ps'isZs'is^p's 

(5) 

(7.18) 

Y^fkl  =  k,  SGS\  s' G  S“(5) 

m 

(7.19) 

Vm,SES',s'ES"W 

(7.20) 

Vz,  1  (^  <  V,.)U(fc  >  "V^),  5G  S‘ 

(7.21) 

All  variables  binary  except  v“  g  IV . 

Specialized  constraints: 

=  0  Vi,  k,mi^i,  SG  s‘ 

(7.22) 

E  E  “.?»  ^  E  E  “  vt ,  m  <  „  - 1,  .  E  s' ,  /  s  S”  (^) 

t  i  k'>k 

(7.23) 

Y^ik'^Y^Uk'  vz<zz-i, 

/:  '>/:  k  '>k 

(7.24) 

> 

— 1 

K~ 

II 

(7.25) 

i 

(7.26) 

ZEk-^’L^K  V.ES' 

(■  k 

(7.27) 

xi=0  Vi,k<v, 

(7.28) 

With  these  restrictions  on  the  data: 

1 

c 

VI 

> 

+ 

Al 

(7.29) 

^^m+M  V/n<zz-l,  V^gS’us". 

(7.30) 
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This  model  contains  some  variables  that  can  be  substituted  out,  namely  , 

and  V,..  Defining  these  variables  in  constraints  (7.5),  (7.9),  (7.14),  and  (7.25) 

generates  “branching  constraints,”  however,  and  branching  on  these  variables  accelerates 
the  branch-and-bound  solution  process  for  the  integer  model;  see  Appleget  and  Wood 
[2000]. 
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APPENDIX  B  -  SINGLE-PERIOD  MIPS 


Following  are  the  formulations  of  the  two  auxiliary  models  used  to  set  some  of 
the  parameters  used  in  RFFAM.  RFFAM-sp  is  a  single-period  model  using  the  same 
assumptions  as  RFFAM.  RFFAM-lb  is  a  relaxation  of  RFFAM-sp.  Both  models  use 
the  same  notation  as  RFFAM,  and  we  provide  it  here  for  the  reader’s  convenience.  The 
equation  numbers  reference  the  first  appearance  of  these  constraints  in  Chapter  II. 


Indices 

i  e  I  ships 

K  level  (number)  of  missiles 

niG  M  missions 

5  e  s'  scenario  5  in  period  I 

s'g  s“(5)  scenario  /  in  period  II 

Subset  of  period- II  scenarios  that  may  occur 
following  period-I  scenario  5 


Parameters  [units] 

d  demand  associated  with  mission  m  in  scenario  s  [missiles] 

(p^  probability  that  period-I  scenario  5  occurs 

conditional  probability  that  period-II  scenario  s'  occurs,  given  that 
period-I  scenario  5  occurs 

probability  threshold  for  period  I  (probability  that  the  realized 
scenario  must  be  satisfied) 

probability  threshold  for  period  II,  if  scenario  s  occurs  in  period  I 

V .  discretionary  lower  bound  on  the  number  of  missiles  that  may  be 

allocated  to  ship  i  [missiles] 

i7.  physical  upper  bound  on  the  number  of  missiles  that  may  be 

allocated  to  ship  i  [missiles] 
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Decision  Variables  [units] 


V*  number  of  missiles  allocated  to  ship  i  in  period  I  [missiles] 

1  if  ship  i  has  k  missiles  before  the  first  combat  period,  and  0 
otherwise 

number  of  missiles  allocated  to  ship  i  in  period  II  following  period- 1 
scenario  5  [missiles] 

x^^^  1  if  ship  i  is  replenished  to  level  k  missiles  following  period-I 

scenario  5,  and  0  otherwise 

2}  1  if  period-I  scenario  5  is  satisfied,  and  0  otherwise 


I  if  ship  i  successfully  covers  its  assigned  mission  in  period-I 
scenario  s  ,  and  0  otherwise 

I  if  period-II  scenario  s'  is  successful  following  period-I  scenario  5, 
and  0  otherwise 

I  if  ship  i  successfully  covers  its  assigned  mission  in  period-II 
scenario  s'  following  period-I  scenario  5,  and  0  otherwise 

I  if  ship  i  is  assigned  mission  m  in  period-I  scenario  5,  and  0 
otherwise 


^  assigned  mission  m  in  period-II  scenario  s'  following 

period-I  scenario  5,  and  0  otherwise 


A.  RFFAM-sp 

minVv’  (2.16) 

u,v,x,z 


s.t. 


> 

II 

(2.17) 

z'  <- 

^IS  ~  — 

/  A 

Z  ^4  -  dm/ims  + 1  m,  5  e  s’ 

(2.18) 

l_l 

bT 

VI 

1— 1  >1 

\/i,  s’ 

(2.19) 

(2.20) 

II 

—1 

V/ 

(2.21) 

4-0 

Vi,  k\[k  <v_.^[j[k>v.) 

(2.29) 
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(2.41) 


ul=l  V/,  5eS‘ 
u-m  =  0  V/,  mj^i,  5 e 


k  '>k  k  '>k 


All  variables  binary  except  v/ . 


(2.42) 

(2.43) 


v"e{vpv,+l,..-,v,}  Vi,  5eS' 

xjl  =  0  Vi,  k\[k  <y^^^\/{k>v^),  5eS' 

m“',  =  1  Vi,  5e  s\  /e  s“(5) . 

All  variables  binary  except  v” . 


(2.50) 

(2.31) 

(2.32) 

(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 

(2.39) 

(2.40) 

(2.51) 
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APPENDIX  C  -  CASE  SPECIEICATIONS 


Following  are  tables  listing  the  case  specifics  for  the  RFFAM  cases  run  in 
Chapter  11.  In  each  case,  ^j=y|s’|V5,  =  ly/ |s“ |  V/ ,  p^=l-y|s'|,  and 

=  l-y|s“|  V5  .  The  capacity  bounds  are  set  at  y_.  =  2  Vz  and  v.=S\/i.  For  each 

case,  the  length  of  the  demand  vectors  equals  the  number  of  ships,  and  we  use  the 
following  indexing  scheme:  s\  5'e  s“  . 


/ 

/ 

/ 

^'1 

^'2 

^4 

^5 

m. 

6 

5 

6 

8 

7 

6 

6 

3 

5 

2 

2 

5 

1 

1 

3 

0 

0 

4 

Table  7.  Parameter  Specifications  for  Case  2a. 


^2 

^3 

^4 

/ 

^5 

/ 

^6 

/ 

Ay 

^8 

m. 

3 

7 

6 

5 

8 

7 

8 

5 

2 

7 

5 

5 

7 

6 

6 

4 

m3 

2 

3 

2 

4 

5 

5 

4 

3 

0 

2 

1 

4 

4 

2 

1 

1 

Table  8.  Parameter  Specifications  for  Case  2b. 


/ 

/ 

/ 

/ 

/ 

^'1 

^'2 

^'3 

^'4 

^5 

^6 

Ay 

^8 

59 

^10 

m. 

6 

5 

8 

6 

8 

7 

8 

5 

7 

7 

m2 

5 

5 

6 

5 

4 

7 

7 

4 

2 

5 

m3 

3 

5 

5 

4 

4 

4 

3 

2 

1 

4 

^4 

1 

3 

3 

3 

2 

2 

2 

2 

0 

2 

m^ 

0 

3 

1 

1 

2 

1 

2 

2 

0 

1 

Table  9.  Parameter  Specifications  for  Case  2c. 
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^2 

53 

^4 

^5 

^6 

/ 

^7 

^8 

/ 

59 

^I'o 

■^12 

mj 

5 

7 

7 

8 

7 

8 

5 

8 

8 

8 

8 

7 

5 

7 

7 

7 

5 

7 

4 

7 

7 

8 

5 

6 

4 

6 

7 

6 

5 

7 

2 

5 

4 

8 

5 

6 

4 

4 

6 

6 

5 

6 

2 

1 

1 

7 

3 

4 

3 

4 

3 

5 

3 

6 

0 

1 

1 

6 

1 

3 

mg 

0 

3 

0 

0 

0 

5 

0 

1 

0 

5 

0 

1 

Table  10.  Parameter  Specifications  for  Case  2d. 


/ 

/ 

/ 

/ 

/ 

^'1 

^'2 

^'3 

^'4 

^5 

^6 

$1 

^8 

59 

^10 

m. 

7 

8 

8 

8 

6 

7 

5 

8 

7 

8 

m2 

4 

7 

7 

7 

5 

7 

4 

6 

7 

7 

m3 

4 

7 

5 

6 

1 

7 

3 

5 

5 

6 

^4 

4 

1 

2 

5 

1 

6 

3 

5 

5 

1 

m^ 

0 

1 

0 

4 

0 

4 

0 

3 

3 

0 

Table  11.  Parameter  Specifications  for  Case  2e. 


/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

^'2 

^'3 

^'4 

^'5 

^6 

Sj 

^8 

59 

^10 

^11 

■^12 

5i3 

^14 

^15 

^16 

mj 

8 

8 

7 

7 

7 

8 

1 

8 

7 

8 

8 

8 

8 

7 

7 

5 

7 

5 

7 

6 

5 

7 

1 

7 

6 

6 

8 

6 

7 

7 

6 

4 

m3 

7 

4 

6 

6 

3 

6 

4 

5 

5 

5 

7 

4 

6 

6 

5 

3 

m4 

6 

3 

5 

5 

3 

6 

4 

4 

5 

4 

7 

3 

4 

6 

4 

1 

mg 

5 

2 

4 

5 

3 

6 

3 

4 

5 

4 

6 

3 

3 

6 

3 

1 

mg 

3 

2 

3 

4 

3 

4 

3 

3 

4 

4 

5 

3 

3 

2 

2 

1 

m^ 

0 

2 

2 

2 

0 

3 

2 

3 

2 

3 

5 

2 

3 

0 

2 

1 

m3 

0 

0 

1 

0 

0 

1 

2 

2 

0 

0 

1 

1 

1 

0 

1 

0 

Table  12.  Parameter  Specifications  for  Case  2f. 
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APPENDIX  D  -  CCNIM-e  ALGORITHM 


In  this  appendix  we  discuss  the  viability  of  a  simple  enumerative  algorithm  to 
solve  CCNIM-MAP  for  cost  ratios  Cj/cj  <  1 .  Such  an  algorithm  can  overcome  the  no¬ 
safety-stock  effect  and  find  a  point  in  x(f)  ,  as  discussed  in  Chapter  IV.  We  refer  to  this 

algorithm  as  CCNIM-e,  where  “e”  stands  for  enumeration.  Note  that  the  validity  of 
CCNIM-e  is  based  on  results  derived  in  Chapter  V. 

We  define  an  atomic  operation  as  the  procedure  required  to  find  an  optimal 
allocation  and  depot  inventory  assuming  that  =  1  and  that  Cj  =  C2 .  CCNIM-e 

consists  of  multiple  repetitions  of  the  atomic  operation.  Let  denote  the  q''^  p^- 

feasible  subset  of  period-I  scenarios,  and  let  Q*  denote  the  index  set  of  all  such  subsets. 
For  each  scenario  s^S^,  there  is  a  (potentially  different)  set  of  period- II  scenarios 
5“  (5)  and  probability  requirement  .  Let  q“  denote  the  index  set  of  the  -feasible 
subsets  of  *S'”(5).  The  atomic  operation  is  performed  once  for  each  combination  of 
period-I  subset  and  period- II  subsets  corresponding  to  the  scenarios  S^.  The  total 
number  of  operations  N  is  given  by 

|q'| 

V  =  Xn|Q"|,  (7.31) 

q=l  ‘J 

which  can  obviously  be  quite  large.  The  costs  of  all  solutions  are  compared,  using  the 
actual  cost  coefficients,  to  find  the  optimal  solution. 

To  assess  the  practicality  of  CCNIM-e,  we  need  to  obtain  an  estimate  of  the 
length  of  time  required  to  solve  the  atomic  operation,  and  calculate  the  number  of  times 
that  operation  must  be  repeated.  The  atomic  operation  is  a  special  case  of  CCNIM- 
MAP,  and  can  be  solved  by  the  following  simple  algorithm.  It  is  described  here  in 
pseudo-code,  given  in  vector  notation,  with  vectors  always  being  n-dimensional  columns. 
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(1)  algorithm  atomic 

(2)  begin 

(3)  v’ :=  maxjdd' j  ;  //  period-I  allocation  that  satisfies  every  scenario 

(4)  for  5:=1  to  |5^|  do 

(5)  r^=(v’-d')  ;  //  remaining  inventories  following  scenario  5 

(6)  r^  :=  SORT  (r^ ,  'descend ') ;  //  by  inventory  size  (see  Lemma  5) 

(7)  w“  =  maxj(dd“^ -r^)  |;  //  number  of  missiles  to  supplement  each 

ship 

(8)  end; 

(9)  :=  max|l^w“j ;  //  depot  inventory 

(10)  end; 


The  complexity  of  the  atomic  operation  can  be  analyzed  as  follows.  The 
operation  in  line  (3)  requires  one  comparison  for  each  element  of  each  dd' ,  and  so  is 

The  operation  in  line  (5)  performs  one  subtraction  and  one  comparison  for 

each  ship,  and  is,  therefore,  0(n) .  Sorting  the  vector  in  line  (6)  can  be  accomplished  in 

O(nlogn)  operations.  The  number  of  operations  required  in  line  (7)  is  (5)1) . 

Because  the  loop  of  lines  (4)-(8)  is  repeated  |s’|  times,  its  complexity  is  given  by 

o|n|s’|  logn  +  max||s“  (5)11  |,  but  in  practical  cases  we  expect  the  second  term  to 

dominate.  Finally,  in  line  (9)  we  also  require  operations.  Therefore,  we  expect 

the  computation  time  of  the  atomic  operation  to  be  governed,  in  practical  cases,  by 
o|n|s’|max||5“  (5)||j .  Note  that  we  do  not  consider  the  computation  of  the  vectors 

dd;  and  dd",  as  part  of  the  atomic  operation,  because  these  can  be  computed  once  for 
each  instance  of  CCNIM-e. 

We  implement  the  atomic  operation  in  MATLAB™  version  6.5,  and  run  it  on  an 
Intel  2.8  GHz  Pentium  IV  personal  computer.  In  order  to  obtain  a  fair  comparison  to  the 
computation  time  required  by  CCNIM-dc,  we  draw  and  solve  random  instances  of 
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CCNIM-MAP  according  to  the  same  demand  distribution  we  use  in  Chapter  V.  All 
demands  for  missiles  are  randomly  generated,  uniformly  distributed  integers  on  [0,8],  and 
every  scenario  is  equi-probable,  with  probabilities  5^  and 

=y|5“(5)|  V/e  5'”  (5)-  We  vary  the  value  of  n,  |5’|,  and  |>S'“(5)|  between  3 

and  30  in  steps  of  size  3,  for  a  total  of  1000  observations.  Because  the  timekeeping 
function  in  Matlab  is  inaccurate  at  measuring  very  short  periods,  we  measure  the  time 
required  to  solve  each  instance  1000  times  in  succession.  We  derive  the  following 
regression  function  for  the  computation  time  of  a  single  atomic  operation,  measured  in 
micro- seconds. 

=-23.13  +  1.45n  +  18.95|5"|  +  1.48|5"|  +  0.1n|5’|-0.08n|5“|  +  0.06n|5’||5“|  (7.32) 

This  model  includes  every  interaction  term  except  for  |5’||5“| ,  and  appears  to  fit  the  data 
best,  achieving  =  0.9987  and  F  =  122,920 . 


The  estimate  the  computation  of  the  enumerative  algorithm,  we  simply  need  to 
multiply  the  computation  time  of  the  atomic  operation  of  the  correct  size  by  the  number 
of  such  operations  required.  When  the  scenarios  are  equi-probable  and 
|5“  (5)!  =  |5“|  \/s^  S\  the  number  of  p-feasible  subsets  from  a  specific  set  is  given  by 

'  |5|  ^ 


|Q|  = 


risip 


(7.33) 


Therefore,  the  number  of  times  the  atomic  operation  must  be  repeated  to  solve 
CCNIM-MAP  is  given  by 


N  = 


5’ 


5"  p" 


(7.34) 


Although  the  atomic  operation  requires  less  than  a  millisecond  to  solve  in 

problems  of  up  to  20  ships  and  20  scenarios  in  each  period,  it  is  evident  that  this 

procedure  can  only  be  used  if  the  probability  threshold  is  set  such  that  every  period-II 

scenario  must  be  satisfied.  Otherwise,  the  number  of  iterations  is  too  large  to  solve  in 

practical  time.  For  example,  in  a  case  involving  10  scenarios  in  each  period,  and 
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probability  threshold  of  0.9,  we  obtain  N  =  10*°  and  =  229.5  fj,  sec  .  The  time  required 

to  solve  that  case  using  the  enumeration  technique  is  on  the  order  of  26  days  of 
computation. 

Obviously  we  can  do  much  better  by  identifying  scenarios  which  are  never 
candidates  to  be  dropped,  or  scenarios  which  may  be  dropped  every  time.  KPPs  MSP 
algorithm,  which  we  integrate  into  CCNIM-dc,  uses  these  principles  to  avoid 
enumerating  most  of  the  period-II  subsets.  Furthermore,  CCNIM-dc  does  not  enumerate 
the  p’ -feasible  subsets  at  all;  rather,  it  uses  theoretical  results  based  on  PEPs  to  minimize 
the  number  of  period-I  candidate  allocations  we  must  explore. 
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APPENDIX  E  -  CCNIM-dc  ALGORITHM 


Following  is  a  pseudo-code  representation  of  CCNIM-dc.  When  describing  a 
particular  function,  the  variables  in  square  brackets  represent  the  output  of  the  function 
and  the  variables  in  parentheses  following  the  function  name  represent  the  input  to  the 
function.  All  vector  notation  refers  to  n-dimensional  vectors  of  ship  properties,  and  all 
operations  are  component-wise.  Functions  are  listed  in  order  of  appearance.  Variables 
are  local. 


Algorithm  [cost.  A,  x“]  =  CCNIM-dc 
DATA 

/e  {!,...,«}  ships 

=|5p...,5|^,|J  period-I scenarios. 

5“  =|5j,...,5|'^„|J  period- II  scenarios. 

dl  demand  i  in  period  I  scenario  5.  dl  >  d’+j  ^  V/  <  n  -1 

d“  demand  i  in  period  II  scenario  s' .  d"  >  d“j  V/  <  n  - 1 

p’  probability  threshold  for  period  I. 

p“  probability  threshold  for  period  II,  following  period-I  scenario  s. 

(p^  probability  of  period-I  scenario  s. 

conditional  probability  of  period- II  scenario  s'  on  period-I  scenario  5. 
vj.  physical  upper  bound  on  missile  capacity  for  ship  i.  v.=v  V/ 

V;  discretionary  lower  bound  on  missile  capacity  for  ship  i. 
y_.  >  V/  <  n-1 

Cj  cost  of  allocating  a  missile  to  a  ship. 

Cj  cost  of  allocating  a  missile  to  the  depot. 

OUTPUT 

cost  cost  of  the  proposed  procurement  plan 
A  set  of  potential  period-I  allocations 
x”  set  of  accompanying  depot  inventories 
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//  set  of  period-I  scenarios 
//  set  of  period- II  scenarios 


begin 

D'-Up:}; 

d":=UK}; 

for  5  :=  I  to  |5'|  do 

dd;  :=  max|d’,v|;  //  scenario  requirements 

end; 

DD‘  :=Uddl; 

P']  =  getPEPbounds (dd\  (p^,  p')  //  bounds  on  minimum  allocation  vector 

if  Cj  >  Cj  then 

A:=  genAll^DD',  (p^  \/s,  p’,  a’,  p’)  //  generate  period-I  allocations 

x“  :=  depofWeight^A,  d\  d",  pf,  y) ;  //  and  resulting  depot  inventory  for  each 
for  j:=l  to  |a|  do 
costj  :=  v'  -I-  CjX" ; 

end; 

j  :=  argmin  I  costal ;  //  possibly  multiple  solutions 

A:=v,eA;  x":=v“ex”; 

j  j 

else 

C  :=  V  ;  //  initialize  candidate  allocations 

5:=n;  1 1  initialize  index  of  successor 

:=  depotWeight(c,  d',  d",  pf ,  y) ;  //  calculate  resulting  depot  inventory 

while 

[CN,  5n]  :=prevLevel(c,  5,  p’,  DD\  (p^,  p^); 
x“  :=depotWeight(cN,  d\  d“,  (p^\  p“) ; 
for  j  :=  |cn|  to  1  by  -1  do 

if  Xj  >  Xj  then 

CN  :=  CN\|u^.| ;  //  delete  allocations  that  increase  depot 

end; 

end; 

if  CN  =  0  then 

A  :=  C ;  //  record  minimum  period-I  allocations 

end; 

C:=  CN; 

end; 

cost  :=  Cj  1^  Aj  -I-  Cj ;  x"  =  ; 

end; 

end; 
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function  [a,  p]  =  getPEPbounds(D,  \/s,  p,  n) 

DESCRIPTION: 

The  function  calculates  the  lower  and  upper  bounds  on  the  PEPs  of  the 
distribution  of  D  .  The  lower  bounds  are  calculated  through  the  marginal  distribution,  and 
the  upper  bounds  from  the  largest  value  possible. 

INPUT: 

D  set  of  scenarios 
(p^  scenario  probability 
p  probability  threshold 

n  number  of  ships 

OUTPUT: 

a  vector  of  PEP  lower  bounds 
P  vector  of  PEP  upper  bounds 

begin 

for  i  :=  1  to  n  do 

let  tt'  be  a  permutation  of  the  scenarios  such  that  > ...  >  ; 

n 

f  :=  argmax  ^  >  p  ;  //  permuted  index  of  threshold-passing  scenario 

f  7r'(s)=t 

a.  :=  d.^, ;  //  PEP  lower  bound  based  on  marginal  distribution 

fdj  :=  max{  ;  //  PEP  upper  bound  -  highest  demand  possible 

end; 

end; 
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function  [a]  =  genAll (dd\  (p^  V5,  p^,  a’,  p’) 

DESCRIPTION: 

This  function  generates  all  the  period  I  allocations  that  need  to  initially  be 
considered.  By  theorem  5.1,  if  Cj  >C2,  the  candidate  allocations  are  the  PEPs  of  the 

distribution  of  dd’ ,  the  allocation  vector  required  for  period  I.  Otherwise,  we  load  the 

ships  to  capacity  so  we  can  calculate  the  minimum  depot. 

INPUT: 

DD*  period-I  scenario  requirements 
(p^  period-I  scenario  probabilities 

p'  period-I  probability  threshold 
a'  ship  inventory  lower  bounds 
P’  ship  capacity  upper  bounds 
OUTPUT: 

A  set  of  period-I  allocation  vectors 

begin; 


C:=[a']; 

//  initialize  set  of  candidate  PEPs 

A:=0; 

//  initialize  set  of  known  PEPs 

5:=1; 

//  initialize  set  of  predecessor  indices 

while 

for  j:=l  to  C  do 

//  check  p-feasibility  of  each  candidate 

u  =  C.; 

mk  :=  checkPE^u,  DD',  (p^ 

V5,  p^) ; 

if  mk  ==  1  then 

A:=  Au{u}; 

//  include  in  set  of  PEPs 

C:=C\{u|; 

//  remove  from  set  of  candidates 

end; 

end; 

if  then 

[c,  5]  :=NextGen_Pl(c,  5, 

A,  p',  n) ;  //  generate  next  set  of  candidates 

end; 

end; 

end; 
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function  [bool]:=checkPF[y,  D,  \/s,  p) 

DESCRIPTION: 

The  function  calculates  the  cumulative  probability  of  scenarios  which  are 
successfully  covered  by  v ,  and  checks  whether  the  probability  threshold  is  exceeded. 
INPUT: 

V  candidate  allocation  vector 

D  set  of  scenarios 

(p^  scenario  probability 
p  probability  threshold 
OUTPUT: 

bool  is  set  to  1  if  v  is  p-feasible 

begin 

bool:=0; 

for  5  :=  1  to  |d|  do 

z,  :=  0 ; 

if  V  >  dj  then 

z,  :=  1 ; 

end; 

end; 

if  Y.^s(Ps^P  then 
bool:=l; 
end; 
end; 
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function  [CN,  5n]  :=NextGen_Pl(c,  6,  A,  p,  n) 

DESCRIPTION: 

This  function  generates  the  next  generation  of  candidate  PEPs.  Eower  bounds  are 
based  on  the  marginal  distribution  of  each  ship’s  demands.  This  enumeration  scheme 
follows  that  described  by  Beraldi  and  Ruszczynski  [2002],  with  additional  conditions  that 
guarantee  the  candidates  are  monotonic  and  the  upper  bound  is  maintained. 

INPUT: 

C  set  of  failed  candidates 

5  their  predecessor  indices 

A  current  known  PEPs 

P  PEP  upper  bounds 

n  number  of  ships 

OUTPUT: 

CN  next  set  of  candidate  PEPs 
6N  their  predecessor  indices 


begin; 

CN  :=  0  ;  //  initialize  set  of  candidates  for  next  generation 

5N  :=  0 ;  //  initialize  set  of  predecessor  indicators 

for  j:=l  to  |c|  do 
for  f  :=  S.  to  n  do 
u:=C,; 

Uj,  :=  Uj,  + 1 ;  //  increase  one  component  only 

if  [uj-  <  u.  V/  <  /-l)  and  [u.  <  then 


CN  :=  CNu{u}; 

5N:=5Nu{/} ; 

end; 

end; 

end; 

for  j:=l  to  |cn|  do 
u:=CN,.; 

mkj  :=  0 ; 

/:=1; 

while  /  <  |a 
v  =  A..; 

if  u  >  V  then 

mkj  :=  1 ; 

/  :=  l-^l  + 1 ; 

else 


//  if  conditions  are  met,  add  to  list  of  next  generation 
//  and  store  its  predecessor  index 


//  initialize  domination  indicator 


//  compare  each  candidate  to  list  of  known  PEPs 
//  record  those  which  are  dominated  by  any  PEP 

//  exit  internal  loop 
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/:=  /+1; 


end; 

end; 

end; 

if  mk  ==  1  then 


CN  :=  CN\{u|; 


5N:=5N\{/}; 

end; 

end; 


//  index  next  known  PEP 


//  delete  the  dominated  candidate  from  the  list 
//  and  its  predecessor  index 
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function  [x“]  =  depotWeight  ( A,  d',  d“,  V5,/,  pf,  v) 

DESCRIPTION: 

This  function  calculates  the  remaining  inventories  on  each  ship,  reorders  the 
remainders,  and  calculates  the  number  of  additional  missiles  each  ship  needs  to 
successfully  cover  each  period-II  scenario.  It  then  calls  minPEP()  to  calculate  the 
minimum  period-II  p-efficient  point  on  the  distribution  of  j .  minPEP()  is  a  generic 

name,  and  actually,  we  call  either  the  function  MSPA()  or  minPEPenum()  to  obtain  the 
minimum  PEP. 

INPUT: 

A  set  of  possible  period-I  allocations 
d’  set  of  period-I  scenarios 
d“  set  of  period-II  scenarios 

conditional  probabilities  of  period-II  scenarios 

pf  probability  threshold  for  period  II,  following  period-I  scenario  5 
y  ship’s  minimum  inventories 
OUTPUT: 

x”  set  of  depot  inventories  calculated  for  each  period-I  allocation 

begin 

for  j:=l  to  |A|  do 
for  5  :=  1  to  |d‘|  do 

r^^.  =  (v'-d')  ;  //  remaining  inventory 

r^^.  :=SORT(r^^., 'descend');  //  see  theorem  5.2 

for  s'  :=  1  to  |d“  I  do 

dd",  j  =  (max |d“ ,  y|  - r^^. )  ;  //  period-II  requirements  from  the  depot 

end; 

dd;:=  U  {dd;,„,,};  //  set  of  period-II  requirements 

s\(p%  >0 

:=  minPEP  (dd“.  ,  >  0,  ) ;  //  optimal  requirements  given  period-I 

//  scenario  5  and  allocation  j 

end; 

:=  max  I  w"  I ;  //  the  depot  required  for  allocation  j 

end; 

x"  =UUh  //  set  of  depot  inventories 

end; 
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function  [A:]  =  minPEPenum(D,  cp^  \/s,  p) 

DESCRIPTION: 

This  function  generates  the  value  of  the  minimum  level  PEP  for  period  II 
requirements  through  the  PEP  enumeration  technique  described  in  Beraldi  and 
Ruszczyhski  [2002]. 

INPUT: 

D  period-I  scenario  demands 
(p^  period-I  scenario  probabilities 
p  period-I  probability  threshold 
OUTPUT: 

k  sum  of  elements  of  the  minimum  PEP 


begin; 

[a,  p]  =  getPEPbounds(D,  (p,,  p) 


k:=fa; 

C:=[a]; 

A:=0; 

6:=1; 

while  A  =  0 
for  j:=l  to  |c|  do 


//  initial  level 

//  initialize  set  of  candidate  PEPs 

//  initialize  set  of  known  PEPs 
//  initialize  set  of  predecessor  indices 

//  check  p-feasibility  of  each  candidate 


mk  :=checkPP(u,  D,  (p^  \/s,  p) ; 

if  mk  ==  1  then 

A  :=  A  u  {u} ;  //  include  in  set  of  PEPs 

end; 

end; 

if  A  =  0  then 

[c,  5]  :=NextGen_PIl(c,  5,  A,  p) ;  //  generate  next  set  of  candidates 

end; 

end; 

end; 
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function  [A:]  =  MSPA(D,  (p^  \/s,  p) 

DESCRIPTION: 

This  function  generates  the  value  of  the  minimum  level  PEP  for  period  II 
requirements  by  using  the  MSP  algorithm  developed  by  Kress  et  al.  [2004].  The  minimal 
p-feasible  subset  is  created  in  a  process  of  elimination.  In  the  preliminary  phase, 
scenarios  that  require  too  many  missiles  in  total  to  be  included  in  the  minimal  p-feasible 
subset  are  eliminated.  In  the  second  stage,  p-feasible  subsets  are  enumerated,  and  the  one 
requiring  a  minimum  number  of  missiles  is  selected. 

INPUT: 

D  period-I  scenario  demands 

(p^  period-I  scenario  probabilities 

p  period-I  probability  threshold 

OUTPUT: 

k  sum  of  elements  of  the  minimum  PEP 


begin; 

for  5  :=  1  to  |d|  do 
end; 

let  be  a  permutation  of  the  scenarios  such  that  > . . .  >  ; 

n 

:=  argmax  ^  >  p  ;  //  permuted  index  of  threshold-passing  scenario 


for  i  :=  1  to  n  do 

5'^:={5l;r^(5)>t^}; 
Pi  ■=  max {4}  ; 

S€.S 

end; 

A::=rp; 

b:=uR}; 

seS'^ 

Pd\B  ^  ’ 

siS” 

'P:=0; 

for  i  :=  1  to  n  do 


//  PEP  lower  bound  based  on  marginal  distribution 

//  nominal  allocation 

//  number  of  missiles  in  nominal  allocation 
//  potentially  optimal  p-feasible  scenarios 

//  probability  of  eliminated  scenarios 

//  initialize  suspect  scenarios 


let  tt'  be  a  permutation  of  the  scenarios  5^  such  that  >...>d] 


t'  :=  argmax  <l-(p  +  p„^J; 

if  t'  0  then 

'P:=Tu{5l;r'(5)<d}; 


120 


end; 

end; 

if  'F  0  then 

^:=min|  !-(/)  + ’I'^lj’  // number  of  scenarios  we  may  eliminate 

for  j  :=\  io  K  do 

let  be  the  set  of  all  subsets  size  j  of  scenarios  from  'F  ; 

end; 

0:=U{0^};  //  set  of  all  possibly  eliminated  subsets 

cm  ^  '  J  ;  //  cardinality  of  O 

j=^  \  j  ) 

for  j:=\  to  cm  do 

if  +  then  //  if  the  subset  can  be  eliminated 

se^  ^ 

k'  :=  V  max  |d.  1 ;  //  calculate  allocation  cost  of  remaining  scenarios 

I 

if  k'  <k  then 

k:=k' ;  1 1  replace  if  improving 

end; 

end; 

end; 

end; 

end; 
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function  [CN,  5n]  :=NextGen_PIl(c,  5,  p,  n) 

DESCRIPTION: 

This  function  generates  the  next  generation  of  candidate  PEPs  in  period  11.  It 
differs  from  the  period-I  generation  in  the  fact  that  there  is  no  monotonicity  condition. 
We  also  need  not  eliminate  dominated  candidates,  because  this  function  is  only  used  until 
the  first  PEP  is  found. 

INPUT: 

C  set  of  failed  candidates 

5  their  predecessor  indices 

P  PEP  upper  bounds 

n  number  of  ships 

OUTPUT: 

CN  next  set  of  candidate  PEPs 
5N  their  predecessor  indices 

begin; 

CN  :=  0  ;  //  initialize  set  of  candidates  for  next  generation 

6N  :=  0 ;  //  initialize  set  of  predecessor  indicators 

for  j:=l  to  |c|  do 
for  f  :=  S,  to  n  do 

u:=C,; 

Uj,  :=  Uj,  + 1 ;  //  increase  one  component  only 

if  [uj,<  then 

CN  :=  CNuju} ;  //  if  conditions  are  met,  add  to  list  of  next  generation 

5N:=5Nu{/};  //  and  store  its  predecessor  index 

end; 

end; 

end; 

if  mk  ==  1  then 

CN  :=  CN\{u} ;  //  delete  the  dominated  candidate  from  the  list 

5N  :=  5N  \  {  /}  ;  //  and  its  predecessor  index 

end; 

end; 
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function  [CN,  5n]  :=prevLevel(c,  5,  p',  DD\  cp^  V5,  ,  n) 

DESCRIPTION: 

This  function  generates  previous  level  of  candidate  period-I  allocations.  This 
enumeration  scheme  follows  the  reverse  PEP  enumeration  described  by  Beraldi  and 
Ruszczyhski  [2002].  Eewer  points  are  enumerated  because  of  the  monotonicity 
requirement.  The  PEP  upper  bounds  are  used  as  a  quick  method  to  verify  p-feasibility. 
INPUT: 

C  set  of  current  known  allocations 

5  their  predecessor  indices 

P’  PEP  upper  bounds 

DD'  Set  of  period-I  required  ship  inventories 

(p^  period-I  scenario  probabilities 

period-I  probability  threshold 
n  number  of  ships 
OUTPUT: 

CN  next  set  of  candidate  allocations 
6N  their  predecessor  indices 

begin; 

CN  :=  0  ;  //  initialize  set  of  candidates  for  next  generation 

6N  :=  0 ;  //  initialize  set  of  predecessor  indicators 

for  j:=\  to  |c|  do 

for  /:=max|5^. -l,l|  to  bj  do 
u:=C.; 

u.,  •=Uj, -I',  1 1  reduce  one  component  only 

if  >M.-  then  //  verify  monotonicity  of  allocation  is  maintained 

mk  :=  0 ; 

if  u  >  p‘  then 

mk  :=  1 ; 
else 

mk  :=checkPp(u,  DD\  (p^  \/s,  ; 

end; 

if  mk  ==  1  then 

CN  :=  CNu{u} ;  //  if  conditions  are  met,  add  to  list  of  candidates 

5N:=5Nu{/};  //  and  store  its  predecessor  index 

end; 

end; 

end; 

end; 

end; 
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